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ABSTRACT 

In this thesis, generahzed Dyson boson-fermion mappings are considered. These are 
techniques used in the analysis of the quantum many-body problem, and are instances 
of so-called boson expansion methods. A generalized Dyson boson-fermion mapping, or 
a Dyson mapping for short, is a one-to-one linear but non-unitary operator that can be 
applied to vectors representing the states of a many-fermion system. A vector representing 
a fermion system maps onto a vector that is most naturally interpreted as representing 
a state of a many-body system that contains both bosons and fermions. The motivation 
for doing such a mapping is the hope that the mapping will reveal some property of the 
system that simplifies its analysis and that was hidden in the original form. The aims of 
this thesis are 

1. to review the theory of generalized Dyson boson-fermion mappings, 

2. by considering a tutorial example, to demonstrate that it is feasible to implement 
the theory and 

3. to find a useful application for a generalized Dyson boson-fermion mapping, by 
considering a non-trivial model, namely the Richardson model for superconductivity. 

The realization of the first two aims mainly involve the collecting together of ideas that 
have already appeared in the literature, into one coherent text. Some subtle points that 
were treated only briefly due to space restrictions in the journal publications where the 
theory was first expounded, are elaborated on in the present work. On the other hand, the 
analysis of the Richardson Hamiltonian that uses a Dyson mapping, goes beyond what has 
already appeared in the literature. It is the first time that a boson expansion technique 
is implemented for a system where the roles of both collective and non-collective fermion 
pairs are important. (The Dyson mapping associates bosons with Cooper pairs, while 
the fermions not bound in Cooper pairs result in fermions being present in the mapped 
system as well.) What is found is that the Dyson mapping uncovers non-trivial properties 
of the system. These properties aid the construction of time-independent perturbation 
expansions for the stationary states of the system, as well as time-dependent expansions 
for transition amplitudes between states. The time-independent expansions agree with 
results that other authors obtained through methods other than boson expansions. The 
time-dependent expansions, that one would be hard-pressed to develop without a Dyson 
mapping, might in future prove useful in understanding aspects of the dynamics of ultra- 
cold fermi gases, when time-dependent magnetic fields are used to vary the atom-atom 
interaction strenght. 
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Molecule, n. The ultimate, indivisible unit of matter. It is dis- 
tinguished from the corpuscle, also the ultimate, indivisible unit 
of matter, by a closer resemblance to the atom, also the ultimate, 
indivisible unit of matter. Three great scientific theories of the 
structure of the universe are the molecular, the corpuscular and 
the atomic. A fourth affirms, with Haeckel, the condensation or 
precipitation of matter from ether - whose existence is proved by 
the condensation or precipitation. The present trend of scientific 
thought is toward the theory of ions. The ion differs from the 
molecule, the corpuscle and the atom in that it is an ion. A fifth 
theory is held by idiots, but it is doubtful if they know any more 
about the matter than the others. 

Ambrose Bierce The Devil's Dictionary (1911) 
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Introduction and Orientation 

This thesis is concerned with a class of hnear transforms that can be apphed to the 
state vectors of some many-body quantum systems. Transforms belonging to this class 
are known as 'generalized Dyson boson- fermion mappings'.^ I make no attempt in the 
introduction to define or explain precisely what such a mapping is. This is done in the 
main body of the text (in Section ITTI) . Here I simply try to indicate the position occupied 
by the content of this text, in relation to other works on the same or similar subjects, and 
in relation to larger fields of study. 

Dyson mappings have been around in various incarnations for roughly forty years. 
These mappings form part of a larger field of study called 'Boson Expansion Methods'. 
The main motivation in the development of boson expansion methods were the need for 
approximations to analyse the nuclear many-body problem. A very complete review of 
this program can be found in |3]. Also, a detailed introductory exposition is contained 
in Chapter 9 of [S]. However, the space dedicated in these texts to mappings of the 
generalized Dyson type is relatively little. To get a better idea of the specific foundation 
of generalized Dyson type mappings, as well as their relation to other branches of the field 
of boson expansion methods, one can consult [Hj. 

The work reported on and extended in this thesis, concerning generalized Dyson boson- 
fermion mappings, were published in the decade after the last of the previously mentioned 
texts appeared. As a result, the theory described in these pages has not yet undergone 
the same 'canonization' as the earlier results. This leaves me with the following tasks: In 
Chapter One the theory of generalized Dyson boson-fermion mappings is presented. The 
argument is developed in a self-contained manner that closely follows that described in 
j7]-PI]. Because the space constraints inherent in journal publications are not an issue 
here, the argument can be given in full, with no or few open questions and unproven 
results. 

Since I have not yet defined what a generalized Dyson boson-fermion mapping is, it 
is hard at this point to explain why these mappings are studied or where they are useful. 
(This problem is only partially resolved after these mappings are properly defined.) The 
course I therefore take is first to develop the theory properly, as if for its own sake, rather 

^This nomenclature is an example of a variant of the Matthew effect In the Gospel According to 
Matthew (Chapter 25, verse 29) there is written: "For unto every one that hath shall be given, and he 
shall hath abundance: but from him that hath not shall be taken away even that which he hath." In 
scientific research, the Matthew effect refers to the tendency to attach an eminent name to every result, 
even when the eminent person to whom the name belongs has little to do with the result. In this case, 
Freeman Dyson cannot be said to have had much to do with the development of the theory expounded 
in these pages. (In a gesture of East- West reconciliation the name 'Maleev' is sometimes appended after 
Dyson's name. Maleev's contribution to the theory expounded in these pages is on par with that of 
Dyson j3j.) 
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than for its usefulness. Afterwards we can then look for properties of the mapping that 
might make it useful in some situations. What we will find is that, unlike other boson 
expansion methods, a general Dyson mapping does in itself not constitute an approxima- 
tion method. Rather, it is a transformation that may cast a system into an equivalent 
form, that is more convenient for certain types of exact or approximate analyses. 

Chapters Two and Three are concerned with applications. In Chapter Two, a nu- 
clear toy model with 5*0(5) symmetry is considered. Here the emphasis is mostly on 
demonstrating that it is feasible to replace the abstract symbols of the general theory 
with concrete numbers, vectors and operators, and that the resulting mapped system can 
be analysed in much the same way as the original. Readers will however probably not 
take from this chapter a sense of what makes a general Dyson mapping a useful analytical 
tool. 

In Chapter Three, I try to remedy this by analysing a highly non-trivial 'realistic' 
model that describes Cooper pairing in a small superconductor. After applying a gen- 
eralized Dyson mapping, I do a perturbative analysis of the mapped system, using ordi- 
nary perturbation theory. (Here I mean the standard techniques of Rayleigh-Schrodinger 
perturbation theory to find the approximate stationary states of the system, and Dirac 
perturbation theory to determine its approximate time-evolution.) It turns out that per- 
turbation expansions are easier to make for the mapped system than for the original. 



Chapter 1 



The Generalized Dyson Mapping: 
Theory 

1.1 Bird's Eye View 

This chapter is concerned with the process of translating the mathematical description of 
a quantum system in a Hilbert space Ti into an equivalent description in another Hilbert 
space 7i. To loosen our limbs and set our terminology, we start by mentioning a few 
things that every reader of this text is probably familiar with already. It is convenient 
to distinguish three components in the description of quantum systems, namely (1) state 
vectors belonging to a Hilbert space 7Y, (2) linear operators that map Ti into itself and 
(3) linear functionals defined on Ti. All information about the system can be derived by 
combining these components to form so-called matrix elements of operators. 

The following concerns the construction of such matrix elements: In 7i, each vector 
\'ip) can be mapped onto a linear functional : Tip ^ C, where, for each in Ti, the 
action of is defined as |0) = (-01 (p). Here (p) is the inner product between 1-0) and 
10) in Ti. As everyone knows, we often write ("01 instead of f^. The mapping 1-0) \—>- {ipl 
is called the canonical map. An arbitrary matrix element of an operator O : Ti ^ Ti 
is a complex number (■01 O The process of mapping should preserve these numbers. 
It is accomplished in a straight forward manner if a one-to-one and onto linear operator 
T : — > can be found. T is called an isomorphism. For the various components, the 
mapping is implemented as follows: 

1. Vectors: Every vector in Ti is mapped onto a vector T in Ti. 

2. Linear Operators: Every linear operator O : Ti Ti is mapped onto a linear 
operator TOT^^ : ii. H.. Note that the mapping O ^ TOT~^ is linear and 
invertible. 
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Section 1.1: Bird's Eye View 
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3. Linear Functionals: Every linear functional f : Ti. ^ C maps onto a linear functional 

The objects T 1-0), TOT^^ and f T^^ are said to be the images of respectively O and 
/ under the action of T. Since O |0) = T'-^) (T0T~^^ (T it is immaterial 
whether we use l-^), |0) and O or their images to calculate matrix elements. The reason 
for insisting on the linearity of T is to ensure a duplication of the linear structure of 
the original system in its mapped counterpart. For instance, if the system is described 
by a superposition alifj) + (3\(/)) then the mapped system is described by an equivalent 
superposition aT \ip) + /3T|0). Operators and functionals also preserve their linearity 
during mapping. 

We will deal with mappings T that are not unitary. On this topic a small aside: It 
has to be remembered that the hermiticity of operators is not preserved by non-unitary 
mappings. Also, consider a vector in 7i that maps onto a vector 1-0) = T in 7i. 
{ipl is the functional obtained from {ip) through the canonical map in Ti, and is the 
functional obtained from li/j) through the canonical map in Ti,} When T is not unitary, 
then in general {ifj\ ^ \ T~^. 

Let us now become more specific. In this chapter we will derive one particular type of 
mapping T, called a generalized Dyson boson-fermion mapping, or a Dyson mapping for 
short. A Dyson mapping T has as its domain the many-fermion Hilbert space Tip (or a 
subspace Vp thereof), that is constructed from the one-particle Hilbert space 7ii,F, where 
Bi^p = {lu)}^^^ is an orthonormal basis for Ti-i^F- With the states that comprise this 
basis are associated fermion creation operators a'^; i/ = 1,2, . . . ,Q, that are conjugated to 
annihilation operators = {a")'^, where + denotes hermitian conjugation in Tip. The 
usual anti-commutation relations between these operators hold. Tip is finite dimensional. 
The vacuum in Tip is denoted |0). A general vector in Tip is denoted by \ifj). The 
functional related to this vector through the canonical map of Tip is denoted {ifj\. The 
subspace Vp oiTip that serves as the domain of the mapping T might also be the whole 
Tip. 

The mapping T maps vectors in the domain Vp, representing states of a many-fermion 
system, onto vectors in the Hilbert space T^bf, representing states of a many-body system 
where some particles are fermions and others are bosons. The Hilbert space H-bf is defined 
as follows: Let Hb he the many-boson Hilbert space constructed from the one-particle 
Hilbert space Hi^s, where Bi^b = is an orthonormal basis for Hi^b- With the 

states that comprise this basis are associated boson creation operators h^; k = 1,2, . . . , M 
that are conjugated to annihilation operators hk = {b^)~^, where we use the same symbol as 

^This notation will be used throughout the text: In Ti., the space from which we map, vectors and 
functionals will be indicated with triangular brackets. In the space Ti into which we map, round brackets 
will be used. 
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before to denote hermitian conjugation, but this time in Ti^. Tlie standard commutation 
relations between tliese operators liold. Ti^ is infinite dimensional. The vacuum in Ti^ 
is denoted |0)^. The infinite dimensional many-body Hilbert space Ti^i? is constructed 
by setting TIbf = 'Hb ® 'Hp- In this space, the operators = ®\ create bosons. 
Conjugated to these are Bk = hk ®\ = {B'^Y which annihilate bosons. Here f denotes 
hermitian conjugation in T-Cbf-^ The usual commutation relations between these operators 
still hold. Similarly a'^ = l^a'^ and a,^ = l^a,^ = {a'^Y respectively create and annihilate 
fermions in T-Cbf with the standard anti-commutators in place. The vacuum in T'Cbf is 
|0) = |0)^ ® |0). A general vector in TIbf is denoted and the functional related to it 
by the canonical map of Ti^i? by (-01. In Ti^F we refer to bosons and fermions as ideal 
bosons and ideal fermions. Note that the ideal boson operators commute with the ideal 
fermion operators. 

In short then, a Dyson mapping T is a non-unitary, one-to-one, linear operator that 
maps vectors representing states of a many-fermion system onto vectors representing states 
of a many-body system, where some of the particles are bosons and some are fermions. 
(Do not read this as a definition. These are just properties Dyson mappings happen to 
have.) Since the mapping T maps vectors in the domain Vp, that is finite dimensional, 
onto vectors in the space TYbf, that is infinite dimensional, the range of T must be a 
finite dimensional subspace Vbf = TiVp) of the infinite dimensional space T~Ibf- Vbf is 
called the physical subspace. 

In the notation of the introductory discussion above, this means that is embedded 
in a larger Hilbert space. For the mapping of operators O : Tip —* Ti-p and functionals 
{ipl : Tip —* C complications then arise. Let us first deal with the functionals. The 
image of {ipl is (i/jI T^^. It is only defined on the physical subspace Vbf- However, the 
most natural way to express linear functionals acting on states in TYsi? of which Vbf is a 
subspace, is as functionals related to vectors in T-Cbp through the canonical map oiTiBP- 
Such functionals are then defined on the whole TIbf- In order to express the mapped 
functionals in this way we therefore have to extend the domain of the image of each 
functional to the whole TiBP- Stated in symbols, for each linear functional (^| on Tip, 
we have to find a linear functional on T-Lbf such that ("^AHy^^ = ("01 T~^. Here the 
suffixed ■\vgp indicates that the object it is preceded by has its domain restricted to the 
physical subspace Vbf- 

A similar problem is encountered for operators. In the Hilbert space TIbf want to 
express all operators of importance in terms of ideal bosonic and fermionic creation and 
annihilation operators. Such operators that are functions of the ideal boson and fermion 

different symbol is used as a precaution to avoid confusion, because we study isomorphisms between 
Vp and subspaces of Ti-BF, that are not unitary. In other words, we work with vector space isomorphisms, 
rather than Hilbert space isomorphisms. 
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operators are defined on the whole T-Cbf- However the image of any operator O : Tip ^ Tip 
is TOT~^, which is only defined on the physical subspace Vbf- Therefore, we have to 
find a linear extension Obf • 'Hbf ~^ T^bfi expressed in terms of ideal boson and fermion 
operators, for the image of every operator O : Hf ^ 'Hf, so that Obf = TOT^^. 

Vbf 

Note that Obf leaves the physical subspace Vbf invariant. 

In this discussion we have already touched upon the answer to the following question: 
Is it not trivial to establish an isomorphism between the many-fermion Hilbert space 
Tip (or a subspace Vp thereof) and a subspace Vbf of the space TLbf"^ After all, we 
could just do the following: If we have a basis Bp = for Hp, (implying Hp is 

Q-dimensional) , we choose any Q-dimensional subspace Vbf of Hbf- For this subspace 
a basis Bvgp = {\ipi)}f^i can be found. A mapping T is then specified by defining 
the image of the basis elements in Bp under T as T \(pi) = What makes the 

Dyson mapping special? A significant motivation for mapping fermion systems onto 
boson-fermion systems is the following. For an operator like the Hamiltonian H of a. 
fermion system, we want to find an operator Hbf that is a function of ideal boson and 
fermion operators and is equivalent to the original fermion Hamiltonian in the sense that 
^bfIvbp ~ THT~^. In other words, the mapping T must be of a type that allows us to 
derive extensions to the whole space Hbf oi the images of the fermion operators, expressed 
as functions of ideal boson and fermion operators. A mapping as arbitrary as the above 
T suggests no clear way of doing this. In contrast, the procedure that is used to construct 
the Dyson mapping allows one to find extensions for the images of all relevant operators 
of the fermion system and explicitly express them as functions of the ideal boson and 
fermion operators. 

It is my aim in this chapter to go through the procedure by means of which a Dyson 
mapping is constructed in detail. However, the danger of such a detailed analysis is that 
the reader looses sight of the forest for the trees. Therefore I try now to give an intuitive 
overview of the argument that comprises the next few sections. 

The argument hinges on the properties of a certain type of superalgebra. For reasons 
that become clear only in hindsight, this type of superalgebra and the properties common 
to certain classes of its representations are introduced in Sections 11.21 and 11.31 In Section 
11.41 we construct a representation of the superalgebra using the fermion operators defined 
in Hp. Then, using the commutation rules inherent to the superalgebra, matrix elements 
of operators of physical interest can be computed. In these computations the details of 
the specific representation of the superalgebra we are dealing with is unimportant. A 
whole class of representations would have given the same answer. This inspires us to try 
and build a representation of the superalgebra in the boson-fermion space Hbf that also 
belongs to this class of representations that all give the same matrix elements. The way we 
go about constructing this boson-fermion representation is as follows. We define a linear 
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operator T that we suspect might be a one-to-one mapping between Tip or a. subspace Vp 
thereof and a subspace Vbf oi the boson- fermion space Hbf- We do not know, however, 
if the operator T is indeed invertible. The operator T has a very convenient property 
though. For operators Op acting in Tip, that are functions of operators in the fermion 
representation of the superalgebra, the nature of the definition of T allows us to find 
operators Obf that are defined on Hbf and are expressed as functions of the ideal boson 
and fermion operators, such that TOp = ObfT. If T is indeed invertible, then Obf is 
the sought-after extension of the image of under the mapping T. We find operators 
Obfi expressed as functions of the ideal boson and fermion operators, corresponding to all 
operators in the fermion representation of the superalgebra. The operators Obf are 
then shown to form a new representation of the same superalgebra, and specifically one 
that belongs to the same class of representations as the fermion representation, for which 
matrix elements of operators come out the same no matter the specific representation. 
This fact is then used to construct a simple proof of the invertibility of T. This all happens 
in Sections 11.51 and 11.61 In short then, the reasons for introducing the superalgebra are, 
firstly that the operator T must be defined in terms of operators belonging to the fermion 
representation of the superalgebra and secondly, the superalgebra is a substantial enough 
structure that, if the operator T takes it as input on the fermion side and reproduces it 
as output on the ideal boson-fermion side, the invertibility of T is guaranteed. 

The reason that I denote this mapping T instead of T is that, for reasons mentioned 
in Section 11.71 the mapping T has some undesirable properties that we want to get rid 
of. Two strategies are proposed, one in Section fl~7l and the other in Section IT!H| that are 
applicable in different circumstances. They both involve acting on T with another linear, 
invertible operator, that for the moment, we will just denote O. The Dyson mapping is 
then given by T = OT. In Section IT^ we discuss why and how the subspace Vbf of 
TiBFi that is the physical subspace, is identified. Finally in Section [1.101 we look back on 
the derivation of the mapping. It is pointed out how this derivation fits into the larger 
literature on the subject. Also, some features of and problems with the mapping are 
mentioned. 

The work that Sections 11.21 - 11.71 are concerned with, first appeared in [Zj - ^D] and 
was reviewed in [TT]. The methods for finding the physical subspace, that are discussed 
in Section Hini are also set out in [12] and jT3]. 

1.2 The Abstract Superalgebra 

A type of superalgebra is defined in this paragraph. For the purposes of this text, a 
superalgebra may be viewed as a vector space, the elements of which are operators, that 
can be partitioned into two subspaces, called the even and odd sectors respectively. The 
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commutator of any two operators belonging to the even sector must be equal to an operator 
that belongs to the even sector. Similarly, the anti-commutator of any two operators in 
the odd sector must be equal to an operator in the even sector. The commutator of an 
operator from the even sector with an operator from the odd sector must be equal to 
an operator from the odd sector. No further knowledge of the theory of superalgebras is 
required for the reader to follow the arguments presented. In fact, the definition given 
here exhausts the author's knowledge of the general properties of superalgebras. A general 
representation is considered of the superalgebra, so that we may be sure that properties 
derived from the definitions are not specific to only one representation of the superalgebra. 
This will be useful when we construct a new representation of the superalgebra from a 
given one, because properties derived for a general representation will then hold for both 
representations. To indicate that we are working with a general representation of the 
superalgebra, over-bars will be superimposed on operators. 

Let 7i be a Hilbert space on which 2M operators Ai and , i,j = 1,2, . . . , M 
are defined, such that operators with upper indices commute amongst themselves, and 
similarly operators with lower indices commute amongst themselves: 



A' , A^ 



Ai , Aj 



0, 
0. 



Set 



Ai , A'^ 



:i.2) 



We assume that the operators Ai, A' and Kl are generators of a representation of a Lie 
algebra, the closure conditions being 



Ki,Ak 



A' 



3 



C^Au 



[1.3) 



where repeated indices, one upper and one lower, imply summation. (When we say that a 
set of operators are generators of a representation of a Lie algebra, we mean the following. 
The representation of the Lie algebra in question is a vector space, the elements of which 
are operators. The set of generators of the representation is just a basis for the vector 
space of operators. Hence, every element of the representation can be expressed as a 
linear combination of the generators.) The Jacobi identity can be used to show that Cjl 
is symmetric in its lower indices and Cjf is symmetric in its upper indices. Using these 
relations it is possible to compute the commutator 



:l4) 



so that we see the assumed commutation relations in eq. 11.31 do indeed ensure that Ai, 
A* and Kj close to a Lie algebra under commutation. Let there further be defined 2Q 
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operators and a'^, fi^iy = 1,2, ... on Ti that obey the anti-commutator relations 



If the commutator relations 



and 



{a^,a^} 


= 0, 




{a^j,,a,y} 


= 0, 




{a^,a^} 


= 


(1.5) 



:i.6) 



hold, then it follows that 



Xi 



:i.7) 



In the derivation of the last two expressions, use has to be made of the Jacobi identity 
and the commutators in eq. 11.71 

To summarize, if we assume that the commutation relations in eq. 11.11 and eq. 11.31 
hold then the commutation relation eq. II. 41 holds automatically. If the anti-commutation 
relations eq. II. 51 and the commutation relations eq. 11.61 and eq. 11.71 are assumed, then the 
commutation relations eq. 11.81 hold automatically. Therefore, under the above assump- 
tions, the operators ^Ai,A^,Kl^ and {ap,^,!} are respectively the generators (read 
basis elements) for the even and odd sectors of a representation of a superalgebra. The 
properties of the superalgebra are determined by the structure coefficients C/^, Cjf , x^^u 
and ■ We therefore indicate the abstract superalgebra as S{C, C, x, X)- 



1.3 The Calculation of Matrix Elements 

Group the generators that were discussed in the previous section into three sets 







-.1 = 1,2,.. 


.,M-,iy = 


1,2, 


= 




-.1 = 1,2,.. 


.,M-,u = 


1,2, 


JC = 




J = l,2,.. 


.,m}. 





Any two elements of either commute or anti-commute. The same holds for the elements 
of OK Assume there is a vector \i) G Ti, called a foundational vector and a linear 
functional (/| on H, called the foundational functional that have the following properties: 
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• 6^ \i) = for all 6^ G 0\ 

• {f\d^ = for all 0^ G 

• Coefficients exist so that Kj \i) = \i) and (/| Kj = (/| for all Kj G /C. 

This definition seems rather arbitrary. 1 therefore mention in passing that it is motivated 
by the concept of the highest weight state of an irreducible representation of a semi-simple 
Lie algebra.^ Beyond this justification of the definition, 1 can only say that it will turn out 
to be a useful definition in the arguments that are to follow. The set of all representation 
of iS(C, C, Xi x) that, for a given set of coefficients possess such a foundational vector 
and a foundational functional is denoted lZs{C, C, x, x^X). A very useful property of the 
representations in 7ls{C, C, x, X, A) can now be stated: 

For any representation in 7ls{C, C, x, Xj Ist be any operator in and let Ol 
be arbitrary elements of 0\ then there exist Olj G and complex coefficients Wj such 
that 

o^fnoi.)iz)=E^. fnoL|K)- (1-10) 

\fc=l / j \k=l J 

The Wj can be chosen such that they are valid for all representations in 7ls{C, C, x^ Xi -^)- 
The proof can be found in Appendix A. 

This result has an important corollary: Let Ok, k = 1,2, . . . , N he operators from the 
set O'^ U U IC. Then complex numbers qj and operators Oj j^ G exist, such that 

/ N \ 



I[Ok]\^) = T.<lAI[OlA\^). (1.11) 



^fc=l / j \k=l 



Again, the coefficients qk can be chosen so that they are equally valid for all representations 
in 7^5(C,C',x,X,A). 

To prove eq. 11.111 one simply has to note that, in the expression (jlk=iOkj \i) the 
operators belonging to /C are commutators (and therefore linear combinations of prod- 
ucts) of elements belonging to and O^, and that one can systematically eliminate all 
operators Ok E O^, by working from right to left using eq. 11.101 

If we take the inner product of eq. 11.111 with |/) we find 

{f\(I[0^\^)= E (1-12) 

\k=l / j(z{k:Nk=0} 



■^In fact, I was tempted to use the term 'highest weight' instead of 'foundational', but resisted because 
of two reasons: Firstly 'highest' suggests uniqueness, which does not have to be the case here and secondly, 
the structure here is a superalgebra, not a semi-simple Lie algebra, in which context the term 'highest 
weight' was coined. 
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since the operators Oj^, when present, yield zero if acted on by (/|. This is the main 
result from this section and might be restated as follows: 

The matrix elements with respect to the foundational states of products of generators 
in any representation in 7^5'(C, C, x, x, A) can be calculated solely from the commutation 
and anti- commutation relations of the superalgebra. These matrix elements are the same 
in all representations in TZs{C, C, x-, '^)- 

1.4 The Single- plus Bifermion Representation of the 
Superalgebra 

Up to this point we have not encountered any concrete representation of iS(C, C, 
but only postulated the existence of a representation. We now construct a concrete 
representation carried hy Tip- This representation will be used to define an operator T 
that maps state vectors in the fermion space Tip onto state vectors in the boson-fermion 
space T-Lbf- As was mentioned in the introduction, we can, with the aid of the operator 
T, construct a new representation of the superalgebra 5(C, C,x,x), this time in TIbf, 
where the generators of the representation are explicitly expressed as functions of the ideal 
boson and fermion operators. Both the fermion and boson-fermion representations lie in 
the same class 7ls{C, C, x, X, A) of representations. This fact will be used to demonstrate 
the invertibility of the mapping T. 

Let x^j,, i = 1, . . . M; /U, = 1, . . . , f2 be complex numbers such that xl^ = 
Define Xi'^ = (^m^^) ■ adopt the following notation: to indicate that the generators 
we construct belong to the representation of S{C,C,x,x) that is carried by Tip, the 
generators will be placed inside round brackets subscripted by F. Define 



It is assumed that 



IX^''X^. = 9SI (1.14) 
where g isa. positive real constant independent of i and j. With the aid of the orthogonality 



"^It will shortly become clear that this is not the only restriction we have to place on these numbers 
for the subsequent construction to work. 
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condition eq. II. 14^ (^Kj^ ^ when calculated explicitly, yields: 

{Ki)^ = g6l-x^%X<'^. (1.15) 

The reason for assuming the orthogonality condition eq. 11.141 is to ensure the linear in- 
dependence of the operators {Ai)p, i = 1, . . . ,M. Strictly speaking we could therefore 
have proceeded on the weaker assumption that the Xi'^ cire linearly independent rather 
than orthogonal. However, if they are linearly independent then there is an invertible 
transformation that will make them orthogonal. On the level of the algebra this merely 
represents a transformation of basis for the linearly independent operators Ai. All we 
therefore assume in enforcing the orthogonality condition, is that we have chosen a par- 
ticular kind of basis for the algebra, which represents no real loss of generality, since, for 
any algebra such a choice is possible. Note that since there are — 1) linearly inde- 

pendent operators a^a^, in Hp the linear independence of the {Ai)p operators places the 
upper bound M < ^fl{Q — 1) on M. We now have to investigate under which conditions 
eq. II.IHI is a representation of iS(C, C,x,x)- The only real work that is required per- 



tains to the closure condition 



Kj j ^ , {Ak)p = CJi^{Ai)p. We must check under which 



conditions coefficients Cjl exist that satisfy this equation. Since 

, {A,)p\ = xrxUl\a,)p{a,)p, (1.16) 
holds, enforcing the closure condition is equivalent to demanding that CH exist that solve 

This is a matrix problem of the form: find x such that Ax = y where a; is an M dimensional 
column vector, A is matrix with M columns and — 1) rows, and y is an |f2(f2 — 1) 

dimensional column vector. This system of linear equations may be over-specified and 
does not always have a solution. Thus there is a further restriction on the x's: they should 
be chosen such that the system of linear equations eq. 11.171 has a solution. For our work 
it turns out that it is not necessary to have a more explicit formulation of this restriction. 

When eq. 11.171 does have a solution, we can multiply the equation by ^xjj,, sum 
over fj, and ry to find^ 

cil = -gxrxui'xU- (1-18) 

From the fact that, in the present representation, operators with subscripts are trans- 
formed into their superscripted partners through hermitian conjugation, it follows that 
Cj'l = (Cil^ holds, and from eq. 11.181 that the C's are symmetric in both upper and 
lower indices: 

C^ik = = Cll = Cfcj. (1-19) 

^It must be stressed that eq. II .181 follows from eq. II. 171 but the reverse is not true. Therefore, simply 
choosing the C^l as in eq. II. 181 does not automatically take care of the closure condition. 
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It is straight forward to check ah other commutation and anti-commutation relations 
that the operators in eq. 11.131 should satisfy. We conclude that the operators in eq. 
II. 131 are indeed the generators of a representation of S{C, C, x, x)^ with Xi'^ = (x^i/) ^"^^ 
Cj'l = (Ci[j . We denote this representation J^. For ease of reference we define a set ^i;' = 

|(ajy)^; {a'^)p; (A)^?; {^^)f''' (-^Of)' "^^^ vacuum |0) in Hp has all the properties that are 
required of the foundational vector The functional (0| has all the properties required 
of foundational functional (/|. In particular = g6j. Therefore the representation JF is 
a member of TZsiC, C\ x, X^Q^D- 



1.5 What is Required of a Mapping Procedure? 

In the representation JF of the superalgebra iS(C, C*, x, x*), the operators {A^)p create 
fermion pairs. We shall call these operators collective pair creation operators, and their 
conjugates collective pair annihilation operators.^ We earmark these operators for asso- 
ciation with ideal boson operators. 

Let us consider what is required to map a fermion system in which the operators in 
Qf and {A^)p and {Aj)p in particular, are the building blocks of the physically important 
states and operators. Note that 

1. any state in Hp can be written (in more than one way) in the form 

l'^) = E'7.n(O,>)^|0), (1.20) 

i k=i 

where the rjj are complex numbers and are {Oj^kj ^ elements oi Qp- 

2. any linear operator {O^ ^ : 'Hp —* 'Hp can be written (in more than one way) in the 
form 

(o)^=E^.n(o,>)^, (1-21) 

j k=i 

where the rjj are complex numbers and are {Oj^f^j ^ elements oi Qp. 

3. any linear functional on Tip can be written (in more than one way) in the form 

(^l = (0|E^.n(O,>)^, (1.22) 

j k=l 

where the r]j are complex numbers and are {Oj^kj ^ elements oi Qp. 

^The term 'collective' is used because, when many x coefficients are non-zero, states such as {A^^ ^ |0) 
consist of a linear combination of many Slater determinants. However, we do not reserve this terminology 
for the case where many Slater determinants contribute. 
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Firstly, we must find a suitable isomorphism T : Tip Vbf C Ti.BF, that maps states 
in the fermion space onto states in the boson-fermion space. Then we have to find the 
image of the foundational vector |0) and an extension of the image of the foundational 
functional (0|. If we can then find extensions of the images of the operators in Qf and 
express them in terms of ideal boson and fermion operators, we see from the expressions 
for general vectors, operators and functionals in Tip, eq. 11.201 eq. 11.211 and eq. 11.221 that 
we are done. 

A remark on notation: An extension of the image under T of an operator (jDj ^ & Gf 
will be denoted (O^ It is related to (O^ ^ through 

{0), = T-^i0),X,J- (1-23) 

This is the only restriction on (cij and does not fix it uniquely, as the condition in 
eq. 11.231 onlv specifies the properties of (O^ in the domain Vbf- Outside this domain 
it is undetermined save that it should be linear.^ We therefore indicate the relationship 



between (O^ ^ and (pj as 



O), — (O),,. (1-24) 



to stress that, while there exists a linear mapping from the extensions (O^ to the 
operators 
direction. 



operators \0j^ stated in eq. I1.23| there need not be a linear mapping in the opposite 



1.6 Concerning the Usui Operator 

I now describe a method for finding an isomorphism T between Tip and the physical 
subspace Vbf C H-bf- We start by defining several operators in the space Hu = 'Hbf ® 
Tip: 

& = S'(g)I, 
B, = S,®I, 
a" = a" 

S; = a^(g)I. (1.25) 

Further linear operators on Tiu are defined via their action on a specific type of basis for 
Tiu: It is possible to find a basis Bu for Tiu such that each element of Bu has the form 
\ip,n) ® \4>) where \ip,n) is a vector in Hbf with a well-defined number of ideal fermions 

^The extensions will have to be linear if we want to express them in terms of ideal boson and fermion 
operators. 
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n, i.e. a'^ai,\ip,n) = n\ip,n), and \(f)) belongs to Hp- Define the images of (a^)^ and 
{au)p on such basis vectors as 

(o^^ n) ® \ct>) = (-1)" n) ® a-' |0) , (1.26) 

and 

ia:j}p \^|J,n) ® |0) = (-1)" |^,n) ® a, |0) . (1.27) 

Because these operators leave the number of ideal fermions unaltered, any product of 
an even number of them looses the (—1)" that appears in eq. 11.271 Thus, for example, 
{a^^)p{a'')p = \® {a'^)p{a'')p. Hence, if we define 

{A)p = '^xrMFi^)F- (1-28) 



then it is also true that 



i^F = ix;.I®(a'^)^(a^)p, 



(A,)^ = lxn®ia,)^{a^)^. (1.29) 



We conclude that the operators {Ai)p, (^*)p, {Ai)p , {A^)p , ia'u)p and {a'^)p obey the 
same commutation and anti-commutation relations as their counterparts without tildes. 
Furthermore, from eq. I1.29l it is clear that {Ai)p and {A^)p commute with Bj, , and 
a'^ . Since Bi and 5* do not alter the ideal fermion number, {ai,)p and {a^)p commute 
with these operators. However, {ay)p and {a^)p were defined in precisely such a way 
that the fact that 5^ and change the ideal fermion number by one results in the 
anti-commutation of {a)p^s and (a)^'s. 

The Usui operator U : Tijj Hu is defined as 

U = exp{B^{Ai)^ + a-{(^p) . (1.30) 

The Usui operator allows us to specify a linear operator T : Tip ^ TisF the following 
way: For each \(f)) in Tij? let |0) = T |0) be the unique vector in Ti^^ such that, for all 
IV') in Hbf 

(V.|0) = (^|®(O|/7|O)®|0), (1.31) 

holds. For the operator T to be an isomorphism between Tip and the physical subspace 
Vbf = TiTip) we need to show that T is 1 : 1. This matter is deferred for a while. 
Until then we must be careful not to construct arguments that assume the existence of 
T^^ : Vbf Ti-p- 

The advantage of defining T along this route lies in the exponential form of U, which 
suggests a systematic way of deriving (extensions of) the images of operators under this 
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F 



mapping. Suppose there exists an operator {O^ : TIbf ^bf such that [O^ 
(jDj where (O^ ^ G ^i? is a generator of the fermion representation of the superalgebra 
S{C, C, X, x)- Then it holds that 



from which follows 



F 



(^^1(6)^^10) = (v^|®(o|/7|o)®(d; 

= (^|®(O|/7{l®(6)^}|O)®|0) 

= (^|®(O|/7(6)^|O)®|0). (1.33) 

with and arbitrary vectors in Tip and TisF respectively and where = T \ . In 
the last line of eq. 11.331 I ® (O^ ^ could be replaced by (^O j ^ without worrying about the 
(— 1)" that appears in the definition eq. II. 271 because the vacuum |0) in Hbf contains no 
ideal fermions. The BCH formulas 

exp(P)Q = (^f;l(P,g),^exp(P), (1.34) 

/ CO /-^^k \ 

Qexp(P) = exp(P)(^5:^-^(P,Q)fcj. (1.35) 

where {P,Q)o = Q and {P,Q)k+i = [P , {P,Q)k] can now be used to find an operator 
O : Hu Ti-u such that OU = u{0^ ^ where O is expressed in terms of the {A)p, {a)p, 
B and a operators. From there on a few tricks are pulled out of the hat to find an operator 
[O^ : l^BF Ti-BF expressed in terms of the ideal boson and fermion operators, such 
that 

(V-l ® (0| {(O)^^ ^1}U |0) ® 10) = (^1 ® (0| OU |0) ® 10) . (1.36) 

These tricks involve, among other things, using the fact that some terms in O annihilate 
vectors like |0) ® |0) and functionals like ("01 ® (0| because they respectively contain the 



vacua in H-bf and Tip- The operator [Oj^^ that, thanks to this procedure, is expressed 
in terms of the B and a operators, satisfies the requirements contained in eq. 11.321 

A slight complication in this process is that it sometimes works better in reverse, so 
that for a given generator {O^ ^ in G we first have to guess the corresponding operator 
(O) and then work our way backwards from eq. 11.361 to eq. I1.33| thereby testing 
whether our guess is correct. 

In Appendix B the images of all the generators in ^i;' are calculated explicitly. The 
final result is: 

(o^i/)p iflu) BF ~ 
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J F 



(A, 



J) BP 



B 



31 



IBF 



a 



Kl 



BF 



BF 



B'lCi-\CilB'B^Bi + A\ 



2 

4k ^ 



alB'^Bi. 



;i.37) 



Here A\ Ai and /C] refer to the ideal fermion counterparts of {A^)p, (^i)^? and (Kj 



Ai 


1 




A' = 


2X.fiu^ ^ ) 




n = 




(1- 



1.38) 



Furthermore, ahhough the mapping is not unitary, the identities 

r|o) = |o), 

(0| ^ (0|. 



;i.39) 



do hold. In Appendix C, I show that these images comprise the generators of a represen- 
tation, denoted B, of S{C,C,x,x) in the whole Hbf, which justifies the notation (Obf' 
by showing that the desired commutation and anti-commutation relations are satisfied. 
This result is non-trivial, since the definition of the mapping eq. 11.321 of the generators 
says nothing about their properties outside the physical subspace Vbf, which we know to 
be a finite dimensional subspace of the infinite dimensional T-Cbf- 

Furthermore, even in the subspace Vbf, we have not yet proved that T is invert- 
ible. This has the following implication: Say that in the representation JF the following 
commutator holds: 

[(00..(0^)J=N.> 

where (Oij ^, (02)^ and (Osj ^ are generators of JF. Then because T is a linear operator 
and 



F \ F 

O] T, we know that 

BF ' 



BF 



BF 



T 



O 



3 1 T. 

•^1 BF 



;i.4i) 



However, we do not know if T is invertible in Vbf so that we can multiply with the inverse 
from the right, to establish that the commutator holds in Vbf- 

The fact that we do have a representation of iS(C, C, carried by Vbf therefore 
strongly suggests the invertibility of T. In fact, as a consequence of this, the invertibility 
of T is proved without much effort: Firstly, it is easy to establish that the foundational 
states |0) and (0| in TIbf have exactly the same properties as those possessed by the 
foundational states of the representation JF. This implies that both T and B belong to 
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TZsiC, C, X, X, gSi)- This in turn implies that for (Oi^ ^, i = 1,2, . . . , N generators of JF, 
and (Oij their extended counterparts in B the equality of matrix elements 

N N 

(0in(O.)^|0) = (0|n(O.)^^|0), (1.42) 

i=l i=l 

holds, as we have seen in Section fl.3l Let |0) be any vector in Tip- Then |0) can be 
expressed as 

l'/') = E^.n(O,>)^|0), (1.43) 

j k=l 



where r]j are complex constants and [Oj^k] are generators of JF. If 



j k=i 





+ 0, (1.44) 

where the order of the factors in the first product is now reversed. The representation 
T is such that if ^ is a generator, so is its conjugate ^ ■ Hence it follows that 
there are operators (Oj^k) '■ T^bf T^bf and (Off^ : Hbf T^bf such that 
ipj^kj^^T = Tipj^k)^, (O^^fc)^^r = T(dj^k)p holds, for which it is true that 




^ (oi(E^; n_ {o.k)^ I (E^.n(o,>)^ 

E^.nfe)«Jio)- (1-45) 



k=l 




In the last line of eq. 11.451 we can replace Vjllk=i {^j,k) ) |0) with |0) = T 



BF 



Furthermore (0| Vj Ilk=Nj (Ptk) bf^ some linear functional (in general not (0|, since 
T is not unitary) on TisF- Therefore eq. 11.451 savs that there is a linear functional on 
T-Cbf that does not map T |0) onto zero. Thus T |0) cannot be the zero vector in TisF- 
Thus if 10) 7^ then T \ (f)) ^ and T is invertible. 

Define the set Gbf = {Mbf'^ {o!')bf'^ {^^bf'^ {^')bf'^ (^Obf)' ^^^^^^^ ^ invert- 
ible, the operators in Qbfi together with the foundational states (0| and |0) comprise the 
essential elements of the mapping of systems in Tip to Vbf- 

However, consider the operators (A)i? = \xT i.(^u)F{'^tJi) f ^"^^ {'^i^)f ^^'^ ^^e exten- 
sions of their respective images: {Ai)^p = Bi and {a^)BF ~ '^t^- 

Clearly {Ai)^p = ^ Ai = ^Xi'^ {(^u) bf('^m-) bf^ even though, in the fermion represen- 
tation, it is true that {Ai)p = |xf''^(ai,)p(a^)^. Is this not in contradiction with the 
requirement that the images of operators should be related to the operators through a 
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similarity transformation? There is no contradiction: and {au)p need not be related 

to {Ai)^p and {ai,)^p through a similarity transformation. It is the restrictions of the 
latter to the physical subspace that are related to the former through a similarity transfor- 
mation. In other words, while Bi is not equal to Ai in TIbf, they are equal when restricted 
to the physical subspace, v4.i|y^^ = Bi\y^^. Consider for instance their respective actions 
on the state 

\i^) = {A^)^^\0) = gB^\0) + A'\0). (1.46) 

It is easily seen that B^ 1-0) = Ai = g |0) as it should be for a state in the physical 
subspace. 

Although there is no contradiction, the above discussion illustrates an undesirable 
feature of the mapping T. The association between the collective fermion pair annihilation 
operator {Ai)p on the fermion side and the ideal boson annihilation operator B^ is not 
as strong or exclusive as might be expected. The nature of the physical subspace is such 
that ideal boson annihilation is equivalent to ideal collective fermion pair annihilation. In 
the next section we attempt to address this problem. 

1.7 A Similarity Transformation 

Consider again the action of (A^)^p on the vacuum: 

{A^)^^\0) = {gB^+A^)\0). (1.47) 

So, the collective state {A^)p |0) is mapped by T onto a state, that although it contains a 
boson, retains a collective fermion component A^ |0). If we cannot get rid of the collective 
fermions through the mapping, then there is little profit in mapping, speaking from a 
computational perspective. Also, the intuitive notion that some fermion pairs are similar 
to bosons demand that we get rid of the collective fermion component in the mapped 
states. We have jumped the gun by settling on T as our isomorphism. 

Instead of the images in eq. let us rather postulate a mapping that results in 

extended operator images of the form 

(A^)^ ^ W =E^K;l-\CgB^B'Bi 

{k[)^ ^ {Kl)^^ = ^-CiiB^B,. (1.48) 

From the calculations in Appendix C, it is easily seen that altering the above operators 
in this way does not alter their commutators with each other. (What does change is 
the commutators of these operators with (a'^)^^ and {a,^)^p.) This inspires the notion 
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that there is some similarity transformation X so that the operators in eq. 11.481 are the 
transformed operators of the original mapped operators, 

X'^aA X = W = (A^) -A\ 

\ J BF \ J BF 



because similarity transformations preserve commutators (and anti-commutators). This 
can be viewed as a new isomorphism between 1-Lf and a subspace of TiBFi which is 
established by the operator X~^ oT, where 'o' indicates the conjunction of operators, i.e. 
X~^ o T = X~^ {T Under this new mapping, the images of the single fermion 

operators become X~^{a'^)^pX and X~^{ai,)^pX. These operators, together with those 
in eq. ll.49l form a new representation of S{C, C, x, x) in TisF- This representation will be 
denoted X. It is also a member of TZs{C, C, X, Xi 9^l)i this time with foundational states 
X~^ |0) and (0| X. The generators of this representation will be denoted [O] = 

X~^ {O^ ^^X . The new physical subspace is Vbf = ° T{Hf)- (We use the same 
symbol as before to denote the physical subspace.) 

I now prove the existence of an operator X such that 

BjX = XBj, 
(^R'+A^)X = XRK (1.50) 

If this X is invertible, it will perform the duties required of the similarity transform 
described above. We start by defining Cf = A^Ak-, which is an hermitian operator in 
TIbf- It is easily verified that Cf commutes with all the operators /C^. Since is 
hermitian, we can find a complete basis for 1-Lbf consisting of eigenvectors of Cf- A 
generic member of this basis will be denoted I'i/', A), with Cf I'V^, A) = A IV^, A). A linear 
operator X is defined by stating its action on this basis: 

oo 

X 1^, A) = E {(Cf - A)"^ A'B.y A) . (1.51) 

n=0 

Note that we do not have to worry about the convergence of the infinite sum. Since 
each term |(Cir — A)~^ ^'-Bjj creates 2n ideal fermions, and the fermionic sector oH-Lbf 
is saturated after the creation of VL ideal fermions, terms in eq. 11.511 with 2n > VL are 
identically zero. 

However, the definition is still problematic: since Cf has an eigenvalue A, Ci? — A has 
a zero eigenvector, which means that it is not invertible on the whole I-Lbf- However, 
Ci? — A is invertible on the orthogonal complement in TIbf of the A-eigenspace oi Cf- 
This orthogonal complement is denoted Aj^. It is therefore necessary that the operator 
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A^Bi maps every vector it acts on in eq. ll.51l onto A_l, in which case we take {Cp — A)~^ 
to be the inverse of C^? — A on the invariant subspace A_|_. I do not know of a result that 
shows under which general conditions A''Bi possesses this property, so that it has to be 
treated on a case by case basis. In Chapter Three we discuss an instance where it can be 
demonstrated that the similarity transformation is well-defined. 

Assuming then that eq. 11.511 does represent a bona fide definition of the operator 
X, 1 prove that it satisfies eq. 11.501 First notice that, if X) is an eigenvector of Ci;' 
associated with eigenvalue A, then so is Bj {ip, A). Therefore 



[X,B,]\i,,X) 



oo 

Y^{{Cf-X)-'A'B,^^^ ,5, 



.n=0 



l^,A) 



0. 



;i.52) 



Because l?/^. A) is an arbitrary member of a basis for T-Cbf^ this proves the first part of eq. 
11.501 The second part is proved by showing that 



Define, for convenience, operators 

XJz) 



X,W =A^X. 



{{Cr-zr'A'B,^''^ 



;i.53) 



;i.54) 



where 2; is a complex number and n = 0, 1, 2, . . . Now define linear operators Xn by means 
of their action on the basis vectors X): 



Xn\^,X)=Xn{X)\ij,X). 



;i.55) 



From this definition it follows that X = J2'^=oXn. By induction 1 now show that 



for 77, = 0, 1, 2, . . ., so that 



;i.56) 



X^W = Y.[^n+l.R^' 
n=0 

00 

= E^'^n 

n=0 

= A^X. 



;i.57) 



We shall often refer back to the following two intermediate results: Firstly 



Cf.R' =0, 



;i.58) 
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so that if I'i/'jA) is an eigenvector of Cp, with eigenvalue A, then so too is W \ip,X)- 
Secondly 



Now we can start the proof proper, by setting n = 0, then 

'X,,R^] \ij,X)=\x^{X),R^] |^,A) 



;i.59) 



;i.6o) 



where the expression on the right follows because both X) and R^ A) are eigenvectors 
of Cp with eigenvalue A. Continue by substituting eq. 11.591 into eq. I1.60l to find 



MX),R^]\tp,x) = {Cp-xy' 



Cf,A' 



\i',x) 



{Cp-xy' iCF-x)A'\ij,x) 

A^ 1^, A) . 



:i.61) 



We conclude that [Xi , R^] = A^ holds, which is what we wanted to prove, since Xq = 1. 
Take the next step, by assuming eq. 11.561 is true for n = k — 1, i.e. 



Xh , -R"' 



:i.62) 



Then 



Xk+,,W\ \ij,X) 
Xk+i{X),R'] |^,A) 
X,{X)Xk{X),R'] \ij,X) 



{[x,{X),R^]Xk{X) + Xi{X) [Xk{X),R^]) \ij,X). 



Now we take the induction step by substituting from eq. 11.631 



Xk+i , R'' 



A) = ( MX) , R' Xi(A) + Xi{X)An Xk-iiX) A) . 



:i.63) 



:i.64) 



The indicated factor can be simplified. Firstly, substitute into it the explicit expression 
for the commutator, that was derived in eq. I1.59| and also the explicit expression for 
Xi{X). This yields 

Xi(A) , R^] li(A) + Xi{X)A^ = {Cf - Xy^ {[Cf,A^] {Cf - A)"' + A') A'S^. (1.65) 

Here we need a little sleight of hand: we write the commutator appearing on the right as 
follows: 



Cf,A' 



{Cf- X)A' - A'{Cf- X), 



;i.66) 
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from which it follows that 

{Cf - A)-^ [Cf , A^] {Cf - A)-i = A^{Cf - X)'^ - {Cf - A)^M^'. (1.67) 

If we substitute this back into eq. 11.651 we find 

Xi(A) Xi{X) + Xi{X)A' = A'XiiX). (1.68) 

This, substituted back into eq. 11.641 shows that if eq. 11.561 is valid for n = /c — 1, it is also 
valid for n = k, which completes the proof by induction. 

Next we must show that X is invertible. In order to do this, define recursively 

Yo = I, 

n+l 

= -Y.X,Y^+,_i. (1.69) 

1=1 

The operator X„ creates 2n ideal fermions when 2n < Q. If 2n > f2 then X„ is the 
zero operator. From this follows by induction that Yn creates 2n ideal fermions (and is 
therefore zero when 2n > Q): Say it is true for all Yk, k < n, then from the definition eq. 
II. 691 it follows that Yn+i creates 2{n + 1 — 1 + 1) = 2{n + 1) ideal fermions, which completes 
the proof by induction. 
This means that 

oo 

Y=Y.Yn (1.70) 

n=0 

contains only a finite number of non-zero terms, because all terms with 2n > Q are 
identically zero. Y is therefore well-defined (provided of course that the operators X„ are 
well-defined, a matter that was touched upon previously). A similar argument, based on 
the fact that A„ annihilates n ideal bosons, shows that Yn also annihilates n ideal bosons. 
Now we calculate 

Coo \ / °^ \ oo n 

^ ^ Xm I I ^ ^ Yfi j ^ ^ ^ ^ XYnYfi—m 

m=0 ) \n=0 / n=Om=0 

oo / fc+1 \ 

= 1 + ^ n+i+^A,n+i_^ , (1.71) 

A:=0 V m=l / 
" V ' 

=0 

SO that Y = X~^, i.e. X is invertible. In the last line of eq. 11.711 we took the m,n = 
term XqYq = I out of the summation, set A; + 1 = n and used the fact that Aq = / so 
that AoYfc+i = Yfc+i. 

Utilizing the recursive definition of Yn, one finds the various terms of X~^ to be 



Yo = /, 
Y, = -Ai, 
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Y2 — —X2 + XiXi, 

I3 = —^3 + X1X2 + X2X1 — XiXiXi, 

n 

Yn = E(-r E X,,X,,...X,,^, (1.72) 

m=l pePmn 

where Pmn = {p = (pi,P2, • • • ,Pm) : 1 < < n + 1 - m; pi + p2 + • • • + Pm = ?^}- 

The foundational states in this mapping are |0) = |0) and (0|X = (0|. To 
summarize then, if X as defined in eq. 11.511 is well-defined, we have an isomorphism 
X~^ o T between Tip and a subspace Vbf C TY^f of the ideal boson-fermion space. The 
essential elements of this mapping read 



|o) - 

(0| - 


- |0), 

- (o|, 








: Bj, 




V Jx,BF 


-- B'lq - \CilB'B''Bi 


(^). - 


V ' J X,BF 


= iq-cilB'Bi, 




— {ttiy)x,BF ~ 


X~^auX, 




~ ('^ )x,BF ~ 


X~' {a" + x^B'a^) X. 



(1.73) 

It is hard to find more explicit formulas for the extended images {0'u)xbf ^^"^ i^'^)xBF 
in this general setting, firstly, because X is defined partially by stating its action on a 
complete basis, rather than by completely expressing it in terms of other known operators, 
and secondly, because of the cumbersome formula for X~^. In some particular cases 
however, more concise expressions may be derived. 

When calculating the image X~^ N) where {ip, N) is a vector with a well-defined 
number of ideal bosons A^, a simplification does occur, thanks to the structure of X^^. 
Since the operator annihilates k ideal boson, the only terms in X~^ that contribute 
are those with such that k < N. In particular, since T{a'^)p |0) = {ay)gp |0) = |0) 
contains no ideal bosons, we have 

K)i.|0)''-^^«nO). (1.74) 

The image of {a^) p{a'^) p |0) contains one ideal boson. It is therefore not to hard too cal- 
culate what X~^ does to this image. Only the terms Yq = \ and Yi = —Xi contribute. Let 
us therefore calculate it in detail. Firstly, note that T{a^)p{a'^)p |0) = {a^a'^ + Xj'^ B^) |0). 
Therefore o T(a^)^(a'^)^ |0) = {a^'a" + Xj'' B^) \0)-xfXiB^ |0). Since CpB^ |0) = 
we have from the definition of Xi that XiB^ |0) = Cp^A'BiB^ |0) = Cp^A^ |0). We 
know that Cp = A^Ak so that CpA^ |0) = A^ \Ak,A^] |0) = gA^ |0). Hence Cp^A^ |0) = 
^A^ |0). The final result then is 

{a^pia^F |0) ""-^^ («^«^ + XTB' ' "Xf^^) |0) • (1-75) 
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Note that when we multiply by ^xfiu sum over /i and u we find ^xjiui'^^) pi'^'^) f I^) ^ 
gB^ |0) as we should. 

It is also not hard to compute (0| {ay)^pX\ for an arbitrary vector 1-0), consider 
(0| {au)QpX 1-0) = I]^o (0| '^uXn Every term ayXn, with n > 1 creates 2n — 1 ideal 
fermions so that a^Xn I'i/'), n > 1 has zero overlap with the vacuum. The only surviving 
term has n = 0. Hence (0| (aj,)^^X 1-0) = (0| It follows that an extension for the 

image of the functional (0| {ay)p is 

(0|(a,)^-^^^(0|a,. (1.76) 



1.8 Alternatives to the General Similarity Transfor- 
mation 

The similarity transformation of the previous section has the draw-back of being compli- 
cated. As a result it is not clear how to express the mapped single fermion operators as 
functions of the ideal boson and fermion operators. Also, and perhaps more alarmingly, 
we could not prove the general validity of the definition of the similarity transformation, 
as has already been discussed. Because of this state of affairs, I now show an alternative 
route that leads to the mapping of eq. I1.7HI in certain restricted circumstances. 

These circumstances occur when, instead of mapping the whole space Tip we only 
want to map a subspace Tip^coh which is spanned by vectors of the form 

fn M^)K)i.|0), (1-77) 

\m=l / 

and 

fn M^)|0)- (1-78) 

\m=l / 

We refer to T-Cf,co1 as the collective fermion subspace. States of the form eq. 11.771 are said 
to belong to the odd subspace of Hf,coI because they contain an odd number of fermions 
and similarly states such as eq. 11.781 are said to belong to the even subspace. Often, a 
system can effectively be described by a Hamiltonian that leaves the collective fermion 
subspace invariant. Furthermore it is often then the case that the states associated with 
the lowest part of the Hamiltonian's spectrum lie in the collective fermion subspace. The 
argument in this section is developed with such a situation in mind. 

Under the mapping that we denoted by T and that results from the Usui operator, 
the states in eq. 11.771 and eq. 11.781 are mapped onto 




(1.79) 
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and 

I f[ W"^ + A^A \0) , (1.80) 

\m=l J 

where we define W = B^}C{ — ^CjlB^B''Bi as previously. Note that does not change the 
number of ideal fermions. What we want to achieve, and for which we had to introduce 
the similarity transformation X in the previous section, is removing the collective ideal 
fermion pair creation operators A^ from the above expressions. However, in the present 
restricted case, we see that this can be achieved very simply. We just take the projections 
onto the subspace with less than two ideal fermions, of the states in eq. II. 791 and eq. 11.801 
We denote this projection operator P. 

Perhaps the reader is troubled by the use of a projection operator. After all, the 
mapping from TYf.coI into T-Cbf must be invertible and we know that projection operators 
aren't invertible. However, although P is not invertible on the whole Ti-BP, it is invertible 
on the range TiTip^coij- I demonstrate here the equivalent statement that PoT is invertible 
on the whole collective fermion subspace Tip^coi- The proof is similar in spirit to the proof 
of invertibility of T. 1 only show the proof for odd states. For even states the proof follows 
the same route. ^ 

We construct an arbitrary non-zero state |0) from the odd subspace of T^f.coI by choos- 
ing appropriate indices Ij^k ^ {li 2, . . . , M}, uj G {1,2, ... ,Q} and complex coefficients 
r]j and setting |(/)) = EjVj (llfii i^^''") p) (^''Of |0)- Then ^ (0| (p) holds. The inner 
product (01 0) is explicitly given by 

(01 4>) = E^.v (oi in KO.l I n (^''Aj JO) • (1-81) 



The properties of the mapping T allow us to replace the right hand side of the previous 
expression with 




E V* (oi n Bi,,, E n R'^'^' + "'^^ io) 



:i.82) 



where we have also separated the T-Cbf inner product into a bra and a ket with curly 
brackets. Since the bra state contains only one ideal fermion, the collective fermion pairs 
A^ may be omitted from the ket state. The ket, with these pairs omitted is equal to 



^Why do I not simply say "When is restricted to T{'Hf.co\) it is identical to P, the projection 
operator onto the space with fewer than two ideal fermions, restricted to the same domain." Answer: 
Because X and have definitions that may be problematic, the cleanest view to take, is that when 
X can be defined as in the previous section, X~^ is an extension of the restricted P operator to the 
whole Hbf- The definition and properties of P restricted to T{Hf,coi) should not rely on those of the 
sometimes trustworthy X. 
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P o T\(f)), while the bra is just some linear functional on T-Cbf, which we'll call / for 
brevity's sake. We have then deduced that 7^ (0| 0) = / (P o T This means that if 
10) is not the zero vector in Tip^cox then P o T |0) is not the zero vector in T-Cbf and hence 
P o T is invertible on the odd subspace of Hf,co1- The same result can be obtained for the 
even subspace, so that we can conclude that P o T is invertible on the whole TYf.coI- 

Now that the isomorphism between states in T-Cf,coI and states in a subspace of TisF 
has been established, we can find the extended images of operators under the mapping. 
Obviously only those operators that leave Hf,coI invariant may be mapped. The operators 
that leave Hf,coI invariant are {A^)p, {Ak)p and (^Kl^ ^, i.e. the even sector of the super- 
algebra so that only these operators can be said to have images under PoT. A calculation 
that is too trivial to reproduce here give these images as 

{k{)^ ^ iq-CgB'Bi, (1.83) 

the same as that in the previous section. 

In passing it may be noted that, had we worked with only the even subspace of 7iF,coi, 
we could have used the projection operator onto the subspace of T-Cbf that contains no 
ideal fermions, instead of P. The mapping would then have simplified to 

ICllB^B'Bi, 

CgB'Bi. (1.84) 

This same mapping could also have been derived by a much simpler route, that starts 
with an Usui operator U = exp (^B^{Ai)p?j, without ideal fermions. 

1.9 Finding the Physical Subspace 

PoT 

In the preceding sections we discussed two situations, firstly a mapping TYf.coI — ^ T^bf 
from the collective fermion subspace into the ideal boson-fermion space and secondly a 
mapping Tip ^ — T^bf with a larger domain. It was noted that oT, when it exists, 
is an extension of P o T from the domain T{T-Cf,coi) to the a larger domain TiVp), where 
Vp is a subspace oi Tip and might in fact be the whole Tip. In what follows we will 
employ a new notation that covers both these situations. We denote the subspace of Hp 
that is being mapped (and which may be the whole Tip) by Vp. The mapping itself (i.e. 
the operator X^^ oT or PoT) is denoted T. The range of T, which is always a subspace 
of Hbf, is called the physical subspace and is denoted Vbp = TiVp). 
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As was previously explained, for operators Op '■ Vp Vp that are expressed in terms 
of products of the generators of the fermion realization of a certain superalgebra, we were 
able to find corresponding operators Obf '■ T^bf 'Hbf, expressed in terms of ideal 
boson and fermion operators, and are such that 



O 



BF 



= TOFr-\ (1.85) 

Vbf 



where |vbf indicates a restriction of the operator's domain to the physical subspace. An 
operator that satisfies eq. 11.851 should properly be called a physical boson- fermion opera- 
tor. For brevity's sake we will use the less descriptive term of physical operator. Physical 
operators should not be confused with the extended images of observables, which, it is 
true, form a proper subset of the set of all physical operators. (In the original fermion sys- 
tem, an observable is any hermitian operator acting in the fermion space.) Note that, as 
an immediate consequence of eq. 11.851 any physical operator leaves invariant the physical 
subspace Vbf- When we analyse a fermion system by investigating its mapped coun- 
terpart, then the fermion operators are replaced by physical operators. Although these 
operators are defined on the whole I-Lbf-, only their action on the physical subspace is of 
relevance. An investigation of the mapped system must therefore always be accompanied 
by an identification of the physical subspace Vbf- 

In this section we discuss techniques for identifying the physical subspace. The most 
obvious of these is the following. If we have a basis for Vf, say then by applying 

the mapping T to each we find a basis {|0i) = T|</)j)}^^ for the physical subspace 
Vbf- However, it turns out that this procedure leaves us with a basis for the physical 
subspace that is complicated to work with. As 1 hope to make clear by considering 
applications of the Dyson mapping in the next two chapters, it is impractical to work 
with the above basis. It entails dealing with structures that so closely resemble the 
immediately evident structures of the original system, that it is hardly worth the trouble 
to do a mapping. If one works in the above basis, it is fair to say that the mapping changes 
the analysis of the system under consideration only cosmetically. Essentially the same 
objects are manipulated, they are only represented by different symbols. The rest of this 
section will therefore be concerned with more sophisticated procedures of characterising 
the physical subspace. 

A basis such as the one above is called a physical boson-fermion basis. It is to be 
contrasted with the ideal boson-fermion basis that consists of states the form 

Ha^'lO). (1.86) 

I 

When analising operators that are defined in T-Lbf-, it is natural to use this ideal boson- 
fermion basis, and only afterwards identify the physical subspace. An objection might be 
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that whereas Vp and therefore also Vbf are finite dimensional, T-Cbf infinite dimensional. 
Since working in infinite dimensional spaces is hard, it seems a good idea to identify a 
finite dimensional subspace that still contains all the physical information before analysing 
the system. This is true. However, it is possible to find such a finite dimensional subspace 
that contains the physical subspace, but has the advantage that ideal basis states may still 
be used rather than physical basis states. This is possible because the fermion systems 
we consider may only contain a finite number of particles. A property of the mapping we 
consider is that if a fermion state \ifj) contains fermions then its mapping S \ ifj) is an 
eigenstate of 2Ni) + Nf with eigenvalue N, where Ni, counts the number of ideal bosons 
and Nf counts the number of ideal fermions. Therefore, if the fermion subspace Vp only 
contains states with Ni, i = 1,2, ... ,n fermions, instead of working with the whole ideal 
boson-fermion basis, we can select only those basis elements for which the twice number 
of ideal bosons plus the number of ideal fermions equals one of the Ni. 

When Vp is the collective fermion subspace TYf.coI we can use an even stronger cri- 
terium to select the physically relevant ideal boson-fermion basis states. Because mapping 
involves projection onto the component oH-Lbf with fewer than two ideal fermions, the 
ideal subspace is further reduced to contain states with at most a single ideal fermion. In 
this instance it is often found that the ideal boson-fermion basis states thus selected in 
fact spans just the physical subspace Vbf and nothing more. The space spanned by ideal 
boson-fermion states that are selected with the aid of the above particle number consid- 
erations is called the ideal subspace. It contains the physical subspace as a subspace. 

Typically, analysing a physical system means diagonalizing one or more physical op- 
erators. The procedure we expound here is the following: 

1. Diagonalize the physical operators using the ideal boson-fermion basis states that 
fulfill the correct particle number requirements, i.e. that lie in the ideal subspace. 

2. From the eigenstates thus obtained, select those that span the physical subspace. 

How the second step in this procedure is implemented, depends on the application we 
have in mind and on the nature of the fermion space Vf- I discuss two procedures. 

The first procedure utilizes the fact that physical operators leave the physical subspace 
invariant. The argument 1 present here implements an idea that can be traced back to 
a paper ^2] by Geyer and co-workers. What 1 present here does not contain any new 
fundamental insight that goes beyond the content of ^Hl- However, 1 present in greater 
detail the actual algorithm that one would have to implement to utilize the idea of Geyer 
and co-workers. The procedure is based on an algorithm for solving the following problem: 

Problem: Given are two linear operators A : V V and B : V V where V is 
a finite dimensional vector space with dimension M. A is assumed fully diagonalizable 
with a non-degenerate spectrum. Find a non-trivial subspaces (if one exists) of V that A 
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and B both leave invariant. (Non-trivial here means a subspace of V that is neither V 
itself, nor does it only contain the zero vector.) 

Solution: Let \(pi), i = 1, . . . , M be the eigenvectors of A, with corresponding eigen- 
values Aj. We have assumed no degeneracy so if 2 7^ j then Aj 7^ Xj. Let (■i/'jl, i = 1, . . . , M 
be the left eigenvectors of A such that {ipi\A = Aj {ipi\. It is easy to prove that {i/ji \(f)j) 
equals zero if and only if z 7^ j. By choosing the left eigenvectors to be properly normalized 
we may therefore, without loss of generality take 

(^, 10,) = 5ij. (1.87) 

This gives rise to the unity resolution in V 

M 

1=1 

Now suppose that W is an A^-dimensional subspace of V that A and B both leave invari- 
ant. Because A leaves W invariant, there exists a basis for W that consists of of the 
eigenvectors of A. Let Cw be this basis. Then we have the following theorem, the proof 
of which is trivial: 

Theorem 1 If lij) is a vector in W, which has the unique expansion 

M 

|^) = Ecfc|0A:), (1.89) 
k=l 

in terms of the eigenstates of A, then for each Ck 7^ 0, \(pk) is a vector in Cw- 

Define a matrix Bji = [i/jil B \ (pj). Thanks to the unity resolution eq. I1.88l it holds that 

M 

i? 10,) = E 5,, !</>,). (1.90) 

i=l 

A corollary of Theorem Q can now be stated. 

Corollary 1 // \ belongs to a subspace that A and B both leave invariant and Bji 7^ 
then also belongs to that subspace. 

Again the proof is trivial. 

Using the tools we have gathered we can formulate an algorithm that enables us to 
find some non-trivial subspaces that A and B both leave invariant, if any exists. If (and 
only if) none exist, the algorithm will find only the trivial subspace V: 

Assume there exists a non-trivial subspace W that both A and B leave invariant. W 
has to contain an eigenvector of A, because A leaves W invariant. Therefore, guess that 
W contains If this guess is wrong, the algorithm 1 am about to describe will return 
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W = V . If this happens, we simply try another eigenvector of A as belonging to W . If 
we have exhausted all options and the algorithm still returns W = V, then we'll know 
that no non-trivial subspace of V exists, that A and B both leave invariant. Define sets 
of indices recursively as follows: Ij'^^ = {j}. For k = 1,2, . . ., let I^^^^^ be the set of all 
/ G 1, . . . , M such that there exists an i G /] for which Bu ^ 0. To see what 1 am up to, 

it holds that Bji 



note that, since | (;/),•) belongs to W ^ if z G it holds that Bji ^ and therefore it follows 



from Corollary ^ that is also an element of W . Then, if z G if^ and / G it holds 
that Bii 7^ 0. Therefore, because |0j) is in W ^ from Corollary ^ follows that so too is 
Thus far we conclude that if z G if^ U if^ U if^ it holds that |0j) G W . For any k define 
Uf = Uto-^f • Inductively we arrive at the result: For any /c, if z G Uj'^^ then belongs 
to W. Furthermore Uj'^^ is a subspace of Uj''^^\ or in other words, if z belongs to Uj^^ 
then it also belongs to Uj''^^^ , but the number of elements in any Uj^^ is at most M. This 
means that at some k, the set f/j*^^ must stop growing. More formally put, there always 

exists a kj such that t/j^^ = uf^'^ for all k > kj. If we define Cj = : z G 

then there exists a basis for the invariant subspace W that contains all the elements of 
Cj. The result we have arrived at may also be stated: If W is a subspace that A and B 
both leave invariant, and W contains the eigenvector \(pj) of A, then W contains all those 
eigenvectors of A that form the set Cj. 

Now we show that A and B both leave span{Cj) is itself invariant. Since span{Cj) 
has a basis Cj that is comprised of eigenvectors of A, we know that A leaves span{Cj) 
invariant. To prove that B leaves span{Cj) invariant, assume that it is not so. Then there 
exists an index i G L^j'^^^ such that B |0j) does not lie in span{Cj). This is only possible if 

(fc) (k) 

there is an index I such that Bu ^ but I ^ Uj . Since z is a member of Uj , it must 
be a member of some jj'"'' with m < kj. Suppose m < kj, then because of the recursive 
definition of the Ij^^ and because Bu ^ 0, it holds that / G /j'"^^^ Since m + 1 < kj, 
it holds that /j"^^^^ is a subset of Uj^^^ and therefore I is an element of Uj^''\ which is a 
contradiction. Suppose on the other hand m = kj. But then, since Bu 7^ 0, the index / is 
an element of which is a subset of uf'^^^\ But kj was defined in such a way that 

_ jjif^j) _ 'p]-^^g again we are led to the contradiction that / G Uj''^\ 

The above results imply that both A and B leave span{Cj) itself invariant and that 
span{Cj) is minimal in the sense that every subspace of V that contains \(j)j) as an 
element, and that A and B leave invariant, contains span{Cj) as a subspace. If we run 
this algorithm, starting with j = 1 and continuing until span{Cj) ^ V or j = M, we will 
find a nontrivial subspace of V, that A and B leave invariant, if such a subspace exists. 
If none exists, and only then, will we reach j = M and span{Ci) = V for z = 1, 2, . . . , M. 

I shall now explain how we can use the above algorithm to find the physical subspace. 
The key insight is that all physical operators leave the physical subspace invariant. To 
implement the algorithm, we set V equal to the ideal subspace. For the operator A 
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we choose a physical operator that has a nondegenerate spectrum in V. If degeneracy 
is a problem, we may simply use a complete set of commuting physical operators that 
we simultaneously diagonalize. Instead of using only one operator B we choose a set of 
physical operators B^, a = 1,2, . . . , N that is large enough that they leave no subspace 
other than the physical subspace invariant. Using a set of operators {Ba}'^^^ simply 
requires that in the algorithm we change the recursive definition of the index sets jj'^'* to: 
Ij''^^^ is the set of all indices I such that, for some i G Ij''^ there exists an a such that 
{Ba)ii 7^ 0. Since we assume that the only subspace that both A and all the B^ leave 
invariant, is the physical subspace, this is the only possible subspace that our algorithm 
can yield. 

1 fear that the above description of the algorithm looks very formal and that the simple 
idea behind it is obscured. Therefore 1 roughly summarize. To find the physical subspace 
we first diagonalize a physical operator A that has a non-degenerate spectrum. A sub-set 
of its right eigenvectors spans the physical subspace. Therefore we make a guess and 
take one particular eigenvector of A and assume that it is in the physical subspace. The 
algorithm will let us know if our guess is wrong. To find the physical subspace we now 
take a physical operator B that we assume is such that the only space that both A and 
B leave invariant, is the physical subspace. If a single operator B does not do the job, we 
can, as was explained, extend the procedure to use several operators B^, such that, taken 
together, the only subspace that both A and all the B simultaneously leave invariant, 
is the physical subspace. For simplicity, in the rest of this summary we assume a single 
operator B is sufficient. We now apply the operator B to the eigenvector of A that we 
assumed to lie in the physical subspace, and expand the result in terms of the eigenvectors 
of A. Every eigenvector of A for which there is a non-zero expansion coefficient in this 
expansion must also lie in the physical subspace. We now repeat the process of applying 
B and expanding in terms of the eigenvectors of A to each of these newly found elements 
of the basis for the physical subspace. If non-zero expansion coefficients for eigenvectors 
of A, different from those already taken to be inside the physical subspace, are generated, 
these eigenvectors are also taken to be inside the physical subspace, and the process is 
repeated for them too. We continue in this way, until non-zero expansion coefficients are 
generated only for states that are already taken to be inside the physical subspace. If, at 
this point, what we take to be the physical subspace does not include all the eigenstates 
of A, our original guess has been vindicated, and we have found a basis for the physical 
subspace. If, on the other hand, we are left with a result that says all of the eigenstates 
of A lie in the physical subspace, this either means that the state we initially guessed to 
be in the physical subspace is actually not in the physical subspace, or there simply are 
no ghost states. To eliminate the first possibility we repeat the process starting with a 
different initial guess. If, after all possible choices for the initial guess, we still end up 
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with the algorithm saying that all of the eigenvalues of A belong to the physical subspace, 
then we conclude that this is indeed the case. 

The astute reader may have noticed that several things can go wrong in this procedure. 
The root of the problem is that, assuming the physical operators are obtained through 
mapping, we only have control over their behaviour inside the physical subspace. What 
they do in the rest of the ideal subspace is determined by the specifics of the mapping. 
So, while it is true that we may find physical operators that satisfy all our requirements 
inside the physical subspace, i.e. no degeneracy for A and irreducibility for the Ba together 
with A, whether they still do this in the whole ideal subspace is not within our control. 
There might be ways around these problems, but 1 do not consider them because, in the 
application 1 will discuss, we will not run into any trouble. 

We now come to a second method for identifying the physical subspace, that goes by 
the rather uninformative name of '7^-projection'.^ The way that this method is presented 
below looks superficially different from how it was originally presented by Dobaczewski 
in jB], but the method is in fact identical to the 7^-projection introduced there. We shall 
see that this procedure is not as widely applicable as the first procedure. It can only be 
implemented when the fermion space that gets mapped is the collective subspace Hf,coI- 
To discuss the procedure effectively we first have to define a few operators and elaborate 
on certain of their properties. 

Let Cbf the canonical mapping on T-Cbf- This means that Cbf maps a general 
vector \if)) in Hbf onto the functional (-01, where the meaning of [ipl is fixed through the 
boson-fermion inner product defined on Hbf, as we have explained in a preceding section. 
The operator Cbf is anti-linear. Let r be the linear operator that acts on the space of 
linear functionals defined on Ti.BF, that maps a general linear functional onto {ifj\T, 
where T is the boson-fermion mapping. The operator r maps linear functionals defined on 
TisF onto linear functionals defined on the fermions space Vp which is the domain of T. 
Lastly, let Cp^ be the inverse of the canonical map on T-Cf- Cp^ is an anti-linear operator 
that maps a general functional (0| onto the vector where this time, the relationship 
between (0| and |0) is determined by the fermion inner product defined on T-Cp- 

We now form a linear operator TZ through the conjunction of the above operators: We 
define 

n = Cp^ otoCbf- (1.91) 

Since, of the operators used in the definition of TZ, two are anti-linear, and the other 
linear, TZ itself is linear. It is an operator that maps the space Hbf onto Vp, the domain 
of the mapping T, which is a subspace of fermion Fock space. This statement can be 
verified simply by keeping track of the domains and ranges of the operators used in the 

^Tlie symbol '7?,' has no special significance and, strictly speaking, the procedure does not involve 
projection. 
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construction of TZ. 

It will now be demonstrated that the operator TZ has the property that 7^ I?/;) = if 
and only if \ip) belongs to the orthogonal complement Vbf± of the physical subspace Vbf- 
First we prove that TZ \ifj) = implies that lies in the orthogonal complement of the 
physical subspace Vbf,±- Let |0) be any vector in Vp. Now we compute the inner product 
in Hbf between TZ I?/;) and the arbitrary vector 



7^ 



Cp'rC 



BF 



[1.92) 



In the above equation, by the respective definitions of r and Cbf, the functional tCbf I"^) 
is equal to T . Furthermore (0| Cp^ maps functionals defined on the fermion space Vf 
onto the complex numbers in the following way: Let (x| be any functional defined on Vf 
then 



(0 Ix) , 



Hence we arrive at the equation 



:i.93) 



:i.94) 



in which the right-hand side is an inner product of with a vector T |0) which lies in 
the range of T, i.e. the physical subspace. Hence, if |^) lies in the orthogonal compliment 
Vbf± of the physical subspace {(j)\ 7Z \ip) equals zero, for arbitrary so that 7Z must 
be zero. Now the reverse is proved, namely if 1-0) is not an element of Vbf± then 7Z\ip) is 
not zero: Take a |^) that is not an element of Vbf±- Then it can be written uniquely as 
IV') = "^ll) + \i'±) where -^n) is a nonzero element of the physical subspace Vbf and {i/jx) 
is an element of Vbf±- Since -i/^nj is a non-zero element of the physical subspace Vbf, 
there exists a non-zero vector 



ip\\) in the fermion space Vf, such that T ip^n = t/j^n. If 



we compute the inner product between 



and TZ 



we find: 



i^\l\Cp'TCBF{\i^\l) + \^P±)} 



;i.95) 



On the right-hand side of this expression, the component \iIj±) can be omitted as it is 
killed off by TZ. This leaves us with 



7^ 



T 

(^11 1 V'l 



Cp^tCbf V'll 



;i.96) 



which is nonzero. This completes the proof. 
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We have seen that at our disposal is an operator 71, that has the orthogonal comple- 
ment of the physical subspace as its null space. Before we can use this operator though, 
we must be able to evaluate its action on a generic basis element of the ideal subspace. 
This evaluation, as it turns out, is only feasible when the fermion space that is mapped, 
is the collective subspace 'Hr,co\- In this case, the ideal subspace only contains states with 
fewer than two ideal fermions. In what follows, 1 therefore evaluate 'R-Ylj (B^)'^^ |0) 
and briefly indicate what goes wrong when more than one ideal fermion is present. We 
start from the definition of 7^ in eq. 11.911 

i j 

= C^F'|(0|«.n(^.f^^|- (1-97) 

When only one ideal fermion annihilation operator is present in the last line of the previous 
expression, and only then, can the functional in curly brackets easily be evaluated. The 
key is to realize that when we restrict the functional {0\ ayHj (Bj)^^ to the physical 
subspace, as we may when we stick the operator T in front of it, we have 



a 



j 



a. 



ly J p 



EiA,)/^r~\ (1.98) 



If more than one ideal fermion were present, the similarity transformation that is part 
of the mapping would have wreaked havoc with the ideal fermions, churning out states 
that contain high numbers of non-collective fermion strings (^/^^^(ct^a)^ • • • ('^Mn)^- 
us therefore remain safely inside the collective subspace and substitute the most recent 
result eq. II. 981 into the expression we have for TZ acting on an ideal subspace basis element, 
namely equation eq. 11.971 




Thus, when the fermion space that we map is the collective subspace TYjr ^oi, we have a very 
simple prescription for evaluating the expression TZYlj (B^)^^ |0): Simply replace each 
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with the corresponding collective fermion operator {A^)p and replace the unpaired 
ideal fermion with a real fermion. 

To find the physical subspace using the TZ operator, we go about as follows. Again we 
assume that at our disposal is a physical operator A that is fully diagonalizable in the ideal 
subspace. We assume that A has a non-degenerate spectrum. (If degeneracy occurs, the 
procedure I am about to describe may trivially be generalized to accommodate, instead 
of a single operator A, a complete set of commuting physical operators, that remove any 
degeneracy.) In analogy with what we did previously, we diagonalize A to find say M 
eigenvectors |0j) associated with eigenvalues Aj such that ii i ^ j then ^ Xj. Similarly, 
a diagonalization of A"^ yields M eigenvectors {i/ji) associated with eigenvalues A*. The 
orthonormality condition 

10,) =(5,,, (1.100) 

holds if the correct normalization is applied. 

Let us denote the dimension of the physical subspace N. Then because A is a physical 
operator, there exists an index set X = {^1,^25 ■ ■ ■ jJn} such that C = {\<Pj)}jfzj is a basis 
for the physical subspace. From the orthonormality condition in eq. 11.1001 it follows 
that C_L = is a basis for the orthogonal complement of the physical subspace. 

Because the orthogonal complement of the physical subspace is also the null space of the 
operator TZ, we conclude that is a member of the basis C for the physical subspace if 
and only if TZ {i/jj) 7^ 0. 

1.10 Looking Back 

Now that the generalized Dyson boson-fermion mapping is derived, we are in a position 
to ask a few questions and put the ideas presented into a larger perspective. For instance, 

• To what end was this theory first developed? 

• Where can it be applied, and what do we gain by applying it? 

• To what extent is the presentation of the theory contained in this chapter similar 
to presentations that have appeared in the past? 

• How complete is the presentation? Are there aspects that are not yet fully under- 
stood? Are there unproved conjectures on which the theory relies? 

As was mentioned in the introduction, generalized Dyson boson-fermion mappings are 
instances of what is known as boson expansion methods. Workers in this field will probably 
explain the utility of boson expansion methods as follows: A generic boson expansion 
will map fermion states onto boson states. (Before the development of the generalized 
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Dyson boson-fermion mapping, attempts have been made to include particles that are 
not bosons in the mapped system. However, these 'particles' were neither bosons nor 
fermions, but had more complicated commutation relations amongst themselves or with 
the bosons. See for instance jH p. 428].) A boson expansion is usually such that that 
a bi-fermion operator (^'^)^ in the original system maps onto a function of a few types 
of boson creation and annihilation operators. This function is typically a Taylor series 
with infinitely many non-zero terms. The order of a term in the series is equal to the 
number of normal ordered boson-operators it contains. The disadvantage of this (which 
partially prompted the study of Dyson-type mappings) is that for matrix elements of 
mapped operators between some state vectors, the series converge very slowly. Among 
the advantages are that, whereas in the original system there were say Q different one- 
particle states that a fermion could occupy, the mapped system consists of M coupled 
bosonic oscillators, with typically M '^Q. Also, for some states or regions of the spectrum 
of the system, high-order terms in the expansions of operators have a negligible effect and 
may be forgotten. Thus one is sometimes able to approximate a complicated fermion 
system with many one-particle states, with a system consisting of a few coupled bosonic 
oscillators. In the mapped system, the truncation of high order terms often result in 
mapped Hamiltonians with no more than, say, four-body interactions. One might be a 
little skeptical about the claim that the approximate mapped bosonic system is simpler 
to work with than the original fermion system. To my mind, it is not clear that an 
interacting boson system, with M one-particle states and on the order of bosons to 
occupy these states, is necessarily simpler than a fermion system with on the order of 2N 
fermions that are arranged among Q one-particle states, even when Q is much larger than 
M. Presumably, the simplification arises as follows: When one throws away high order 
terms in an expansion, one is effectively taking a limit and reducing the description of 
the system to include only what is pertinent to that limit. In the end, I think one has 
to take the word of authors who found it possible to extract useful information about 
systems using this technique, that the approximate mapped systems are simpler to work 
with than the original ones. 

A generalized Dyson boson-fermion mapping, on the other hand, results in images of 
the fermion operators that are functions of boson operators, with only a few non-zero 
terms. Our instinct to throw away high order terms in power series expansions is of no 
use.^'' In other words, a Dyson-type mapping does not wear on its sleeve a sign that 
suggests how to proceed after the mapping is done. Perhaps the reader is puzzled by the 
negative tone of the last statement: After all, Dyson type mappings give us boson-fermion 
systems with relatively simple Hamiltonians. As is the case for other boson expansion 
methods, fermion systems with 2N fermions arranged among Q one-particle states are 
^"Sometimes this is a good thing. 
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replaced by boson-fermion systems where the states we are typically interested in, are 
mostly bosonic. There are of the order bosons, arranged among M <^ Q states. It 
seems we have all the advantages of the other mappings, but without having to throw 
away terms. In my experience though, a Dyson-mapped system is as non-trivial to analyse 
as the original. It is just that some techniques of analysis are better suited to the mapped 
system than the original. Other analytical tools might be easier to use if one sticks to the 
original system. I would venture to say that, whereas other mapping techniques separate 
out a certain limit of the problem and introduce simplification by disregarding details 
unimportant to that limit, Dyson-type mappings simply constitute a way of looking from 
another angle on the original problem, with all its complicated detail still in tact. Why 
then have people developed the mapping? I think there are four reasons: 

Firstly, Dyson-type mappings can be used as a stepping stone to deriving many of the 
other boson mappings that comprise the field of boson expansion methods jE] . Secondly, 
there was a need to accommodate systems with an odd number of particles. To many 
workers in the field, the idea seems natural that a bosonic description should be useful for 
fermion pairs. However, if there are an odd number of fermions in a system where bound 
fermion pairs are energetically favoured, one would guess the low energy states to consist 
of configurations with all but one of the fermions bound in collective pairs. Ideally, in the 
mapped system, one would like the one unpaired particle to be represented as a fermion 
independent of the bosons. A Dyson-type mapping has proved to be the only mapping 
'simple' enough to accommodate such fermions in the mapped system. At reading the 
previous sentence, the reader might have felt a little uneasy. I previously mentioned that, 
in my experience, the mapped system is as non-trivial to analyse as the original system 
it derives from. So, if it does not lead to an obvious simplification, what does it help to 
have a mapping, even if it can accommodate fermions on the mapped side? The answer 
to this question is the third reason: To build nuclear models one can follow one of two 
approaches. The first is a microscopic approach. Here we start with a number of protons 
and neutrons, whose properties and interactions we think we know something about. We 
then write down a Hamiltonian that is based on this knowledge about the properties of 
the nucleons, and check if indeed nucleons with these properties result in nuclei with the 
properties we measure in the laboratory. The other approach is phenomenological. Here 
we make very few assumptions about the properties of nucleons themselves but focus 
on the known properties of the nucleus as a whole. Our goal is to write down (almost) 
any Hamiltonian that reproduces (if interpreted sympathetically), the properties of the 
nucleus. If we are successful in this task, we try to deduce properties of the nucleons 
from the form of the Hamiltonian. It so happens that phenomenological models where 
the quanta are bosons, or for systems with an odd number of particles, bosons and one 
fermion, reproduce observed nuclear spectra very well. To understand what these models 
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say about the fermionic nucleons, and to check that the reproduction of spectra is not 
simply a coincidence, we would like to construct boson-fermion mappings so that we can 
make a translation between the microscopic models and the phenomenological models. 
For this translation it would be nice if the same mapping works for systems with an 
even number of particles and systems with an odd number of particles. After all, on the 
(microscopic) fermion side, we believe that the stationary states of even and odd nuclei 
are eigenstates with different numbers of particles, of the same Hamiltonian. On the 
phenomenological boson side we give expression to this notion by also writing down one 
Hamiltonian that treats both the even and odd case. A single mapping should link the 
two sides. Also, since the phenomenological models contain only few-body interactions, 
it should be investigated whether a mapping is possible where a correspondence can 
be established without infinite series that have to be truncated. This wish list is very 
compatible with Dyson-type mappings. The one problem is the non-unitarity of Dyson 
mappings. Sometimes this can be overcome by a similarity transformation that preserves 
the desirable properties of the generalized Dyson mapping, but still transforms the mapped 
Hamiltonian into an hermitian operator The last reason that generated interest in 
Dyson-type mappings is this: In the calculation of matrix elements of operators, it is 
more pleasant to use the commutation rules of bosons than those of the more complicated 
fermion superalgebra. A mapping such as the generalized Dyson mapping is useful in this 
context, because it allows us to work with bosons and their commutators, without making 
operators so complicated under the mapping that the advantage of working with bosons 
is lost. For mappings in which the images of the original operators are infinite series of 
boson operators, we might have to calculate vacuum expectation values of products of 
arbitrarily large numbers of boson operators. 

If we remove the search for a link between microscopic and phenomenological models 
from this list of reasons, the motivation for studying Dyson-type mappings becomes a 
little thread-bare. One aim of the this thesis then, is to find new situations in which 
Dyson-type mappings can be useful. Also, in the type of application mentioned above, 
systems are considered that have Hamiltonians that can be written in terms of only the 
operators (^'^)^ and (^4^)^ of the even sector of the superalgebra. In such systems, the 
mapping is then usually only carried out for the collective subspace, so that the mapped 
systems contains at most one ideal fermion. Another aim of this thesis is then to push 
the generalized Dyson boson-fermion mapping beyond these restrictions. 

I conclude by making a few remarks about how the development of the theory, as 
presented in this chapter, compares with the original presentation published in [7j-|llj. 
As mentioned in the introduction, the present text contains more detail than the original 
publications. Therefore, for instance, when I derive the properties of the superalgebra, I 
use a general representation, to be doubly sure that the original fermion representation and 
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the 'induced' boson-fermion representation share these properties. This is in contrast to 
what is done in ^T]. Also, I take care to show a way in which meaning can be given to the 
demand that, in the definition of the Usui operator, the ideal and real fermion operators 
(or more precisely a and a) anti-commute. Another example is the demonstration of the 
invertibility of the mapping. However, these are small matters. A person, like myself, 
who is new to the field, might struggle with them for a while, but the solutions are simple 
enough that after some thought anyone can see the light. A further general remark is that 
I shied away from constructions that rely heavily on group-representation theory. Rather, 
I stuck to more or less elementary linear algebra of finite dimensional vector spaces. Were 
I more proficient in representation theory, I might have been able to derive generalized 
Dyson boson-fermion mappings rigorously from the (super)-coherent state representation 
of the states of fermion systems, in a way that is analogous to the construction in jH]. 
The disadvantage of not taking the coherent state route is that the definition of the Usui 
operator seems mysterious, as if plucked from thin air. On the other hand, I find the 
directness of the pedestrian approach appeahng. 

The one result in this chapter that I do consider a new contribution of some importance 
is the fact that the similarity transform X of Section 11.71 is not necessarily well-defined, 
but that trouble is avoided if we work inside the collective subspace. One may ask how 
serious a problem this definition of X is. Could it be that X is well-defined, but the proof 
of this fact eludes us? In this case, one may gain confidence in the similarity transform by 
constructing concrete mappings that use X, checking that in these specific situations, no 
problems occur. On the other hand, if X is sometimes not well-defined, we might hope 
to demonstrate this by an example as well. In Chapter Three we will do just that. We 
will look at a very special case where the superalgebra is such that the action of X and 
X~^ can be probed carefully. We will see that, for the wrong number of particles in the 
system, X is indeed not well-defined. Therefore, the problem with A is a serious one. 
If we are particularly pessimistic, we might be troubled by the fact that the example we 
consider is so special and ask if the situation is not worse for more generic examples. Does 
our example where the similarity transform is well-defined represent the exception or the 
rule? This is an open question. Because of this state of affairs, I think it is very important 
to re-emphasize that, if we stay in the collective subspace, the mapping is always valid 
and does not depend on whether X is well-defined. 

Appendix A 

In this appendix I derive the result stated in Section 11.31 (Symbols used here have the 
meaning assigned to them in that section.): For any representation in 7^5(C, (7, x, X, A), 
let be any operator in and let Ol be arbitrary elements of 0\ then there exist 
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Ol j G and complex coefficients wj sucli tliat 



;i.ioi) 



\k=l 



The Wj can be cliosen such tliat tliey are valid for all representations in 7ls{C, C, x, X, A). 
Consider a vector 



N \ 

k=l / 



;i.io2) 



Since all the elements of commute or anti-commute, we can, without loss of generality, 
take this vector to be of the form 



N2 



Two cases must now be distinguished: 

1. = di, and 

2. 6^ = Ai. 

We start with the ffist case. It is clear that a vector of the form 



(1.103) 



can be written as 



\i=i 



;i.io4) 



;i.io5) 



by anti-commuting to the right and using the fact that \i) = 0. Now we are ready 
to compute 



Ni 



N2 



^1=1 
' N2 



N2 



' N2 



\i=i 



N2 



;i.io6) 



By anti-commuting the fljy to the right in the second term, we can cast it in the desired 
form, as mentioned above, so that we only have to worry about the ffist term. Therefore, 
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consider tlie commutator tliat appears in the first term: 



iVi 



j=i \i=i 



ay , A"' 



n 



ExjM n ^ 



;i.l07) 



This result, substituted into the first term of eq. 11.1061 results in a sum over terms of 
the desired form and so completes the proof for the case where = a^. Notice that 
the complex coefficients uj were calculated without any knowledge of the representation, 
other than that it belongs to 7ls{C,C,x^X^ '^): and are therefore equally valid for all 
representations in 7ls{C, C, x, '^)- 
Now, for the case where = A,-: 



No 



N2 



:i.io8) 



It is not hard to compute the commutator appearing in the first term. The result is 

( Nx _ 

A, , n 



i=i 



Ni 



Ni 



E n ^M^r- E c^ '^i n 

j=l I \l=l7^j J l=j+l \m=l^j,l 



;i.io9) 



If we substitute eq. 11.1091 into eq. 11.1081 we see that we still need to deal with the the 



Kj^' operators. The following result takes care of that: 



A'2 



,1=1 



N2 



N2 

n 

,1=1 



N2 



U=l 



N2 

n ' 

.m=l^l 



N2 



;i.iio) 



If we substitute this into eq. 11.1091 we are done with the first term in eq. 11.1081 Let us 
turn our attention to the second term of eq. 11.1081 We have to compute 



N2 



\i=i 



N2 



N2 



;i.iii) 



1=1 



\m=l 



Appendix B 



43 



We can now anti-commute the to the right, thereby manipulating the second term of 
eq. 11.1081 into the form we want. This completes the proof for the second case where 
= Ai. Again, note that the numerical coefficients appearing in the derivation are 
calculated without any knowledge of the particular representation and can be expressed 
in terms of the quantities C, C, x, x ^-iid A which are common to all representations in 
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Here I derive the mapping of eq. 11.371 reproduced here for convenience: 

{a'')p < — {a'')j^p, = a'" + xrB'af,, 

(^)^ — {Ki)^^ = iq-CilB'Bi. (1.112) 

As extensive use will be made of the BCH formulas, they too are reproduced here for easy 
reference: 

exp(P)g = (^£l(P,Q),^exp(P), (1.113) 

Qexp(P) = exp(P)(^^^(P,Q),j. (1.114) 

where (P, Q)q = Q and (P, Q)k+i = [P , {P, Q)k]- As has become customary in the text, 
we take {ip) to be an arbitrary vector in Hbf and \(f)) to be an arbitrary vector in Hp- 
We set 10) =T\(f)). 

We start with the image for {a^)p. This is one of the instances mentioned in the text 
where the first step involves guessing the image, and then testing it afterwards. Taking 
{ay)gp = ay OS our guess, we may write 

{ilj\ay\4>) = a^T , 

= {(^|aj®(0|f/|0)® !</.), 

= {ij\0{O\aM\O)®\(t)). (1.115) 

In the last line of eq. \l.ll^\ a^,(I can now be calculated, using the BCH formula eq. 11.1141 
The result is 

aM = U{ay-{c^p}. (1.116) 
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If we substitute this into 11.115^ the ai, term falls away because annihilates the vacuum 
|0), so that we have 

(V'la.l^) = (V'|®(0|[/S;|0)® 10), 

= (^|®(O|/7|O)®{(a.)^|0)} 

= (^|T(a,)^|0), (1.117) 

which proves that (oj,)^ < — a^,. The steps for deriving the image of {Ai)p are exactly 
the same, this time starting with a guess Bi as the image. 

Next we derive the image for (a^)^. The process in this case runs forward, but the 
calculations are more involved. We start by writing 

(^1 {a^F \<P) = m ® (0| Ui^j, |0) ® 10) . (1.118) 

The BCH formula eq. ll.llHI is now employed to calculate 

UiaJ^p = \^{aP)p + xrB^{^p +^ U. (1.119) 

When this is substituted back into eq. ll.llHI term (1) drops away because {d')p annihilates 
the functional (0|. We turn our attention to term (2). Since (a^)^ commutes with U we 
have 

(^|®(O|xr5'^(a3'i.^|O)®|0) = (V^|®(O|xr5'^f/(«3FlO)®l0), 

= (^|xri?^r(a;.)^|0), 

= (v^ixr^^v (1-120) 

where we used the image T{a^)p = a^T which was proved above. We conclude that 

Now for deriving the image of {A')p, which proves to be the most taxing. Again we 
start by writing 

(^1 {A^)^ \4>) = (V^l ® (0| U{Af)p |0) ® 10) . (1.121) 
Now we employ BCH formula eq. II. 1131 to find 



Ai , A^ 



- '^CilB^B'^Ai - xrxixB^c^'K + xU^^c^" 

V ' ' V ' 

4 5 



U. (1.122) 



Upon back substitution into eq. 11.1211 terms (1) and (3) disappear because both {A^)p 
and {a'^)p annihilate the functional (0|. Since (0| [(A)i? , {A^)f\ ~ 9^1 (^1 replace 
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term (2) with gSlB^. In term (4) we can replace {Ai)p with Bi by the same argument 
that allowed us to replace (a^)^ with in eq. I1.12()[ Similarly in term (5) we replace 
{au)p with 5^. Terms (2) and (5) in their altered form may now be combined and the 
factored out to find 

B' {g5l - x'p^a^^^p) = B^. (1.123) 

Term (6) is simply equal to A> . Having said all this we conclude {A^)p < — B^lCl — 
^CflB^B'^Bi + . It is now straight forward to calculate the image of [kI^ ^ from the 
commutator of the images of {Ai)p and [A^)p. 
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In this appendix, I check that the images in eq. 11.371 obey the commutator and anti- 
commutator relations that they should, if they are the generators of a representation 
of iS(C, C, X, X) ill the whole TisF- To evaluate most of the commutators and all of the 
anti-commutators involves only the mechanical use of standard manipulations. I therefore 
concentrate on the following three commutators, the calculation of which I deem to be less 



trivial: 



{A 



qf(A')^^, [{A^)^p,{A^)^p] = and [(A^^^ , (a,) 



i^JBFi 



'BF- 



For the first two of these, the whole calculation is done explicitly, for the last, 
the essential hint necessary to do the calculation is given. We have to use some properties 
of the structure constants C and x that were derived earlier, and are reproduced below 
for convenience: 



^ki — ^ik 

2^ikM 1 

^ fip j Xrj I 
^Xi X.Xp/(.k X,pr]' 



;i.i24) 
;i.i25) 
;i.i26) 

;i.i27) 



Before we start evaluating commutators, it is necessary to derive an identity that the Cjl 
obey, namely that 

(1.128) 



'-^ik'^ln 



^nk^li ■ 



The proof is as follows: Consider the left-hand side CUC]^ of the equation and shuffle 
around the indices of the respective C's, using eq. II. 1251 to find CUC]^ = CUCfJ^. Now 
substitute from eq. ll.l27l for CfJ^^ to find CUCf^ = ^Clix'i'^X^isXZ^XaS- In this expression, 
isolate Cllxf and substitute from eq. |1.126| to find C^Q™^ = ^xrxipXf^X^^pXnXaS- 
Reshuffle the x's to find 



^ik^ln 



\ 



9 



m OLp J 



(1.129) 
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From eq. 11.1261 it follows that the first factor is equal to C^l^xf'^, which when combined 
with the second factor, and with the aid of eq. II. 1271 gives CUC^ = C^^C]^^ as promised. 
Let us now calculate 



J BF \ 1 J BF\ 

cfi {A' + 5^4) - \ [c^tci: + ci'^ci: - c^rcnj) b-b^'b^. 



;i.l30) 



To get to the last line only the familiar commutator identities between boson operators, 
and the commutator identities associated with S{C, C, x, X) ^''^^ required. We now apply 
identity eq. 028l to the last line of eq. 0801 The term Cfr^uv is equal to C^tQ^ so 
that cancellation between the second and third terms takes place. Only Cj^Cf^ remains 



and this is equal to C^C^^. Thus we find 



(A') 



BF ' \ ^ J BF. 



Cfi (a' + BPK}.p - ICZB'^B^B, 



1^ \ / BF 



'i.isr 



Now for the commutator [{A^)^p , {A^)^p]: if we substitute the explicit forms of {A^)^p 
and {A^)^p and distribute the commutator over the various terms thus generated we find 
that we have to compute six nontrivial commutators: 



{A') 



BF \ J BF. 



A'^ , A^p 



B'B^B, , B^B' B. 



1 

- \CilB'B^ [Bi , B^] fCl + B' 



fCj^ , A!^ 



3 4 

ICl^B^B' [B., , B-'] jq + B'B^ \}C{ , JC^ 
' ' ^ . : 

5 6 



;i.i32) 



It is not hard to check that terms (1) and (4) cancel, and that term (2) is zero all by 
itself. Furthermore 



term (3) + term (5) = - (C^/Cf - 0^10) B^B\ 
For calculating term (6), recall that we know 



;i.i33) 



;i.i34) 



(see eq. II. 4|) . But we know that fCl , /C^ 



This implies another expression 



for 



namely 



;i.i35) 



Appendix C 



47 



This means that there are two equivalent ways of writing term (6): 

term (6) = (C^liq - Cl^K^] B^B\ (1.136) 



and also 

term (6) = (C^^/CL - C^lCl) B^BK (1.137) 
If we take half of the first expression plus half the second then 

term (6) = \ (cgiq - CglCl) B^B\ (1.138) 

Thus term (6) cancels terms (3) and (5), so that (^*)_bf ^^'^ '^'^ indeed commute. 

I do not consider it necessary to derive in such detail the commutator 



BF 



1 \"'i^JBF 



The result follows easily when one uses eq. ll.l26l to convert a factor of the form Cx into 
one of the form XXX- 



Chapter 2 



A Toy Application: 
The S0{5) Model 

2.1 Introduction 

In this chapter we apply the theory of Chapter One to a specific example. By this we 
mean that the abstract symbols of the previous chapter are replaced with actual numbers, 
vectors and operators. In doing this, our aim is to show that the mapping procedure we 
advocate is fairly straight forward and that it is feasible to analyse the mapped system 
along the same lines as we are used to doing for any quantum system. We will also be 
afforded the opportunity to compare the original system to its mapped counterpart. This 
will firstly indicate the correctness of the theory, as for instance the physical eigenvalues 
of operators will be seen to agree in the two systems. Secondly, we will develop some 
understanding of what a generalized Dyson mapping does to a system, by examining how 
agreement between the two systems emerges. 

The model to be studied will be properly introduced shortly. Its most important 
feature for our purposes is that it may be expressed in terms of products of the generators 
of the algebra 5*0(5). Although we consider it a toy model, and exploit it only for its 
non-trivial but simple mathematical structure, it has frequently featured in the literature, 
sometimes with a view to understand very 'realistic' physics. I mention a few instances 
here to underline the fact that, despite our opting for the tutorial advantage of simplicity, 
we are within range of the coal-face where actual research is currently being done: For the 
early history of the model, which involved several different applications in nuclear physics, 
a brief discussion and many references may be found in |31 p. 411]. More recently, an 
S0{5) model was analysed to understand the role of proton-neutron pairing and four- 
particle correlations in nuclear wave-functions ^Hj - [T^ . 

Outside the field of nuclear physics, 5*0(5) models are reported to have some relevance 
in non-perturbative Quantum Chromo dynamics |2I], [21]] • In the field of High Temper- 
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ature Superconductivity, there is also significant interest in 5*0 (5) models fHl- For this 
last application, a substantial review has recently appeared [TH]. It is interesting to note 
that in all these applications, boson expansion methods also feature. Especially in more 
recent applications, the expansion method of choice is the generalized Dyson mapping. 
However, to be fair, I must add that boson expansion methods, generalized Dyson map- 
pings and the results obtained with their aid do not occupy a central position in any of 
the fields mentioned. 

This chapter is structured as follows: First we introduce the algebraic building blocks 
of the fermion system to be studied. We then do a mapping to a boson-fermion system 
for the collective subspace of the original fermion system. We explicitly include the 
possibility of having a single unpaired fermion in the mapped system. This extends the 
type of analysis that is done in for instance ^3] , where only the even part of the collective 
subspace is mapped, so that there are only bosons in the mapped system. We do not, 
however, go beyond the inclusion of a single fermion for the mapped system. This would 
have entailed working with the similarity transform of Section [TT7I We rather defer doing 
this until Chapter Three, where we will work with an even simpler algebraic structure 
and the similarity transform can be implemented more transparently. With the mapping 
established, a specific Hamiltonian is introduced. This Hamiltonian is interpreted as a 
nuclear Hamiltonian in for instance [H p. 411] and references therein, which explains 
the terminology we adopt. After the mapped counterpart of the Hamiltonian is written 
down, both the original and the mapped Hamiltonians are diagonalized. It is shown that 
following what seems to be the simplest procedure for diagonalizatoin in each system, 
leads to virtually identical calculations in both systems. This is perhaps not yet fully 
appreciated in the literature, as is briefiy discussed. After diagonalization, ghost states 
are identified, and aspects of their behaviour, as a function of the parameters of the model, 
are analysed. 



2.2 The Algebra and the Mapping 

We consider a one-particle Hilbert space for which an orthonormal basis can uniquely 
be specified using two quantum numbers m G + and a G { — ^il}- 

Here j is an odd multiple of a half. We define an integer Q = ^ + j- The one-particle 
Hilbert space therefore has dimension 4f2. With these quantum numbers we associate 
anti-commuting fermion creation and annihilation operators and ama that operate in 
a fermion Fock space with vacuum |0). We will say that the operators with cr = | create 
and annihilate protons, while the cr = — ^ operators create and annihilate neutrons. With 
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the single-fermion operators as building blocks, we construct the following operators: 

Ei±i = ±t (2-1) 



-2 "'^2 



Elo = E (2-2) 



-J 



1 ^ 

^01 = ^ E (2-3) 

and their hermitian conjugates £"-1^:1 = -£"1^1, -E-10 = and -Eq-i = -E'oi- Furthermore 
we define two hermitian operators 

1 

1 ^ 1 

= ^A^-^, (2.4) 

and 

1 

2 J 



^2 = E E 
1 



. i m= 
2 



2 (N, - iY„) . (2.5) 

In these last expressions, Np counts the number of protons, while A'^ counts the number 
of neutrons and of course = Np + Nn counts the total number of particles. The operator 
Ell creates collective proton-proton pairs; Eiq creates collective proton-neutron pairs and 
creates collective neutron-neutron pairs. The operator Eqi changes one neutron into 
a proton and i?o-i does the reverse, namely it changes one proton into a neutron. From 
this the rule emerges that an operator i?Q^Q2 changes the particle number of the system 
by 2ai and the charge by ai + a2 times the charge of a proton. 

The above ten bifermion operators form a basis for a Lie algebra under commutation. 
In fact, they represent the so-called Cartan-Killing basis for the algebra 5*0(5). This 
sounds rather intimidating but means nothing more than that a simple recipe exists for 
writing down the commutator of any two of the operators. To explain this recipe, we 
take up the two indices ai and a2 that every E operator is subscripted by, into a two 
dimensional vector a = (ai,a2). There are eight possible vectors a, namely (1, 1), (1,0), 
(1,-1), (0,1), (0,-1), (-1,1), (-1,0) and (-1,-1). These are called the roots of the 
algebra. The commutators between the bifermion operators can now be listed. If we 
commute an H and an E we get 



[Hk , Ea] = CUkEa, 



(2.6) 
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so that for instance [H2 , i^o-i] = — -E'o-i- If we commute and E^a the resuh is 

[E^,E.o] = E (2-7) 

k=l,2 

So for instance we have , = —Hi + if2. The operators Hi and commute 

with each other. Lastly, if we commute Ea with Ejj then the resuh is zero unless a + (3 = 
{ai + Pi, a2 + P2) is a root of the algebra. If a + is a root, the commutator is proportional 
to the Ea+/3 operator. This result can be stated as 

[Ea , Ep] = Ma,pEa+/3, (2.8) 

where M^^is is a complex number, which is zero in cases where a + (3 is not a root of the 
algebra. So for instance [En , i?o-i] = — -E-io since (1 + 0, 1 — 1) = (1, 0) is a root of the 
algebra, whereas [i^i-i , -E-i-i] = because (1 — 1, —1 — 1) = (0, —2) is not a root of the 
algebra. The constants Ma^p have the following symmetry properties: 

Ma,p = -Mp^a = -M_a,-0 = Mf3^_a-/3 = M_a-/3,a. (2.9) 
The value of Ma^p can be read off from the table: 



a 



Ma,p 


(1,1) 


(1,0) 


(1,-1) 


(0,1) 


(0,-1) 


(-1,1) 


(-1,0) 


(-1,-1) 


(1,1) 














1 





-1 




(1,0) 











1 


-1 


1 




-1 


(1,-1) 











-1 







1 





(0,1) 





-1 


1 










-1 


1 


(0,-1) 


-1 


1 










-1 


1 





(-1,1) 





-1 







1 











(-1,0) 


1 




-1 


1 


-1 











(-1,-1) 




1 





-1 















(2.10) 

There is an hermitian operator C, called the second Casimir invariant, that commutes 
with all the E and H operators that constitute the Cartan-Killing basis. This operator 
is given by 

C = H^i+ Hi + {^n , E^i^i} + {Eio , ^_io} + {^i-i , ^-ii} + {^oi , ^o-i} , (2.11) 

where curly brackets indicate anti-commutators. Using the commutator identities for the 
generators of the SO{b) algebra, we will later on manipulate this expression for C into 
other, more convenient forms. 



Section 2.2: The Algebra and the Mapping 



52 



The next step is to decide on a subspace of the fermion Fock space on which our 
operators are defined, to use as the domain of our mapping. As can be seen from the results 
in P, expressions involving the similarity transformation of Section [TT7I and the images of 
the single-fermion operators and ama are unpleasant to work with. We therefore stick 
to the collective subspace where the complications arising from the similarity transform 
are avoided. This is the space spanned by states of the form 

{Enr{E,or{E,^,r~^\o), (2.12) 

{Ei^riEioriEi^ir-'aiJO). (2.13) 

All the Ea and Hk operators leave this space invariant as was discussed in Section EHl and 
the general formulas of that section are applicable. We will set up our notation as follows. 
Let = -^Eii, A'^ = -^EiQ and A~^ = and let A^ be the hermitian conjugate of 
A'' so that for instance A_i = -^E^n. It is this set of operators A'' and A^ that we use 
in the definition of the Usui operator in Section Fl.TI ^ The structure constants C*f such 



that 



A' 



A, , A' 



C^^A^ can be derived from commutator identities for elements of 



the Cartan-Killing basis. For the non-zero structure constants we find 



'-'10 ~ '-"01 ~ '-'10 ~ '-'01 ~ '-"00 ~ '-"-10 — '-'0-1 — '-"-10 — n 

r*!-! — r'OO _ -^00 _ ^-11 _ 1 (n ^A\ 

•-"00 — '-'1-1 — '-'-11 — '-'00 — ~fi- l-^'-L^J 

Every structure constant not present in these expressions, is zero. With these structure 
constants the essential part of the mapping reads 



,^11 
'-'11 - 


czlz 


1 2 

1 - n 


r'Oi _ 

'-'01 — 


'-'00 - 


rfO-i 
■ '-'-10 


'-'1-1 - 


r'OO 
'-'-11 


— '-'00 



A' ^ 


- B'- 


A' ^ 


- B^- 


A-' ^ 


- B-' 


Ai ^ 


- Si 


Ao - 


- Bo 


A-i - 


- B_, 



i - B' (No + N,+^^p) - B% 

U-B'B-'Bo -B^(n^, + N, + ^ 



Mp + AC, 



+B^^Sqi — B^£o_i^ 
i - B~' {No + N^i+Afn) - B'So 



(2.15) 

The operators A/'p and A/'n count the numbers of ideal protons and neutrons respectively, 
while Sqi and Sq^i are the ideal fermion counterparts of Eqi and i^o-i- The operator A''^, 
counts bosons of species k G {—1,0, 1}. The images of the operators Eqi, Eq-i, Hi and 
^The factors of in the definitions of the A'^ will be explained after the mapping is written down. 
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H2 can either be obtained as the commutators of the appropriate images in eq. 12.151 or 
by using the general equations of Section ITTHl For Eqi and i?o-i we find 

Eoi < — B^Bo — B^B^i + Soi 
Eo-i ^ B'B,-B-'Bo + So-i. (2.16) 

Instead of writing down the images for Hi and H2, 1 write down the images for the proton 
number operator Np and the neutron number operator Nn, that are linear combinations 
of Hi and H2 (and the identity operator). The result is 

Np < — 2Ni + No+N'p 

Nn ^ 2iV_i + A^o+A4. (2.17) 

The reason for the factor of ^ in the definition of is this. Say we had defined = Ei^ 
instead, and used this definition in the Usui operator. Then the images of for instance En 
and would have been QB^ + - B^ {Nq + Ni + Mp) - B^£m and Bi. The 

state Ell |0) would have mapped onto fiS^ |0) and the functional (0| -E-i-i would have 
had an image (0| This asymmetry, namely that the images of pair creation operators 
stretch states by a factor Vt while functionals undergo no stretching of order VL if the images 
of pair annihilation operators are appended to them, is untidy. With the introduction 
of the factor ^ into the definition of the A'' used to construct the Usui operator, the 
images of En |0) and (0| respectively become v/fi (0| Bi and VUB^ |0) so that the 

fi-scaling of the images of bra's and kets are now symmetric, and the mapping 'only as 
non-unitary as is absolutely unavoidable'. 



2.3 The Model Hamiltonian 

To get a feel for the mapping, we now consider a toy model in both its original and 
mapped incarnations. As the terminology we have used up till now indicates, the model 
is ostensibly a nuclear model. The original Hamiltonian is given by 

H = ^EioE^io + g2EoiEQ_i. (2.18) 

There is no one-body term J2m,a '^tia'^ma so that, were the model to describe anything, 
it would be the residual interaction between nucleons in a single j-shell. Because of the 
absence of a one-body term, the Hamiltonian leaves the collective subspace invariant. We 
are allowed to restrict our analysis to the collective subspace so that the mapping we have 
just written down is applicable. 

Since there are only two independent energy-scales, the Hamiltonian can be charac- 
terised by one dimensionless parameter, which we will call 9. If we are interested in the 
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behaviour of the system at every point in parameter space only up to an arbitrary (possi- 
bly negative) scaling of the spectrum, we may learn all we want to know by parametrising 

91 = GlCOSy, 

92 = G2siny, (2.19) 

and considering the system for all 6 G [—1, 1). (Here Gi and G2 are convenient constants 
whose job is to normalize the two limits where the first or second term in the Hamiltonian 
dominates. We chose Gi equal to one divided by the maximum eigenvalue of ^E'lo-E-io 
and G2 equal to one divided by the maximum eigenvalue of EqiEq^i.) The range of 6 
includes both the limits where the first or second term dominates. Also, both the cases 
where the relative sign of gi and g2 is positive and negative are included. This is the 
real reason why I chose the present model, namely that we can understand the essential 
behavior of the spectrum by varying a single parameter 9. However, the Hamiltonian 
in eq. 12.181 is not such a special Hamiltonian, in the sense that a rather general class 
of Hamiltonians reduce to it. 1 briefly digress to explain this claim. For any nuclear 
Hamiltonian it is a reasonable approximation to insist that (a) the total number of protons 
and neutrons should be conserved separately and (h) the system is invariant under the 
simultaneous transformation of all protons into neutrons and all neutrons into protons. 
If we add to this the unrealistic assumption that low energy collective excitations are 
described by a Hamiltonian that is a quadratic function of the S0{^) generators of eq. 
12.11 -eq. 12. 5^ all allowed Hamiltonians are of the form 

H' = §^10^-10 + If (£^11^-1-1 + £^i-i£^-ii) + ^3 {^01 , £^0-1} + /(^i, H2, C), (2.20) 

where / is an arbitrary function. We can get rid of the i?o-i-E'oi term using the commutator 
[Eqi , i?o-i] = H2. Also, we can manipulate the expression eq. 12.111 for the Casimir 
operator C to find 

-E^ii-E'-i-i + Ei^iE^ii ^ ~ -^i^-^i ~ 3) — H2{H2 — 1)) — EiqE^iq — E'oi-E'O-i- (2-21) 
If we substitute this into the expression eq. 12.201 for H' we find 

H' = ^EiqE^iq + q^EqiEq^i + f'{Hi, H2,C), (2.22) 

where g'( and g2 are numbers (different from the original g'l and g!^) and /' is a function 
(different from /). The assumption that the Hamiltonian describes low-energy collective 
excitations implies that diagonalization takes place in the collective subspace. Further- 
more, since the numbers of protons and neutrons are fixed, diagonalization takes place 
in a subspace of the collective subspace with a specific number of protons and a specific 
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number of neutrons. This means tliat in the subspace where diagonalization takes place, 
the operators Hi and H2, that are linear combinations of the proton and neutron number 
operators, are proportional to identity. Furthermore, because the proton and neutron 
numbers are fixed, we either work in the subspace of the collective subspace that contains 
states with an odd number of particles or the subspace that contains an even number 
of particles. In each of these subspaces, the Casimir operator is also proportional to the 
identity operator. Thus the function f'{Hi, H2, C) is proportional to the identity operator 
in the subspace in which we diagonalize H', and we may therefore neglect it at the cost of 
shifting all energies by a constant. Along this route we then also arrive at the Hamiltonian 
of eq. ITTRl 

The image of the Hamiltonian H of eq. 12.181 under the mapping we employ is 

Hbf = glS^No-^{2N_l + 2Nl + No-l+^^p + ^^n)No 
+ -^B^B^^BqBq + —B^^BqSqi — —B^ BqSo^i^ 

+g2 {SoiSo^i + iVi (1 + iVo) + iVo (1 + A^-i) - B^B'^B.B, - B^B^B^B^^ 

+B^BoSo-i — B^^BqSqi — B^B^iSq^i + B^BiSoij ■ (2.23) 



2.4 Diagonalization 

The obvious thing to do with a Hamiltonian is to diagonalize it. We will consider diag- 
onalization for both the original and mapped incarnations of the system. We will deal 
with the case of an odd number of nucleons, as this is the more interesting case. An even 
number of nucleons can be treated with a boson mapping that does not introduce ideal 
fermions into the picture. This has already been done in the literature - ^Z|; ^H] - 

122! 

Our first task is to characterize the space in which we diagonalize to such an extent 
that we can label and calculate the necessary matrix elements of the Hamiltonian. On 
the fermion side the space in which we diagonalize is spanned by states of the form 

{Eur (Eior {Ei^iT-' |0) . (2.24) 

The Hamiltonian preserves the quantum number m so that we can decide on a fixed value 
for it, and not explicitly indicate it in our notation. If we consider a system with Up 
protons and n„ neutrons, the following constraints hold: 
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Up = 2ni + no + ^ + a, 
1 

Un = 2n_i + no + - - a, 

Up + Hn is odd. (2.25) 

It is not hard to see that these constraints are enough so that, if Up and n„ are given, then 
for each permitted value of no, the number of proton-neutron pairs, the corresponding 
values of n_i, rii and a are uniquely given by 



a = {up + uq) mod 2 
1 / 1 



1 

2' 

ni = -ynp-no---a 

n^i = ^(^rin-no-^ + aj . (2.26) 

The constraint that determines a simply says that if Up and no have the same parity, 
then the unpaired particle must be a neutron, whereas if they have different parity, the 
unpaired particle must be a proton. The permitted values for no are 

no G {0, 1, 2, ... , min(np, n„)} . (2.27) 

We therefore adopt the following notation. The space, in which diagonalization takes 
place on the fermion side, is characterized by specifying n^, n„ and m. For these specified 
values and for every permitted no we understand ni, n_i and a to take on the values 
given by eq. I2.26[ and let |no) represent the state 

The set v = {|?t-o)}™=o*'"""'* spans the space in which diagonalization takes place on the 
fermion side. The reason for the factor , ^ and choosing to work with the A 

^(n_i)!(no)!(ni)! 

operators instead of the E operators in the definition of |no) is that, with this choice, the 
set V tends to an orthonormal basis in the limit Q ^ np + Un- However, as Q becomes of 
the order , the set v becomes a linearly dependent set. 

For a fermion space spanned by the set v, the physical subspace under a boson-fermion 
mapping is a subspace of the boson-fermion space spanned by states of the form 

B-)"-(B»)""(B')'%L|0). (2.29) 



The same particle number constraints (eq. 12. 25^ as for the original fermion system apply. 
(This follows from the form of the images of the total proton and total neutron number 
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operators as given in eq. 12. 171 ) Therefore we choose a basis v = {|'T'o)}n'=o where 



and ni, n„i and a are given by eq. 12.26^ for the boson-fermion space in which we 
diagonalize the mapped system. The tilde in the notation |f?^) is meant to emphasize 
that the state |?t^) is not the image under mapping of the fermion state |no). The set 
V always forms an orthonormal basis for the space in which we diagonalize the mapped 
system. 

We can now conclude that ghost states come into play for the mapped system as 
soon as the set v that spans the fermion space becomes linearly dependent. When the 
elements of v are linearly independent, v is a. basis for the collective subspace. In this 
case the dimension of the physical subspace is equal to the number of elements in v. 
The dimension of the boson-fermion space spanned by which is the space in which the 
diagonalization of the mapped system takes place, is equal to the number of elements in v. 
The sets v and v have the same number of elements. Hence the dimension of the physical 
subspace is the same as the dimension of the space span{v) of which it is a subspace. This 
of course implies that the physical subspace is the whole span{v). As soon as the elements 
of V become linearly dependent though, the dimension of the fermion space becomes less 
than the dimension of the space span{v) in which we diagonalize the mapped system. The 
space span{v) in which we diagonalize the mapped system must then necessarily contain 
more than the physical subspace. 

The spaces in which the operators are to be diagonalized are now characterized and 
diagonalization can begin. We first diagonalize on the fermion side. Some authors consider 
this a more difficult task than diagonalizing the mapped counterpart in the boson-fermion 
space. They see the fact that the set v, that spans the fermion space, is not a set of 
mutually orthogonal states, as a stumbling block. It seems the impression exists that an 
orthonormal basis should first be found and that because this is a complicated task, the 
alternative i.e. doing a boson-fermion mapping, is to be recommended. However, the 
simple construction below shows that the complications that arise when one works with 
a non-orthonormal, possibly linearly dependent set of states for the original system, are 
equivalent to the complications associated with the non-hermiticity and the presence of 
ghost states in the mapped system. 

To explain the diagonalization procedure on the fermion side, we first look at the 
following general situation: Let V be a vector space with dimension M. Let v = 
be a set of states that span V. In general we allow the possibility that > M so that the 
vectors in v are linearly dependent. None the less, we want to express arbitrary vectors 
in terms of the elements in this set. We are asked to diagonalize a fully diagonalizable 
operator H : V — > V. To keep things simple, we will take it that H has a non- degenerate 




(2.30) 
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spectrum. Let w = {\k)}i^^^ be the eigenvectors of H, sucli tliat H \ k) = Xk\k). We do 
not know tlie eigenvalues A. Nor do we know how to express the eigenvectors in w in 
terms of the elements of v, or even what the dimension of V is. Our aim is to find an 
expansion for each eigenvector \k) in terms of the elements of v. 

The set w is then linearly independent and a basis for V. We assume that at our 
disposal is an X matrix h with entries ha^p such that for every element \a) of f it 
holds that H \ a) = Y,^=i ha,i3 Note that when H is hermitian (which we don't have 
to assume), the matrix h will in general not be a hermitian matrix, the exception being 
the case where the elements of v form an orthogonal basis for the space V. To keep the 
discussion as simple as possible we assume that the matrix h is fully diagonalizable and 
has a non-degenerate spectrum, which, it turns out, is consistent with the assumption 
that H is fully diagonalisable and has a non-degenerate spectrum.^ 

Now, label the left-eigenvectors of /i by 7' = 1, 2, . . . , and let cy^a be the component 
a of the left-eigenvector of h that is labeled by 7', i.e. 

N 

Cy,aha,l3 = X'yCy,i3, (2.31) 

0=1 

where Ay is the eigenvalue associated with the left-eigenvector labeled by 7'. Now set 
1 7') = J2a=i cy,a I a)- Then 

TV 

a=l 

N N 

= Cy,„/l„,/3 1/5) 

(3=1 a=l 
N 

= a;, E 1/5) 

/3=1 

= A'yiy). (2.32) 

This means that if I7') is not the zero vector then Ay is an eigenvalue of H (and hence 
equal to one of the A^) and I7') is an eigenvector of H. 

Two questions still need answers. Firstly, does the spectrum of h contain all the 
eigenvalues of H7 Secondly, if Ay = A^, are we guaranteed that I7') is the eigenvector of 
H corresponding to the eigenvalue A^, and not the zero- vector? 

To answer these questions (both affirmative), there is the following argument. Let 
da,k, Oi = 1,2, . . . , N, k = l,2,...,Mbe the expansion coefficients of the states |a) in 

^The simplifying assumptions we make here happen to be vahd for our 50(5) example. They may 
however be omitted to fit more general applications, at the cost of making the argument and results more 
intricate. 
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terms of the eigenbasis states \k): 

M 

\a) = J2d^,k\k). (2.33) 

k=l 

Note that, for any fixed value of k, there exists at least one index a such that da,k 
is unequal to zero. This follows from the assumption that the set {|a)}^^^ spans the 
space V for which {|A;)}^j^ is a (linearly independent) basis. Now construct the vector 
Sfcli ^kda,k \ k) and compute 

M M 

'^hda,k\k) = ^da,kH\k) 

k=l k=l 

= H\a) 

N 

(3=1 
M N 

= T.T.h»,pdfs,k\k). (2.34) 

k=l (3=1 

By the linear independence of the different vectors |A;), the equality between the left- and 
right-hand sides of the above equation must hold term for term, so that J2p=i ha,(3dp^k = 
Xkda,k- This means that all the eigenvalues Xk of H are indeed contained in the spectrum 
of h. Now, take an index k, associated with eigenvalue of H, and let 7' be the index that 
refers to the unique (thanks to the non-degeneracy of the spectra of H and h) eigenvalue 
Ay of h such that Ay = A^. Now consider the associated state I7'), which we wish to 
show is not the zero-vector: 

TV 

IV) = H cy,„ |a) 

a=l 
M N 

= X] X] C7',aC?a,H0 • (2-35) 
1=1 a=l 

We now claim that Y.a=i '^Y ,0^0,1 is zero if and only if / 7^ A;. To prove this, let /„, 
a = 1,2, ... be the components of any right-eigenvector /, excluding dk, of h with 
eigenvalue A. fa, a = 1,2, N could be equal to da,i, a = 1, 2, . . . , A^, for some fixed 
/ 7^ k, but it could also be the components of another right-eigenvector of h that we have 
not yet encountered, since there are only M right-eigenvectors denoted with d^s, while 
there are right-eigenvectors in total. In any event, the set of all possible choices for / 
spans an — 1 dimensional subspace of C^. The non-degeneracy of the spectrum of h 
implies that A 7^ A^. Now consider the expression Y.a,(3=i '^Y,aha,(3f(3- Using the fact that 
c^i is a left-eigenvector of h we conclude that the expression is equal to XkY.a=i'^'y',afa- 
However, using the fact that / is a right-eigenvector, we conclude that the expression is 
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also equal to A Y^a=i CY,afa- Since A 7^ it must hold Y^a=i CY,afa = 0. This means that 
the vector with components c*, ^, a = 1, 2, . . . , is orthogonal to all right-eigenvectors 
of h that have eigenvalues different from A^. If J2a=i '^y,ada,k is zero too, then this vector 
has to be the zero-vector, which it is not, since cy is a left-eigenvector of h. Therefore 
Y.a=i c-y',ada,k is uot equal to zero and Y.a=i c-y',ada,i is zero if and only if / 7^ /c. From 
eq. I2.34l then follows that I7') is proportional to |A;) with the proportionality constant not 
equal to zero. 

The procedure for diagonalizing H if we have the matrix h at our disposal is clear: 
We first find the left-eigenvectors of h. We label a left-eigenvector by an index 7', and 
let its components be denoted cy^a, ol = 1,2, . . . , N. The eigenvalue associated with 
it is denoted Ay. From the components cy o, we construct states I7') = Z)«=i cy,a 
The non-zero I7') comprise all the eigenvectors of H. If I7') is non-zero, and hence an 
eigenvector of H, then it has an eigenvalue Ay. 

For our 5*0(5) fermion Hamiltonian, a valid matrix h of expansion coefficients can be 
found by making use of the commutation rules that define the superalgebra of which the 
5*0(5) generators are the even sector and the single fermion operators are the odd sector, 
as well as the action of generators on the vacuum ket. Quite naturally, we choose to work 
with the set {|no)}",=?'"""^ defined in eq. 12.271 When the Hamiltonian H acts on one of 
these states, we find the following results: 

• If no has the same parity as Hp-. 




(2.36) 




(2.37) 
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• If no and Up have different parity: 

np + Un - no - 1\ . , 
no [I ^ ) |no) 

+ ^\/('^n - no) {no + 1) \no + 1) 

+ ^ - no - 1) (n„ - no) (rio + 1) (no + 2) jrio + 2) , 

(2.38) 

-^\/('^p - no + 1) (n„ - no + 2) no (no - 1) |no - 2) 

+ ^\/("'P - ^0 + 1)"'0 I'^o - 1) 
+ |y {np + n„ - 2no) + yj 1%) 

-^\/(^n - '^o) ('^o + 1) l^^o + 1) 

-^\/(% - no - 1) (n„ - no) (no + 1) (no + 2) |no + 2) . 

(2.39) 

From these expressions the matrix elements /ino,n^ such that H |no) = J2^^o'"""^ ^no,'"o l^o) 
holds, may be read off. We see that /i is a matrix whose only non-zero entries lie on the 
main diagonal, the two diagonals above the main diagonal and the two diagonals below the 
main diagonal. At this point we can use a computer to diagonalize the matrix h numeri- 
cally for given values of np and n„. From the left-eigenvectors thus found, the eigenstates 
of the system can be constructed in terms of the states |no), no = 0, 1, . . . ,min(np,n„) 
and the unphysical eigenvalues of h identified. 

However, before we do this, we will discuss the diagonalization of the mapped Hamil- 
tonian Hbf- In the mapped system we work with an orthonormal basis of states |no), 
no = 0, 1, . . . , min(np, n„), as defined in eq. 12.301 and the conventional diagonaliza- 
tion procedure is used. First we calculate the matrix h whose elements are /imo.no = 
(nol i^Bi;' |nio). Then, with every left-eigenvector of h with components fng such that 
A//no = Emo fmohmo,no holds, correspouds a right-eigenvector IJr) = J2no /no \no) of Hbf 
with eigenvalue Xf. Similarly, with every right-eigenvector of h with components such 
that Xggmo = Y.no^mo,no9no hoMs, correspouds a left-eigenvector Igi) = Eno5'nol^) 
Hbf with eigenvalue Xg. (When we say that Iql) is a left-eigenvector of Hbf "we mean 
that {qlI Hbf = Xg{gL\.) 

As confusion might arise, I point out the following: In the mapped system the numbers 
n„ and n^ do not refer to the numbers of ideal neutrons and protons. They refer, of course, 
to the numbers of neutrons and protons of the original system. 



—EioE^io |no) 



Eni El 



Ol-C/Q-l 



ho) 
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h 



We now simply use the rules for commuting bosons to calculate the matrix elements 
{uqI Hbf |"^o) of Hbf- The same distinction as for states of the original system 



mo ,no 



must be made between the case where uq and Up have the same parity (meaning the ideal 
fermion is a neutron) and the case where uq and Up have different parity (meaning the 
ideal fermion is a proton). In what follows below, we again separate the Hamiltonian 
into two parts. The image of the term ^EiqE^iq of the original Hamiltonian is denoted 
(^^i^io-E-io)^^ and similarly the image of E'oi-E'o-i by [EqiEq^i] ^p. The explicit forms 
of these images in terms of the boson and ideal fermion operators are given by respec- 
tively the contents of the first and second set of curly brackets in the expression for Hbf 
(eq. 1223 • We find the following: 



If no has the same parity as rip: 



^EioE^iQ I \no) 

\l J BF 



— \j{np-nQ + 2) (n„, - no + 1) {no - 1) 



no 



+ ^\/{n'n -no + l)no 



+no 1 



no - 1 

np + n„ - no - 1 



2Q 



no) , 



(2.40) 



{EoiEo^i)^p |no) 



(np - no + 2) (n„ - no + 1) no (no - 1) no - 2 
no - l) 



-^{un-no + l)no 



I TliQ . > I I — \ 

+ 1 y ('^P + - 2no) + — \ |no) 



+ 2 V ("p ^ "o) ('^o + ^) '^o + 1 



Up - no) {jin - no - 1) (no + 1) (no + 2) no + 2) . 



(2.41) 



If no and rip have different parity: 



^^^lo^-io) |?^o) = -h^^^np - no + 1) (n„, - no + 2) no (no - 1) 

\l J BF ZiZ * 



no 



n^^nn- no 



no - 1 
- 1 



Section 2.4: Diagonalization 



63 



(-E'oi^O-l)Bi? |''^o) 



^2 V - ^0 + 1) {rin - no + 2) no (no - 1) no - 2^ 



no 



n^ 



+ 1 TT ('^P + '^n - 2no) + — r Fo 



(2.43) 



Comparison of eq. 12.361 - eq. I2.39l wliere we found tlie entries of the matrix h of tlie original 
fermion system and eq. 12.401 - eq. 12.431 for the entries of matrix h of the mapped system 
reveal that the matrices are related by /i^ = h. (Note that the entries of the matrices 
are real so that is simply the transpose of h.) It is not hard to understand why this 
relationship exists between the two matrices. Instead of talking about the present 5*0(5) 
example, we keep the discussion general. The only restriction is that we work with the 
mapping of Section 11.81 for states in the collective subspace, and not with the mapping 
for a larger subspace of the fermion Fock space that is arrived at in Section ITTTI bv means 
of the similarity transformation. We use the notation 



|n) 



M 

n 

k=l 



■n-k 



-a- |0) , 



(2.44) 



where n refers to the combined indices rii, . . . ,nM and u for fermion states in the collective 
subspace. Here the are the collective fermion pair operators. (We confine ourselves 
to the odd part of the collective subspace. The same argument works for the even part.) 
We assume a Hamiltonian that is a function of the A'', Ak and = [Ai , A^] operators. 
For the mapped system we use the notation 



|n) 



fe=i 



nfc! 



-a"" |0) . 



In the original system the matrix h has matrix elements such that H |n) 
As H is hermitian, this can also be written 



in 



(2.45) 

(2.46) 



In calculating the numbers /i* ^/ we do not have to make use of any properties particular 
to the fermion representation of the superalgebra. Let us suppose that we only use the 
commutation and anticommutation relations, and the action of various generators on the 
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vacuum bra, which are properties possessed by the boson-fermion representation too, for 
calculating /i* „/. (This is what we did in eq. 12.361 - eq. 12.391 ) In eq. 12.461 we may then 
replace all operators by their images under the mapping and the equality will still hold. 
It is here that the key lies. Thanks to the fact that the image of Ak is Bk and the image 
of (0| ttu is (0| ctjy, the bra {n\ is replaced by (n|.^ Therefore we conclude that 



{?i\Hbf = Y.[^' Kw (2.47) 

Hbf 



For the mapped system, the matrix h was uniquely defined such that /i„,^„/ = \n' Mbf n 
from which then follows that hn,n' = h*^' n- Thus the matrix h that we work with on the 
fermion side is the hermitian conjugate of the matrix h that we work with on the boson 
side. 

Bearing this in mind, consider the following claim, quoted from the introduction of 



It has often been recognized and stated that boson mapping is not only rel- 
evant to discussions about the relationship between phenomenological boson 
models and shell model type fermion models, but also constitute an attrac- 
tive many-body formalism in its own right [. . .] The latter pronouncement is 
against the background that one might profit from the use of boson variables 
in the description of many-fermion problems through the potentially simpler 
algebraic structures associated with bosons. This philosophy has its roots in 
the very simple observation that bosons commute to c-numbers, while fermion 
pairs commute to a non-trivial operator. Correspondingly the construction of 
an orthonormal basis is in general then simpler for bosons than for fermion 
pairs. 

In the paper in which this claim is made, boson mappings of models based on the algebras 

of iS'0(12) and S'p(lO) are considered. In ^Hj, where the model considered is a special 

case of the 5*0(5) model of this chapter, it is again claimed that 'the boson mapping may 

provide a simple technique to solve the fermion problem' (footnote 3, p. 241). There 

are two aspects that need to be considered when one claims that it is easier to work 

with the mapped system than the original. (This discussion refers to situations when we 

work within the collective subspace on the fermion side.) The first concerns the algebraic 

manipulations needed to compute h on the boson side, as compared to those needed to 

compute h on the fermion side. The mapped Hamiltonian is a function of a few types of 

boson operators. The result is that the commutations that must be performed to calculate 

the matrix h on the boson side are easier to keep track of than the commutations needed 

■^If we venture outside the collective subspace and consider the bra (0| a^a^, the argument will break 
down because this bra does not have an image (0| a^a^. 
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to calculate h on the fermion side, because bosons commute to c-numbers. None-the-less, 
the structure of the commutation rules of bi-fermion-type algebras are, in my opinion, not 
vastly more complicated than those on the boson side. If the remarks made in ^3] and 
|23] purely refer to a preference for computing the matrix h on the boson side, rather than 
the matrix h = h'^ on the fermion side, which is almost a matter of taste rather than a real 
simplification, one cannot object. However, especially in a second aspect seems to 
inform the opinion of the authors. It seems to be taken for granted that diagonalization 
on the fermion side can only be achieved after an orthonormal basis is constructed for 
the collective subspace. In holding this view they may have been influenced by the 
early papers in which the 5*0(5) model is introduced and analysed without the aid of a 
boson mapping. Consider for instance the following remark from a 1972 paper [21] by 
Chattopadhyay, Krejs and Klein: 

Although not mutually orthogonal, these are complete and form a conve- 
nient basis for the construction of matrix elements of the generators, and by 
extension of the Hamiltonian itself, which can be programmed in the form of 
recurrence relations. Just as for the problem of small vibrations, the diagonal- 
ization of the Hamiltonian is then carried out in two steps, the first involving 
the diagonalization of the metric tensor. 

The states refered to are simmilar to the ones in eq. 12. 441 that we worked with to compute 
h on the fermion side. The diagonalization of the metric tensor that is mentioned is 
nothing but the construction of an orthonormal basis, which precedes diagonalization in 
|24j . Because of the less trivial commutators of the building blocks of the states that span 
the collective subspace, constructing an orthonormal basis by means of a procedure like 
the diagonalization of the metric, is a complicated task. It involves the calculation of the 
inner product between fermion states of the form eq. 12.441 If this were unavoidable, then 
indeed, working with the mapped system would have been preferable by far. However, we 
have seen that we may diagonalize the original fermion Hamiltonian from the knowledge of 
its action on a non-orthogonal, possibly linearly dependent set of states. The calculations 
involved in this procedure mirror those in the mapped system to the extent that we end 
up with essentially the same matrix in both systems. The complications arising from 
the non-orthogonality of states in the original system is captured in the mapped system 
by the non-unitarity of the Dyson mapping. The complications associated with linear 
dependencies between states on the fermion side, are equivalent to those associated with 
the identification of ghost states in the mapped system. 

Our task is now to diagonalize the matrix h and find the ghost states if any, in its 
spectrum. For general values of the dimensionless parameter 9, this is done numerically. 
Later we will return to the specific values ^ = ±1 where (71 = and ^ = where (72 = 
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and, in trying to understand the behaviour of the ghost states, derive some results without 
resorting to numerics. 

In Section 11.91 we looked at two ways of identifying ghost states. The first procedure 
identifies the physical subspace by looking for a subspace of boson-fermion space that is 
left invariant by all physical operators. This procedure is implemented as follows for the 
system currently under consideration: First we diagonalize the matrix h for given values 
of Qi and §2- This matrix, and therefore the Hamiltonian Hbf as well, happen to have 
a non-degenerate spectrum, so that left and right eigenstates can uniquely be labeled 
by their eigenvalues. We adopt notation that denotes the right-eigenvector of Hbf with 
eigenvalue A by \\r). The left-eigenvector of Hbf with eigenvalue A is denoted |Al). We 
normalize these states such that (Al|A'j:j) = (5a, a'- In terms of the usual basis elements 
|?T^), the states \\r) and the functionals (A^l can be expanded as 

min(np,n„) 

|Ar) = E CA,„o|no), (2.48) 

no=0 
min(np,n„) 

(Ail = E KMno,A. (2.49) 

710=0 

If we adopt the usual convention that the first index labels rows and the second index 
columns, then the matrix c with entries CA,no contains on the row labeled by A the compo- 
nents of the left-eigenvector associated with eigenvalue A of the matrix h. The matrix d 
with entries dn^^x is the inverse of c. The column of d labeled by A contains the components 
of the right-eigenvector with eigenvalue A of h. 

For the algorithm we are implementing, we need another physical operator Q such that 
the only subspace left invariant by both Q and Hbf is the physical subspace. An operator 
that should fit the bill is Hbf with a different choice of the parameters gi and We 
take Q to be this operator. We define q as the matrix with entries qmo,no = (jk) \ Q |^)- 
Now we let q be the matrix cqd, so that the entry qx^y is given by 

^X,moQmo,no^no,\' ■ (2.50) 

mo, no 

It is not difficult to see that the following holds: if we expand Q\Xr) in terms of the 
right-eigenvectors of Hbf, we find 

Q|Aii)=E9A,v|A^), (2.51) 

A' 

SO that Q \ Xji) contains components along all the eigenbasis vectors |A'^) with A' such that 

gA,A' is non-zero. It is then by examining the entries in the matrix q in the way prescribed 

in Section imi that the eigenstates of Hbf that span the physical subspace are identified. 

^In general, a set of operators Qa may be needed. However, numerical results show that we are 
fortunate here to get away with only a single operator Q. If more operators were required, we could have 
chosen the mapped Hamiltonian, with several different choices of the parameters. 
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Continuing to use tlie same notation that we used above, the second procedure for 
identifying ghost states can also be discussed. This procedure involves the operator TZ 
whose null-space is the orthogonal complement of the physical subspace. The TZ operator 
as defined in Section 11.91 maps boson-fermion space basis states onto the states in the 
collective subspace of the fermion system according to the simple rule^ 

7^|^) = |r^o). (2.52) 

Ghost states in the spectrum of Hbf are now identified by applying the operator IZ to 
the left-eigenstates to find 

min(np,nn) 

n\\L)= E <o,aI^o). (2.53) 

rao=0 

If and only if this last expression is equal to zero, is the eigenvector associated with 
eigenvalue A a ghost state. Note here that this last result also follows from the analysis 
we did on the fermion side. Compare eq. 12. 531 with eq. I2.32l and the subsequent discussion, 
remembering that the matrix h is the hermitian conjugate of h. 

A problem with this method of identifying ghost states is that we need an economical 
method for determining whether the vector on the right of eq. 12.531 is zero. The most 
economical way that I can see, is to take the the inner product of the vector with itself. 
For this to work, we must be able to calculate inner products of the form (mo| rio), which 
we can do, if we can calculate the value of a general expression 

(oi a^,. {E_,_,r (E-ior {E-ur^^ {Enr (Eior (Ei-iT-^ ai^, |o) = c':;™x\- 

(2.54) 

This expression can be evaluated after commuting all pair annihilation operators past all 
pair creation operators, using the commutator identities of the iS'0(5) algebra. However, 
the number of commutations that need to be made, and hence the number of terms 
generated, is staggering. A better strategy is to commute just one of the pair annihilation 
operators past the pair creation operators, thereby expressing the inner-product between 
two (say) particle states in terms of inner-products between N — 2 particle states. This 
strategy leads to recursion relations for the inner-products. These recursion relations are 
fairly complicated to derive for the S0{5) algebra already, and will no doubt be even 
harder to derive for more elaborate algebras. Since I do not intend to develop further the 
method of 7^-projection for finding ghost states in a general setting, I do not write down 
these recursion relations here. It is far more practical to use the preceding method for 
finding ghost states in general. The 7^-projection method is of limited value in numerical 
^Recall that the operator TZ simply replaces each with and at with a+ in a state like 



Section 2.4: Diagonalization 



68 



work. It is liowever sometimes useful in an analytical as opposed to a numerical setting, 
as I hope to show in a little while. 

Everything has now been done except for actually diagonalizing our mapped Hamil- 
tonian, finding its ghost states and interpreting the results. Therefore we do this now. 
In Figure 12.11 we display the results of diagonalization and finding the ghost states. We 
chose a system with 9 protons and 20 neutrons. When ghost states have real eigenvalues, 
the eigenvalues are plotted as dashed lines. When ghost states have complex eigenval- 
ues, these eigenvalues are no longer displayed. When we analyse the results we display 
here, the focus is on understanding the 'ins and outs' of the mapping, not on interpreting 
possible physical phenomena encoded in the results. 



ID 
□1 

in 





Figure 2.1: Spectra for 9 protons and 20 neutrons, for various values of Vt. Real eigen- 
values associated with ghost states are included as dashed lines. 



Let us examine the spectra displayed in Figure ITTl At ^ = — 1 the coupling constant 
Qi is switched off and the operator Hbf is semi-negative definite. As 9 increases, the 



^-E'lo-E'-io 



BF 



term is switched on, with a positive coupling constant, and the eigenvalues 
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generally become less negative. At ^ = 1 the operator Hbf is just the negative of what 
it was at ^ = — 1. At the points 6 = ±1 and ^ = the operator Hbf can be diagonalized 
exactly, using algebraic techniques and the fact that these points in parameter space are 
SU{2) symmetry limits. Alternatively, for the ^ = limit, one can diagonalize Hbf by 
inspecting the matrix h. This matrix happens to be lower-triangular at the ^ = point 
in parameter space, so that its eigenvalues are the diagonal entries, and eigenvectors are 
particularly easy to find. Without going into further details, at 6* = ±1 the spectrum is 

Ekie = ±1) = ± / , {k K -np) + k {k + 1)} , (2.55) 

Up [Un + 1) 

where k = 0, 1, . . . ,np. In this equation, and from here on further, we assume, without 
any real loss of generality, that there are fewer protons than neutrons in the system. Note 
that in the 6* = ±1 limits the mapped Hamiltonian Hbf is hermitian and the eigenvalues 
do not depend on Q. 

In the ^ = limit the spectrum is 

k / 2Q — rir, — n.„ + A; + 1 \ 

where k = 0,1,..., n^, provided there are no ghost states present. This spectrum is 
almost equidistant for Q ^ rip + ?t,„. These results, which can be derived independently 
of all the arguments that lead to the numerically obtained spectra of Figure 12.11 fit the 
numerically obtained spectra perfectly, thus confirming the correctness of our treatment. 

The next thing to consider is the identification of ghost states. We would like to have 
some confirmation that our procedure for finding ghost states works. Also, we would 
like to understand how ghost states behave in dependence of both Q and 6. Looking at 
decreasing Q in Figure ITTl we see the first ghost state appearing at = 14. The eigenvalue 
associated with this single ghost state is real for all 9. Each time we further decrease Q 
by one, two more ghost states appear. Close to 6* = ±1 these ghost states all have real 
eigenvalues. As we get closer to ^ = 0, the eigenvalues associated with all but one of the 
ghost states become complex. The eigenvalues become complex in pairs. 

An elaborate analysis is possible, in which we derive all these results about the ghost 
states from scratch. However, in the interest of not boring readers more than is necessary, 
and of saving myself some writing, I take a less strenuous approach, that still paints the 
above picture emerging from the numerics, albeit in broader strokes. 

Firstly, let us consider the behaviour of the ghost states as we vary 9. The parameter 
9 determines the form of the operator Hbf that is diagonalized, but has no bearing on 
the mapping itself. The same mapping is used for all 9, and hence the number of ghost 
states does not change if we change 9. Furthermore, the fact that the operator Hbf is 
hermitian for ^ = ±1 suggests that there are regions close to 6* = ihl where the spectrum 
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of Hbf, ghost states included, is real. When there is only one ghost state present, its 
eigenvalue must be real for all 6. This is because the remaining eigenvalues are physical 
and hence real, and furthermore the operator Hbf has a real trace. 

Now we consider the behaviour of the ghost states as we decrease Q. It is not hard to 
understand why decreasing Q might increase the number of ghost states. It has everything 
to do with the exclusion principle and the fact that Q is proportional to the number of one- 
particle states in the fermion space from which we map. If we decrease Q while keeping 
the number of particles fixed, we decrease also the number of ways to arrange the original 
fermions among the one particle states they may occupy. No similar reduction happens 
in the boson space into which we map, and thus more of the linearly independent boson 
states must lie outside the physical subspace, whose dimension decreases as we decrease 
Q. A little more concretely we may recall that ghost states start entering on the boson 
side when the set of states v = {|no)}"^^Q that spans the collective fermion subspace 
becomes linearly dependent. Now, the state |no) is defined as 

with rii, n„i and a given by eq. 12.261 in terms of no, n„ and Up. For our present purposes 
it is necessary to know that of the three numbers rii, uq and ri„i, it is ?7,_i that may 
take on the largest value (assuming still that > Up). This value is n_i = ^ (for even 
numbers of neutrons). This happens for Uq = 0. Furthermore, thanks to the exclusion 
principle we have that 

'A^^Y = if n>Q, 



A^j =0 if n>2n. (2.58) 

Thus, we conclude that, when Q is smaller than or equal to half the number of neutrons, 
the state |no = 0) is the zero vector. This means that when Q hits the value ^ the set v is 
certainly linearly dependent and there must be ghost states. For smaller Q even more of 
the vectors |no) are zero and more linear dependencies should occur. However, the value 
r2 = ^ is an under-estimation for the threshold value of Q where ghost states enter the 
spectrum. This means that the set v becomes linearly dependent before such time as the 
state |rio = 0) is the zero vector. To get the correct value of Q below which ghost states 
enter, we dust off the TZ operator. From our numerical results we see that in the 9=1 
limit, ghost states start entering the spectrum from above. If n ghost states are present, 
they are associated with the n largest eigenvalues. (This is one of the results we would 
have derived from scratch if we did the more complete analysis.) So when there is just 
one ghost state present, it is the eigenvector of Hbf{0 = 1) with eigenvalue Ek=no = 1- It 
is not too hard to check that this eigenvector is given by 

np + n„-l -I , 

\9 = l,E=l) = iEo,)l^^{B~') ' (2.59) 
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I) contains one ideal neutron and {Eoi)^p is the image under the mapping 



where the state 

of the operator Eqi. To verify that |^ = 1, = 1) is indeed the eigenstate we are after, one 
uses the fact that (£"01)5^ and {Eo_i)^p are SU{2) ladder operators, and that (-Eo-i)^^ 
annihilates the state (B^^) ^ 



We do not distinguish between left- and right- 
eigenvectors because Hbf is hermitian when 9 = 1. If and only if the state \9 = 1, E = 1) 
is not a physical state, then, when we act on it with the operator TZ, we must get the zero 
vector. Acting with TZ we find 



n\e = 1,E = 1) = (E, 



01 j 



A' 



(2.60) 



where 



contains one real neutron. Looking at the exponent of the A ^ operator, we see 



that ghost states start to appear for values of f2 smaller than or equal to . For a 

system of 9 protons and 20 neutrons, this means that we expect ghost states for Q < 14, 
a prediction that is confirmed by the numerical results. 



2.5 Conclusion 

On this note, we conclude our analysis of the S0{5) toy model, that, in outhne, included 
the following. A mapping of a collective fermion representation of the 5*0(5) algebra 
together with single fermion creation and annihilation operators was derived. The map- 
ping is valid in the whole collective subspace of the fermion system, as opposed to the 
even part of the collective subspace only. This means that a single unpaired fermion may 
be present in the mapped system. A specific Hamiltonian, typical of algebraic nuclear 
models, was introduced. The diagonalization of the original and mapped Hamiltonians 
was considered. Ghost states were identified by finding the subspace of boson-fermion 
space that is left invariant by all physical operators. Numerical results were compared to 
some analytical results for corroboration and also to understand what goes on 'under the 
hood' of the mapping. A result that has emerged from our analysis is that, when we work 
in the collective subspace, there is a closer similarity between the original and mapped 
systems than was perhaps previously realised. 

What is still lacking is an example of how a generalized Dyson mapping may be a 
useful technique to understand the behaviour of certain many-body systems. We try to 
address this in the next chapter. 



Chapter 3 



An Application: The Richardson 
Model for Superconductivity 

3.1 Introduction and Motivation 

In this chapter we apply the techniques of Chapter One to a system with dynamics given 
by the Richardson Hamiltonian. This Hamiltonian describes fermions interacting in such 
a way that Cooper pairs are formed. In the sixties there appeared a set of papers in 
which R.W. Richardson used the model to understand pairing correlations in nuclear 
wave-functions. See for instance [H]. (In the bibliography of [2E1 a complete list can be 
found of the papers in which Richardson analyses the model that bears his name. I do 
not reproduce the list in the bibliography of this thesis, for fear of creating the impression 
that I have personally read all Richardson's papers.) What is notable about Richardson's 
analysis, is that it involves the exact solution of the model. In the context of super- 
conductivity in metals, the Richardson Hamiltonian is also thought relevant. Although 
the famous BCS paper |2ZI deduces the phenomenon of superconductivity from a more 
general Hamiltonian, for some calculations the problem is simplified by rather using the 
Richardson Hamiltonian. (In the notation of |2Zj, the interaction matrix elements Vk,k' 
are set equal to their average V, for |k|, |k'| < hu, and zero elsewhere. This simplifying 
assumption is made in several important calculations in the BCS paper.) Lately, authors 
have argued that there are fundamental reasons for considering the Richardson Hamilto- 
nian to describe very accurately the dynamics of electrons in mesoscopic metallic grains 
|28j . A thorough review of experimental and theoretical developments in this field can 
be found in [21]. Furthermore, the Richardson Hamiltonian has a simple but non-trivial 
algebraic structure: The superalgebra, the generators of which are the building blocks of 
the Hamiltonian, has an even sector that is nothing but the SU{2) algebra. 

On the one hand, applying a generalized Dyson mapping to a system described by 
the Richardson Hamiltonian, shows that the mapping can survive in the 'real world', 
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particularly when the mapped Hamiltonian cannot be expressed as a function of only 
the operators in the even sector of the superalgebra. On the other hand, aspects of the 
generalized Dyson mapping, such as the validity of the similarity transform of Section ITTTl 
for which our previous 5*0(5) example was too intricate, can be dealt with. 

There is another reason to do a mapping of the Richardson model: The Richardson 
model is exactly solvable, in the sense that expressions for the eigenstates and eigenvalues 
of the system can be derived in terms of a set of parameters that are solutions to a system 
of equations known as Richardson's equations. To find these solutions one has to resort to 
highly nontrivial numerics. This can only be done if the system does not contain too many 
particles. Also, in the limit of essentially infinitely many particles, the model simplifies 
to a form that is tractable [20]. If neither of these circumstances apply, or if one wants 
analytical rather than numerical results, there are two ways to proceed. The first more 
purely mathematical option is to focus on Richardson's equations and develop analytical 
techniques to solve these equations to a given degree of accuracy, with the minimum nu- 
merical work required. This was done with success by Yuzbashyan and co-workers in [21] . 
The other option is to forget about the Richardson equations and simply apply familiar 
quantum many-body approximation techniques to the Richardson Hamiltonian. The BCS 
variational treatment of superconducting systems in the thermodynamic limit (where the 
number of conduction electrons tend to infinity) is an example of such an approach [2Zj . 
The approximation method we will consider is conventional perturbation theory in both 
its time-independent and time-dependent forms. This approximation method can poten- 
tially yield information about systems with more than a few but less than infinitely many 
particles. The quantity with respect to which we expand is the inverse of a coupling con- 
stant, from there the term 'Strong Coupling Expansion' in the title of [211 • This quantity, 
which we will denote ^ is not very small in a superconducting metal, but the eigenstates 
and energies of the system are analytical functions of it, so that in principal we can get 
accurate answers by going to high enough orders in the expansion. Obtaining high orders 
however involves huge effort. In practice we therefore settle for truncating expansions at 
low orders. At best (if the properties of the system at small ^ are not too vastly differ- 
ent from the properties at moderate ^), this will give us a qualitative understanding of 
realistic systems. This does not make the second option inferior to the option chosen by 
Yuzbashyan and co-workers per se. In fact their method just gives an order for order solu- 
tion in ^ to the Richardson equations. Hence their method gives results equivalent to the 
usual time independent perturbation theory for the eigenstates and eigenenergies of the 
system. It is just the way in which they obtain these results that is specifically designed 
around the Richardson equations. Their method also runs into the difficulty that high 
order corrections are required for realistic systems. If we settle on the conventional way 
of doing a perturbation expansion, rather than the method of Yuzbashyan, we encounter 
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another problem though. Because of electron-electron correlations in the eigenstates of 
the zero'th order problem, the calculation of matrix-elements of the perturbation with 
respect to these states, which is the stock in trade of any perturbation expansion, is a 
formidable task. Here the boson-fermion mapping that was developed in Chapter One 
might come in handy. It turns out that this mapping transforms away the bothersome 
correlations, so that for the zero'th order mapped system, the eigenstates look like free 
particle states. Matrix elements of the perturbation are then more simply calculated in 
the mapped system. 

We will do both a time-independent expansion to obtain approximate eigenstates and 
eigenvalues and a time-dependent expansion to obtain approximate transition matrix 
elements. We compare the time-independent expansion for the mapped system with the 
results published in j^Hl- We find exact agreement. The method in [22] might however be 
more straight forward to implement when higher order corrections than we calculated, are 
required. The value of this component of our work then, is that it shows that it is feasible 
to do calculations using a mapped system obtained through the boson-fermion mapping, 
rather than that it allows for calculating properties of the system with unprecedented 
ease. More speculatively, the intuitive physical picture behind our method is enticing. 
Firstly, the Dyson mapping high-lights a few previously hidden features of the model in the 
strong pairing limit. Secondly, our results are obtained by means of a transformation that, 
crudely speaking, replaces Cooper pairs with bosons. Thinking about a superconducting 
system in this way might shed some light on the current hot topic of the Bose-Einstein 
condensation to BCS cross-over in cold atomic gases [HB]- Or it might not. This last 
speculation does not receive a systematic analysis in the present text. 

Since it seems that a time-independent perturbation expansion will not yield anything 
that is not within the scope of previous analyses, we will also do a time-dependent ex- 
pansion for the ground-state to ground-state transition amplitude, in the strong coupling 
limit. In doing this expansion I am not so much concerned with what useful information 
or measurable effect is exposed by the result. Rather, I illustrate that this quantity is 
easily calculated with the aid of the mapped system, whereas I suspect that it involves 
quite a complicated procedure to calculate it by any other means. 

Some of the matter of this chapter grew from work I presented in the form of a 
talk at the second international workshop on pseudo-hermitian Hamiltonians in Quantum 
Physics held in Prague (14-16 June 2004). A paper [33] based on the talk is to appear 
in the proceedings of the workshop. In that paper some of the simpler results from time- 
independent perturbation theory are derived. The paper constitutes a brief (and almost 
representative) sampling of the work contained in this chapter. 



Section 3.2: The Model 



75 



3.2 The Model 

We start by introducing the Richardson model, including the fermion Fock space on 
which it is defined. Therefore, let us consider a one-particle Hilbert space for which the 
orthonormal states 

k), = -J, -i + 1, • • • , j - (3.1) 
form a basis, where j is an odd multiple of a half. We also define a convenient integer 

^ = J + l- (3.2) 

The single discrete quantum number u combines an electron spin index and a Bloch wave 
number in such a way that the state {—lY'" |— i^) is the time-reversed state of the state 
lu). Note that j — u is an integer in the range 0, 1, 2, ... , 2Q — 1. This convention, which 
can be seen as a particular choice for the phase in the definition of the normalization 
constant for time-reversed states, is different from the one that is usually employed, but 
leaves measurable quantities unchanged. The reason for the unconventional convention is 
that it makes formulas for the structure constants of the emerging superalgebra simpler. 

With each state {u) is then associated anti-commuting fermion creation and annihi- 
lation operators and a^, that act in the fermion Fock space Tip, and which are each 
other's hermitian conjugates. In terms of these operators, and without the summation 
convention that we often employed in previous chapters, the Richardson Hamiltonian is 

H=f: e^aUu - ^ E (-l)'^-''-^a+a^^a_.a.. (3.3) 

The are the one-particle energies. We will assume time-reversal symmetry so that 
e-u = e^- The one-body term accounts as well as possible for everything that electrons do 
inside a normal metal while the interaction term contains only the mechanism for Cooper 
pairing. This interaction involves the scattering of time reversed pairs of conduction 
electrons among each other and is governed by a single coupling constant G with the 
dimension of energy. In order to apply a boson-fermion mapping to the system, we must 
identify a superalgebra such that the Hamiltonian can conveniently be expressed in terms 
of the superalgebra's generators. With this in mind, we define the collective fermion pair 
operator 

S+ = V {-ly-'ata^,, (3.4) 

that contains all time-reversed pairs with amplitudes differing only by a phase. Conju- 
gated to it is a pair annihilation operator 
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Using summation convention, we may also write 



s 



(3.6) 



where x,iu = ^(-'^V ""^t^-u and x^" = {x,iuT = Xt^v and a^" = (^a+). The operators 
S, and K = [S , S~^] are the generators of a representation of the SU{2) algebra. 
These, together with the single- fermion operators and form the generators of a 
single plus bi-fermion superalgebra of the type discussed in Chapter One. The important 
commutators are: 



S,S^ 



CI,, , 



U=~J 



The Richardson Hamiltonian can be expressed in terms of these generators as 



(3.7) 

(3.8) 
(3.9) 
(3.10) 



H= J2 - GS+S. 



U=-J 



(3.11) 



Since the operator S'^S is semi-positive definite, we take G to be positive so that the 
interaction —GS~^S lowers the energy. If we neglect the one-body term and solve the 
pairing-only problem H = —GS^S for 2N fermions, we find an unnormalized ground 
state! ^g+y^ rj.^-^ 

means that the operator creates Cooper pairs, while S annihi- 
lates them. 



3.3 Properties of the Eigenstates of the Richardson 
Hamiltonian 

Before we map the Richardson model onto a boson-fermion model, I briefly state a few 
facts about its eigenstates. The discussion is not intended to be exhaustive, nor is it 
essential for performing a boson-fermion mapping. (The boson-fermion mapping and 
subsequent analysis of the system do not require anything more than the fermion Hamil- 
tonian to be known.) The results I mention here should rather be seen as a reference 
point which we will refer back to when we interpret the behavior and properties of the 
mapped equivalent of the Richardson model. 

^Throughout this section of the text, unless otherwise indicated, eigenstates will not be normalized. 
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The first thing to consider is the so-called blocking effect, which effectively reduces the 
problem of finding the eigenstates of the Hamiltonian in the whole fermion Fock space 
Tip to the problem of finding the eigenstates in a certain subspace of Tip- 

To explain the blocking effect, the following terminology is introduced. We refer to 
the one-particle states |— z/) and \v) together as the one-particle level v. The fermion Fock 
space Tii? is a many-particle space in which each one-particle level u can be unoccupied, 
singly occupied or doubly occupied. Electrons in singly occupied levels do not take part 
in the dynamics in any other way than that, due to the exclusion principal, they prohibit 
scattering into the levels that they occupy. For instance, consider the two-electron state 
a^a^ |0) with 7^ Since the operator S can only remove particles from doubly occu- 
pied levels, the pairing interaction does not notice the two unpaired electrons, and we have 
S^Sa^a^ |0) = 0. Hence the Richardson Hamiltonian endows the state a^a^ |0) with the 
same dynamics as if the electrons were not interacting: Ha'^a^ |0) = (e^ + e^,) a^a^ |0). 

To make these ideas more precise and extend them to an arbitrary number of particles 
we introduce a specific kind of basis for M-particle fermion Fock space. For each element 
of the basis we are talking about, we specify two sets of indices. The first set B is called 
the set of blocked or singly occupied levels. Each set B contains no more than M of 
the indices v = — j, . . . ,j, with the condition imposed that if the index u is an element 
of B, then the index —u is not an element of B. Let b be the number of elements of 
B, then we insist that b has the same parity as M. The reason is that we want the 
indices in B to denote one-particle levels in Tip that contain only one electron. The 
second set D specifies the doubly occupied levels. It contains ^^^^ indices from the set 
U = . . . , j| \ {\u\ : u G B} of levels that are not singly occupied. These are called the 
unblocked levels. Each basis element is then uniquely specified as 



B,D) = 




(3.12) 



The collection of all possible B partitions the basis into disjoint subsets Pb that are 
defined as follows. For a given set B of singly occupied levels, let Pb be the set of all 
states \B, D) such that the set D of doubly occupied levels is consistent with the choice of 
singly occupied levels in B. The Richardson Hamiltonian leaves each subspace span {Pb) 
invariant. In fact, a little algebra reveals that 

H\B,D) = (y: e.) \B,D) + ( n < I f n io) , (3.13) 

where the operator H' is the Richardson Hamiltonian redefined on the levels that are not 
singly occupied, i.e. R' = J2^.,\^\<.u + G"S"+S", with S"+ = -^^^ J2ueui-y~''o,tO'-u 

and G' = G^^. This means that to diagonalize H on the whole fermion Fock space 
Hp, we only have to diagonalize the various H' (each of which looks just like H but 
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with fewer levels and a different coupling constant) for the electrons distributed among 
doubly occupied levels, and afterwards add the electrons in singly occupied levels and 
their energies. 

The problem considered in this text then, is diagonalizing H for an arbitrary set of 
Q one-particle levels and with an arbitrary coupling constant G, for an arbitrary even 
number of particles {2N), but with no singly occupied levels. The space in which H is 
diagonalized has dimension (^^^. If we solve this problem, the problem of diagonalizing 
H with singly occupied levels taken into account, is solved without further complication. 

The second thing I now discuss is a solution of sorts to the problem of diagonalizing the 
Richardson Hamiltonian when there are no singly-occupied levels. We consider a system 
of 2N electrons. Let E^, k = 1,2, . . . , N he any set of complex numbers that satisfy the 
so-called Richardson equations 

E k=l,2,...,N. (3.14) 



^ 1 Ek — 2eu ,_i /I, El — Ek 



There is a correspondence between solutions of eq. 13.141 and the eigenstates of H with 
no singly occupied levels. With every eigenstate \E) of H with no singly occupied levels, 
such that H\E) = E\E), is associated a solution Ei, . . . , E^- of eq. 13.141 (and vice versa, 
with every solution of eq. 13.141 is associated an eigenstate of H with no singly occupied 
levels), such that the energy of the state is given by 



N 

E=Y.Ek, (3.15) 



k=l 



and the eigenstate has the form 



l^) = n E^^P^IlO). (3.16) 
fc=i \ 1 - Ek 



To derive this is not hard but will not be done here, for two reasons. Firstly, a good 
derivation is available in the literature [IHl Appendix B]. Secondly, we will shortly en- 
counter an operator that is very similar to the Richardson Hamiltonian. Its eigenstates 
have properties very much like those of the Richardson Hamiltonian and the derivation of 
these properties is accomplished using the same strategy as with the Richardson Hamil- 
tonian. Since this new operator has not yet featured in the literature, I rather do the 
derivation for its eigenstates. 

To claim that the above result completely solves the problem of finding the eigenstates 
of the Richardson Hamiltonian with no singly occupied levels, is overly optimistic. The 
reason is that we still need to find all solutions of eq. 13.141 This is equivalent to finding 
all simultaneous roots (of which there are , one for each dimension of the space in 
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which H is diagonalized) , of polynomials of variables and of degree NQ. Since 2N 
is the number of electrons and Q the number of one-particle levels, these two numbers 
are both large in the truly many-body regime. This makes the solving of eq. 13. 141 a very 
daunting task and one, it hardly needs to be said, that can only be treated numerically. 

The last matter to be discussed is the behavior of the solutions to the Richardson 
equations as the coupling constant G in front of the pairing term becomes large when 
compared to the one-particle energies e^. Specifically, we mean 



^0, (3.17) 

so that, if we take the one-particle energies as our yard stick, the left-hand side of eq. 
13. 141 tends to zero. As this happens, some of its solutions have components that decrease 
without bound while there are others all of whose components stay finite [22] • To be 
precise, of the (^^^ solutions {Ei, . . . , En) to eq. 13. 141 there are only (^^^ x ^"^^^^^ whose 
components Ek, k = 1, . . . , N all remain finite as G — > oo. If we set ^ = on the 
left-hand side of eq. 13.141 only these solutions survive, the others having vanished off into 
negative infinity. Since the energy of an eigenstate is given by the sum of the compo- 
nents of its corresponding solution to the Richardson equations, this means that even if 
the coupling constant G is infinite, there are still (J^^ x ^^^^^^ eigenstates with finite 
energy. This concludes the discussion on the properties of the eigenstates of the original 
Richardson Hamiltonian. 



3.4 The Mapping 

We now construct the boson- fermion mapping from fermion Fock space Tip, or a sub- 
space thereof, into the boson-fermion Fock space Ti-BP, that results from the superalgebra 
introduced above. First we must decide which subspace of Tip (which may turn out to 
be the whole Tip) we want to map into H-bf- In Section HTHl we saw that we could avoid 
growing a few gray hairs when we restricted ourselves to a so-called collective subspace 
of Hp- In the present case the (even sector of the) collective subspace is spanned by the 
single state 

(5+)"|0). (3.18) 

The Hamiltonian H only leaves the collective subspace invariant when the one-particle 
energies e^, are all set equal to each other. 

It is evident that we cannot map the collective subspace alone. Therefore we set foUical 
concerns aside and map a larger subspace of Hp into Hep- The program of Section ITTTI 
must be followed, where states in Tip that contain more than one fermion that is not in a 
collective pair are mapped onto states in TCbp using the compounded operator o T. 
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Here T is the mapping of Tip into T-Cbf that results from the Usui operator constructed 
in Section 11.61 and X is the similarity transformation defined in Section 11.71 (Without 
the similarity transformation, the association between Cooper pairs and bosons is not as 
close as we would like.) 

However, the term 'defined' is used loosely, since it was mentioned in Section im that 
the defining equation for X, eq. 11.511 is problematic. Our first task then is to determine 
if and when X is well-defined for the superalgebra we are currently working with. 

We set up our notation as follows. With the boson-fermion Fock space T~iBp are 
associated boson creation and annihilation operators and B, as well as fermion creation 
and annihilation operators aj, and a^. We refer to these fermion operators as ideal fermion 
operators. They commute with the boson operators. We also define ideal collective 
fermion pair operators 

5^ = E {-iy-'c.W-u, (3.19) 

and S = {'S'^y which are the counterparts in T-Cbf of the collective fermion operators 
and S that are defined in Tip. Mb = B'^ B is the boson number operator and Mp = 

Yl,u o^luy is the ideal fermion number operator. 

The problematic defining equation (eq. ll.Slj) for the similarity transformation X now 

reads 

X 1^, np, UB, A) = E {(Cf - A)-' S^b} ' |^/;, np, ub. A) . (3.20) 

fc=0 

In this equation Cp is hermitian and given by Cp = S'^S. The state 1-0, np, hb, A) G Hbf 
is an eigenstate of Cp with eigenvalue A. Since Cp commutes with the ideal fermion 
number operator and the boson number operator, we choose 1-0, n^?, n?,, A) to have a well 
defined number of ideal fermions, namely np, and a well-defined number of bosons namely 
ub- ijj refers to any other simultaneous eigenvalue (s) that are required to label the state 
uniquely. Since Cp has an eigenvalue A, Cp — X has an eigenvalue zero, and hence it is not 
invertible on the whole Ti-BP- It is however invertible on the orthogonal complement in 
Hbf of the A-eigenspace of Cp. This space is denoted by Aj^. The operator {Cp — X)~'^ 
refers to this inverse of Cp — X in A±. Hence it is an operator that is only defined on the 
space Aj^. 

The implication for the attempted definition in eq. 13.201 can now be tackled. We have 
to establish whether every term 

\n) = {iCp- Xy' S^b}' np, ns. A) , (3.21) 

makes sense. It turns out that particle number considerations are important in this 
argument. To avoid muddling further down the line, we recall a few facts and define some 
notation. Since the Richardson Hamiltonian preserves the number of fermions, it makes 
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sense to choose a subspace Vp ^ Tip that is left invariant by the fermion number operator 
J2i=-j o.t'^yi the domain of the mapping o T. One can then find a basis for Vp in 
which each basis state has a well-defined number of fermions. Let |0, M) be such a state, 
containing M fermions.^ A property of the operator T, as defined in Section 11.61 is that 
T 10, M) is an eigenstate of + A/f with eigenvalue M. We therefore choose a basis 
for 1-Lbfi the elements of which have both a well-defined number of ideal fermions np and 
a well-defined number of bosons ub- If we then demand that states in Vp have no more 
than Mmax fermions, we are guaranteed that the range of T, namely T(yp), is contained 
in the subspace of Ti^F spanned by states for which twice the number of bosons plus the 
number of ideal fermions is no more than Mmax'- 

np + 2nB < Mmax. (3.22) 



Now we investigate when eq. I3.21l is acceptable. Because of the commutator S , S'^ 

JVf 

n ' 



the operator identity 



CpS^B = S^B (Cp+l- ^\ , (3.23) 



holds. From this it follows that 



CpS^B np, Ub, ^) = + ^ " ) ^ |^/;, np, ub, A) . (3.24) 

B\il),np,nB,\) is an eigenvector oi Cp with eigenvalue A + 1 — In parenthesis we 
note that for \Ti) to be well-defined, nj?, ra^. A) must lie in A_l. From eq. 13.241 

follows that, because Cp is hermitian, S'^ B 1-0, np,nB, A) lies in A_|_ if and only if np ^ Q. 
This condition is enforced by demanding that the states in the subspace Vp ofTip contain 
less than Q fermions. That is to say, we consider a conduction band at less than half filling. 
(States in Tip can in principal contain up to 2^ fermions.) When the conduction band is 
less than half-filled, it follows that np and ri^ < . Since we get the state |Tfc), if 
it is well-defined, from the state np,nBi A) by removing k bosons (and adding 2k ideal 
fermions), |Tfc) is the zero vector in TisF for k > . 

Continuing in the same vane as in eq. 13.241 we can use the identity ()3.23|1 to compute 



CpiS^'B^j \ilj,np,nB,X) 



k 

X + -{n-np-k + l) 



S^BY\ilj,np,nB,X), (3.25) 



or in other words [<S^ B) \'tp,np,nB, X) is an eigenvector of Cp (when it is not the zero 



vector), with eigenvalue X+^{Q — np — k + l) which, because of the restrictions on np and 

^In deriving the mapping, we work as generally as possible and don't make the assumption that 
we work with an even number of particles. This will allow us to derive images for the single-fermion 
operators. 
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k that we have derived above, is different from A. This means that for k < the state 

{<S'' {i/j, Up, riB, A) lies in A_l. For k > , the state np, ub, A) is the zero 

vector, which also lies in A_l. It follows that we may always act on [<S'^ Bj |?/^, nj?, n^. A) 
with {Cp — A)~^, the result being 

{Cp - X)-' [S^bY np, UB, A) = -rip-k + 1) ('^^^)' l"^' ''''' ' ^^'^^^ 

where {^^B^ \ilj,np,nB, X) becomes the zero vector before the eigenvalue becomes sin- 
gular. The result in eq. I3.2fil is used to compute 

\n) = {Cp-X)-'S^B...{Cp-Xr'S^B\ij,np,nB,X) 
V ' 

k times 

- ^^'""^'^^--n'iS^BYlilj^np^nB^X). (3.27) 



{n-npy.k\ 

In the second line we simply allowed the [Cp — X)^^ operators to act from right to left on 
whatever is in front of them. Since the 'whatever' is either an eigenvector of {Cp — A)~^ 
or the zero vector, the {Cp — A)~^ operators were replaced by numbers that were then 
collected in the curly brackets. This last calculation confirms that the similarity transform 
X is well-defined, provided the conduction band is less than half-filled. Notice that eq. 
13.271 is independent of the value of A and so that linear combinations of states for 
which these differ may be taken, with eq. 13.271 remaining valid. This is done implicitly 
by taking up these labels in the generic label ip. Thus, for a general state lip, np) G H-bf 
containing np ideal fermions, the action of X can be stated as 

^ 1^' ^f) = t ^o'^'^Tw? ^' i^'B)' 1^, np) . (3.28) 



fc=0 



{n-np)\k\ 



Note how fortunate we are to have this last result: The source of our good luck is essen- 
tially the fact that the commutator Cp ,S^B is a function of the ideal fermion number 
operator only (eq. 13.23^ . For more complicated algebras, similar result do not hold. 

Having confirmed that the similarity transformation X that was introduced in Section 
11.71 does indeed exist for the superalgebra that we deal with here, we can write down 
the images of the generators of the superalgebra under the mapping. Using the general 
formulas of Section 11.71 adapted to the current instance, we get for the images of 5" and 



S < — B, (3.29) 

+ . 



^ i?t(l- ^^+^^ ). (3.30) 
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Because of the explicit expression for X, that is eq. I3.27| we can derive explicit expressions 
for the images of the single fermion operators and (which is more than we were able 
to do in Section First to be considered is the single fermion annihilation operator 

tty. Its image under o T is an operator {au)x bf '■ '^bf — * 'Hbf defined by^ 



X(a 



i^JX.BF 



(3.31) 



In order for us to use the explicit expression (eq. I3.27j) for X, we have to let the operators 
^i.o,y)xBF ttj/X act on a state \ip,nF) G TIbf that contains np ideal fermions. I 
therefore start to manipulate the state a^X \ip,np): 



(fi-n^)!/! 
^ {Q-nF-l)\ 

h i^-riFV-n 



\i{j,nF) 



(3.32) 



We now treat these two terms separately. Term (1) is simplified if we use the commutator 



a 



to compute [a^ , {S^y\ = l{S^y~^aU. Then 



oo 



°° {Q-nF-k-iy. 



E 

fc=0 



iSl-npyM 



\tp,nF) 



k f-iv-^ 



- (fi-(np-l))!fc! 



|V',?^F) 



(3.33) 



As the vector in curly brackets contains np + ^ ideal fermions, the operator that it is 
preceded by represents the action of the similarity transform X. The final form of term 
(1) is then 



h i^-riFyn 



OLv 



S^B 



\^J,nF) = X 



{-ly-^Bal 



(3.34) 



■^Strictly speaking, because we work with a conduction band at less than half filling, the domain of 
this operator is not the whole Hbf, but this technicality is innocuous and may be forgotten. 
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Term (2) is split up as follows 



5tS 



{Vt-nF){k+l) 



\iJj,nF) \ . 



(3.35) 



We can write |'0, np) as the first term of Xa^ 1-0, nj?), i.e. 



E 



[ - - k)S^Ba^ 



\{Vt-nF){Vt-nF + l){k + l) 
By substituting this expression into eq. 13.351 one finds for term (2) 



|^,n^). (3.36) 



(fi-n^)!/! 



^ ^n-{nF + l)-k)\ ^,f^,^^k J 1 5tga , 



rip 



1-^1 



np- 



(3.37) 



We now combine term (1) and term (2), replacing the number uf with the operator JVf 
to conclude that the image of is given by 



a,y < — a 



1 



1 



(3.38) 



The easiest way to find the image [o.'^^xbf ^'^is point is by exploiting the 

commutator 

(3.39) 



X,BF 



V^(-)^'+'^ 



,(s-) 



X,BF 



which is one of the basic commutators between generators of any representation of the 
superalgebra. Since the right-hand side of the equation involves the images of a„j, and 
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S~^, both of which we have aheady computed, we quite simply calculate 



(3.40) 

The images of the singel-fermion operators of eq. 13.381 and eq. 13.401 can be compared 
to those derived in that were derived by essentially the same method. We see that 
those expressions differ from eq. 13.381 and eq. 13.401 by some fi-dependent factors in 
front of the various terms in the images. The reason for the difference is that in a 
different normalization is used for the opperator that features in the definition of the 
Usui operator. This means that the boson-fermion representation we work with here is 
different from that considered in Both representations are equally valid though. 

All the building blocks required to write down the image {H)-^ of the Richardson 
Hamiltonian 



H= J2 ^yCila^ - GS+S. (3.41) 

u=-j 

are now available. The image of the Richardson Hamiltonian is then written as 

{H)^,,p = K + Ho, (3.42) 
where the operator Hq, the image of the pairing term, is given by 

Ho = -^ATb (n + l-AfF-AfB). (3.43) 

The operator K, the image of the one-body term, can be written as the sum K = + 
Kq + of an operator that annihilates a boson and creates two ideal fermions, an 
operator Kq that leaves both the number of bosons and the number of fermions unchanged, 
and an operator that creates a boson and annihilates two ideal fermions. After a little 
algebra we arrive at the following expressions for these operators: 

K_ = B^i, (3.44) 



njn-Mr- Mb) 

{Q - ATf - 1) {ft - ATf) 



njn-Afp- Mb) 



^ ^ {Q-ArF-l)i^-W i^~AfF + l) 



Section 3.4: The Mapping 



86 



In these expressions we exploited the fact that = e^^ to affect some simplification. The 
operator ^"l" is defined as 

^^ = ^t^ 2e./?t. (3.47) 

Here the operator Pi = {—y^'^a\Q^_y creates time-reversed pairs of ideal fermions. In 
deriving the expressions in eq. I3.46[ we used the symmetry property jS'Ly = (3l. The 
average one-particle energy is denoted e = Z^^^i ff- 

The image of the pairing term is simple. It is only a function of the boson and ideal 
fermion number operators. The diagonalization of this term alone is accomplished by 
simply using a basis for the boson- fermion space Ti-BF the elements of which have well- 
defined numbers of both bosons and ideal fermions. If the original system contained 2N 
fermions, the physical subspace is contained in the subspace of T-Cbf spanned by all states 
with N — s bosons and 2s ideal fermions where s may take on values 0, 1, 2, . . . , A^. Such 
states may be denoted \il^,s), where "0 refers to all other quantum numbers needed to 
specify the state uniquely. In this space Hq has a spectrum 

^(0) = -9.(^N-s)in + l-N-s), (3.48) 

with s = 0, 1, 2, . . . , A^. As mentioned previously, for the mapping to be valid, more than 
half of the 2Q states of the original system must be unoccupied. This brings about the 
restriction that 2A^ < f2 and thus also 2s < il. With this restriction the energy E^'^^ 
is a monotonically increasing function of s. The number s, which in nuclear physics is 
called the seniority, labels the energy levels of the unperturbed system, with larger s 
referring to higher excited energy levels. All but the lowest (s = 0, i.e. bosons only) 
level is degenerate. The spacing between level s and level s -|- 1 is G (^1 — , which is 
almost equidistant for low-lying states. What is harder, is identifying the ghost states in 
this spectrum. It turns out that the original system does have eigenstates associated with 
each E^'^^ of eq. 13.481 The ghost states only affect the degeneracy of the levels. Ghosts are 
present in every energy level, except for the ground state level, which is non-degenerate. 

On the other hand, what was the trivial one-particle contribution in the original 
system, has now become the non-trivial contribution in the mapped system. Where, 
in the original system, we had only fermion number operators a'^a^, we now also have 
contributions that scatter bosons into time-reversed fermion pairs and vice versa. Any 
scattering event can be decomposed into a series of first order processes, where there are 
only three possible varieties of these first order processes. Firstly there is the process, 
where the system loses two ideal fermions and gains a boson. Then there is the Kq process 
where two ideal fermions pair scatter off each other. Lastly there is the process where 
a boson decays into two ideal fermions. This result is so important that I restate it in a 
slightly different way: There are only three first order scattering processes. The first is 
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represented by the operator K^. It takes states from the energy level Ef''^ to states in 
the level directly below i.e. the Ef}i level. (The ground state is annihilated by K_.) The 
second process is represented by Kq. This is an 'elastic' process in the sense that it leaves 
states in any level Ef^ in that same level. Lastly there is the process represented by Kj^ 
which takes states in the level E^^^ to states in the level directly above, i.e. the -E^+i level. 

This structure makes the mapped Hamiltonian well-suited to perturbation expansions 
where we treat the image of the pairing term as the unperturbed Hamiltonian and the 
image of the one-body term as a perturbation. In order to be able to truncate such 
expansions after a finite number of terms, we assume that the coupling constant G in front 
of the pairing term is large compared to the one-particle energies e^. Expansions of this 
type are therefore called strong coupling expansions. We will do both a time-independent 
perturbation expansion for the spectrum of the system, and a time-dependent expansion 
for transition amplitudes. 

3.5 Time-Independent Perturbation Theory 

First on the agenda is time-independent perturbation theory. We are dealing with a 
non-hermitian Hamiltonian. To complicate matters, the unperturbed Hamiltonian has a 
highly degenerate spectrum. Therefore I start the discussion by presenting the general 
formalism for finding perturbatively the spectrum of a non-Hermitian operator if the 
spectrum is degenerate when the perturbation is switched off. The argument is developed 
along the lines of the non-degenerate hermitian treated by Sakurai [SI]. 

We are given an operator H{X) = Hq + XHi. Neither Hq nor Hi has to be hermitian. 
We assume that H is fully diagonalizable in a finite region around A = 0. For A 7^ 0, we 
assume that ^^(A) has a non-degenerate spectrum, so that a single label a = 1,2, . . . ,fl 
may be used to specify eigenvalues and eigenvectors uniquely: 

H{X)\aR)^ = E^{X)\an)^, (3.49) 

with Ea{X) 7^ Ep{X) ii a ^ [3 and A 7^ 0. The subscript R is used to indicate that we 
are dealing with right-eigenvectors, later to be contrasted with left-eigenvectors, that do 
not coincide with the right eigenvectors when H{X) is not hermitian. We take it that 
all the Ea{X) and \c(r)^ are analytical functions of A in a finite region around A = 0. 
Furthermore, if two eigenvalues converge at A = 0, i.e. Ea{0) = Ej3{0), we assume that 
these eigenvalues have different first derivatives at A = 0. This condition is referred to 
by saying that all degeneracy is lifted in the first order. At A = the spectrum of H 
may be degenerate. Our notation takes care of this as follows: At A = we let Q' be 
the number of distinct eigenvalues of Hq. We then indicate these distinct eigenvalues 
of Hq as £"/, / = 1,2, ... , fl'. We partition the labels a = 1,2, ... ,Q into disjoint sets 
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Pi, labeled by integers / = 1,2, .. .Q', by defining Pi to be the set of all indices a such 
that Ea flows to Si as A goes to zero. We can then also define a function p from the 
index set 1,2, ... ,Q to the index set 1,2, ... ,Q' as follows: p{a) = I where / is the 
unique label such that a ^ Pi. We refer to the space spanned by all eigenvectors of Hq 
with eigenvalue Si as the £^/-eigenspace of Hq. It is assumed we know a basis for each 
Si eigenspace of Hq. The basis elements that span the Si eigenspace will be denoted 
\'iI^,Ir) where -0 is a discrete index that runs from 1 to some integer oji which is the 
dimension of the Si eigenspace. The set //j) : I = 1, . . . ,Q';'iIj = 1, . . . ,lji} forms a 
basis for the domain of H{X). It is in terms of these basis states that we want to express 
the eigenstates of H{X) approximately. Corresponding to this basis we can always find a 
left basis {1-0, Z^) : / = 1, . . . , -0 = 1, . . . , ui} such that the equations 

Il\ 4>, rrin) = Si^^S^^^, (3.50) 

hold for all l,m,ip and (j). 

We uniquely define a projection operator 11/ as the linear operator that maps any vector 
in the £^/-eigenspace of Hq onto itself, while mapping any vector that can be written as a 
linear combination of vectors from the other eigenspaces of Hq, onto zero. (This might 
not be an orthogonal projection operator though. The different £^;-eigenspaces are usually 
not orthogonal to each other as Hq is generally non-hermitian.) An operator 11; is defined 
as Hi = I — 11;. In terms of the states and 10, m/j), the identity operator I, the 

unperturbed operator Hq and the projection operators 11; and 11/ can be expressed as 
follows: 

I = EElV',^fl)(^,/L|, (3.51) 

Ho = E^'El^'^R)(^'^il' (3-52) 
1=1 ^=1 

= J2\^^^R)i^^^L\, (3.53) 
ip=i 

n« = E EI^'M(^,^lI, (3.54) 

Note that HlUm = ^Im^m- 

With the above preliminaries out of the way we now expand the eigenvalues Ea{X) 
and eigenvectors \aFi)x H{X) in terms of A: 

Ea{X) = + AA« + A^A(f ) + . . . , 

Mx = \a^S^)+X\Ji^)+X'\a^i^) + .... (3.55) 

This definition fixes the directions of the A-independent zero'th order eigenvectors ctR^^^ 
uniquely as the direction to which the a-eigenvector of H{X) converges as A goes to zero. 
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This direction is well-defined since we assumed that H{X) has a non-degenerate spectrum 
for A 7^ 0. The magnitudes of the eigenvectors are left arbitrary. 

The set of eigenvectors | ct^R^ | ^ also forms a basis for the domain of H{X). We use 

it to define a second basis (sometimes called the contra- variant basis) | ^l'')} ^ through 
the set of equations 



(4°' I 



Eg 13. 561 implies that [a 



„(o) 



6a,f3- (3.56) 

Using the states | tti^^^j ^ and functionals 

I (ql^ I ^ we define another set of projection operators that will come in handy during 
our calculations. Define a projection operator 



■p(a) 



(0) 



a 



(0) 



a 



(0) 



(3.57) 



that singles out the a^'') componet of the decomposition of any vector in the | cir^^^ ^ 
basis. Let tTq, be the complementary projection operator of tTq in the Sp(^a) subspace by 
setting 



E 

/36-Pp{a)\{a} 

Thus, the identity operator may be decomposed as 



n 



I = n, 



p{a) 



7Tn 



(3.58) 



(3.59) 



The next step is to choose normalization for the eigenstate |tt/j)^. The most convenient 



choice turns out to be (^a 



(0) 



Or 



1. With this choice it follows that 



a 



TTq: \ aR)x 

Now we rewrite the eigenvalue equation in the form 

(^Spi^a) - Hq^ \aR)x = (Ai^i - A„(A)) {ur)^ , 

where A,(A) = Er=i ^i^^U^ 

The first thing we do with this equation is multiply it by Ilp{a), noting that 



(3.60) 



(3.6i; 



p(a) 



He 



0, to find 



(3.62) 



which has to hold order for order in A. Note that, by definition, the zero'th order term 
in the A-expansion of \aR)^^, namely a^R^, is an element of the Sp(^a) eigenspace of Hq 

I. Keeping this in mind, we look at the terms in eq. 13.621 of 



and hence 11 



p{a) 



a 



a 



(0) 



lowest order in A to find 



^p(a)Hi 



a 



a 



(0) 



(3.63) 
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This eigenvalue equation enables us to find the zero'th order vectors a^'') as the eigen- 
vectors of Ilp(^a)Hi in the Sp^a) eigenspace. Since we assume all degeneracy to be lifted in 
the first order, no two eigenvalues A^^^ and A^"'^^ with a, /3 G Pp{a) are the same, and eq. 



13.631 is necessary and sufficient to determine the various a^j in the Sp(^a) eigenspace. 
In order to proceed, we define an operator 

Qk = Y.y^^i- (3-64) 
l^k ^' 

We return to eq. 13. (ill and multiply it with Qp[a), noting that Qp(^a) {fp{a) — Hoj = 6^(0). 
Thus we arrive at one of the most important formulas in the present discussion, namely 

%ia) |aR)A = Qp{a) {^Hi - Ac,(A)) Ior)^ . (3.65) 

This equation has, on the left-hand side, the part of the eigenvector {aRj^ that can be 
writen as a linear combination of vectors in Si eigenspaces with I ^ p{o:)- We want a 
similar equation for that part of \(yR)x that lies inside the £^p(Q,)-eigenspace. 
For this purpose, we firstly rewrite eq. 13.621 to read 



Aa{X) - Xn.p(^a)HlIlp(^a)) '^-p^a) l^njx " >^'^p{a)HlIlp(^a) l^njx = 0, (3.66) 

by recalling that 11^(0) + Ilp[a) = I and (llp(^a)J = ^p{a)- Now we substitute eq. 13.651 into 
the second term of eq. 13.661 to find 

(Aa(A) - Xn.p(^a)Hin.p(^a)) '^-p^a) I^RIx ^ >^^p{a)HiQp(a) (A-f^i - Aq(A)) |a/j)^ = 0. (3.67) 

Define an operator 



A A ^ ■ 



Note that Qa y^a'^ — ^p{a)HiIlp(a)) = '^a and that qaT^p{a) = Qa- Hence, by multiplying 
eq. 13.671 with and deviding by A, to find 

oo 

Mx = qaHiQp(^) {\H, - A,(A)) \aR)x - E ^L'"''^A'=g„ \aR)x ■ (3.69) 

k=l 

This is the counterpart of eq. 13.651 that we needed for the part of la^);^ that lies inside 
the £^p(a)-subspace. When we combine the results of eq. 13. 60^ eq. 13.651 and eq. 13. 69^ 
remembering that the identity operator may be decomposed as I = Ilp[a) + + ttq we 
find 

oo 

I««)a = \^r) + (1 + QpW (Ai^i - A„(A)) \an)x - E AL'+'U'g„ |«fl), . (3.70) 

k=l 
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Collecting terms of order + 1 we find, for > 0, 



a 



(JV+l) 
R 



[1 + QaHi) Qp(a) \ Hi 



a 



(TV) 



N+1 
M=l 



a 



R 



N+2 
M=2 



a 



(N+2-M) 
R 



(3.71) 



Note that Qp^c 



a 



(0) 



and Qa \a^li 



0. Thus, for > 1, the upper bounds on the 



two summations in the above expression may respectively be decreased from + 1 to 
and from A^ + 2 to A^+ 1, while for the case where A^ = 0, the two summations may be left 
out entirely. To complete the expansions, we need an expression for This we find 



by multiplying eq. 13.611 with the functional 



a 



recalling our normalization convention 



(o-L^ (^r) ^ = 1; cind collecting terms of order M in A to find 



a 



(0) 



Hi 



a 



(M-l) 
R 



(3.72) 



Our final results are then: 
,,(^+1) 



a 



R 



1 + QaHi) Qp(a) \ Hi 



a 



R 



N 



^ A(M) 



M=l 



a 



(^r+l-M)^ 
R 



N+1 

E Ai^^^go 

M=2 



a 



{N+2-M) 
R 



and 



a 



(0) 



Hi 



a 



(N) 
R 



(3.73) 



(3.74) 



It is understood that, in the case where A^ = and the lower bounds of the summations 



in eq. 13.731 exceed the upper bounds, the summation contains no terms, i.e. 



QaHi) Qp(a)Hl 



a^^^y The formulas of eq. 13.731 and eq. 13.741 are all we need to 
calculate the corrections of order A^ + 1 to both the eigenvectors and eigenvalues of H, 
if the corrections of orders 0, 1, 2, . . . , A^ are known. Together with the initial conditions 
provided by the eigenvalue equation 



a 



^p(a)Hi 



a 



Ad) 



a 



(3.75) 



this allows us to write down in principal the expansion of any eigenvector and its eigenvalue 
to arbitrary order in A. Concerning eq. 13. 741 for the eigenvalue-correction of order A^ + 1, 



note that 



(0) 



HiQo 



0. This is true because the states 



a 



and functionals 



(0) 



with a, P E Pi were chosen to diagonalize Hi in the Sp(^a) subspace, in the sense that 



(0) 



Hi 



a^"*) = A^^^6a,i3- This means that, when we substitute from eq. 13.731 for 



a 



(N) 



in eq. 13.741 the terms that are multiplied from the left with an operator disappear and 
we are left with 



Al^^^^ = {a^^\HiQp,^jHi 



a 



(N) 
R 



N 



E ^'f' 



M=l 



a 



R 



(3.76) 
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which is an expression for the eigenvalue-correction, with all the dead-wood cut away. 



If the an unperturbed eigenstate 



a 



has a non-degenerate eigenvalue, i.e. Hq does 



a 



the situation and hence 



not contain other eigenvectors with the same eigenvalue as 
the formulas are simpler than above. In this case, the proper operator to uses for is 
simply the zero-operator. With this adjustment the general formulas reduce to formulas 
valid for the specific case of no degeneracy. 

Below, we write down explicitly the first order correction to an eigenstate and the 
second order correction to its eigenvalue, for the general case where degeneracy is present 
in the unperturbed system: 



a 



A(2) 



E 
E 



a 



(0) 



p{a) 



E 







(o)^ 






{£-p{a) — £-1^ 









a 



(0) 



a 



(3.77) 



(3.78) 



£p{a) £l 

Note that the formula for the second order correction to the eigenvalue, eq. l3.78l looks the 
same for the degenerate for the non-degenerate case. The only difference is that, 

is a simultaneous eigenvector 



a 



(0) 



in the degenerate case the zero'th order eigenvector 
of Hq and 11^(0,) i^i lying in the £^p(Q,)-eigenspace, whereas in the non-degenerate case it is 
uniquely fixed by diagonalizing Hq. 

We now make a few remarks concerning the results we obtained. Firstly, when we apply 
this perturbation expansion method in conjunction with the boson- fermion mapping that 
we developed earlier, there is still the matter of the physical subspace to address. The 
question is: If we approximate the eigenvectors of H{X) by only taking a few terms in the 
perturbation expansion, do we lose the notion of physical eigenstates that lie exactly in the 
physical subspace? Recall that if we diagonalize H{X) exactly, we can express the physical 
subspace as the space spanned by a certain subset of the eigenvectors of H{X). If we only 
approximate the eigenvectors by the first few terms of their A expansions, can we still find 
a subset of physical approximate eigenvectors that lie exactly in the physical subspace and 
in fact span the physical subspace? Or do the A-expansions for the physical eigenvectors, if 
truncated after a few terms, only lie close to, but not necessarily in the physical subspace? 
Perhaps not unexpectedly, the answer is that the approximate physical eigenvectors still 
lie exactly in the physical subspace. This can be seen as follows. Recall that any operator 
that leaves the physical subspace invariant is said to be a physical operator. Then because 
we assume that H{X) is a physical operator for all A, it follows that Hq and Hi are both 
physical operators. Furthermore, it is a simple matter to show that the operators 11; and 
Qi are physical operators. This means that HiHi is also a physical operator. We find 



the zero'th order eigenvectors a^li] by diagonalizing the operator Ilpfa)Hi in the £. 



■'p{a) 



eigenspace of Hq. Since Iip(^a)Hi is a physical operator with a non-degenerate spectrum. 
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the eigenvectors thus obtained can be separated into a set that spans the overlap of 
physical subspace and the £p(a) eigenspace, and a set of ghost states. Since the operator 
qa is defined in such a way that it has the same eigenvectors in the £p[a) eigenspace as the 
physical operator np(Q,)i^i, it follows that is a physical operator. Furthermore, if we 
start the recursion of eq. I3.fi4l with a state from the physical sector, the corrections that 
the recursion generates remain inside the physical subspace. This follows inductively from 
the fact that higher order corrections in eq. 13.731 are generated by acting with physical 
operators on lower corrections. If the zero'th order eigenvector is physical, this implies 
that the first order correction is physical. If the zero'th order eigenvector and the first 
order correction are both physical, then so too is the second order correction, etc. This 
implies that the exact physical eigenstates, when expanded in terms of A, consist of terms 
that all lie in the physical subspace themselves. Truncating the expansions after a finite 
number of terms therefore still leaves one in the physical subspace. 

ap^l and 



Secondly, as we saw, we know how to express the operators 11;, the states 
as a result also the states c^i}^ in terms of the known states and \iP,Ir). This 

means that all quantities in the results we just derived are expressible in terms of known 
objects. For instance, for the second order corrections to the eigenvalues we have 



A(^)= t ^ {a'?\H,\^M){i^.lL\H,a^^') ^ ^^^^^^ 



Remarkably, it follows from the derivation of the perturbation expansion that it is not 
necessary for the basis to be separable into physical and unphysical states. (The 

basis states enter only as a decomposition of the 11; projection operators, whose properties 
are basis-independent.) Thus we are free to choose the basis in which we express the 11; 
operators such that the calculation of matrix elements of the perturbation is as simple 
as possible, without worrying whether we mix up physical and unphysical states in a 
degenerate subspace. 

We are now ready to start our time-independent perturbation expansion for the 
mapped Hamiltonian. Since all energies but that of the ground state have degenerate 
eigenspaces associated with them, our first task is to find the correct zero'th order eigen- 
states along with the first order corrections to the energies. We take it the original system 
has 2N fermions. Then each degenerate eigenspace of the unperturbed system is labeled 
uniquely by the integer s, and characterized as the space of all states with 2s ideal fermions 
and N — s bosons. In order to find the correct zero'th order eigenstates and the first order 
corrections to the energies of the states in the s-eigenspace, we must diagonalize the op- 
erator UsK in this eigenspace. The operator lis is in this case quite simply the projection 
operator that selects states with 2s ideal fermions and N — s bosons. This makes the 
decomposition K = + Kq + K_ very useful, as already mentioned. More precisely, if 
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|s) is any state in the s-eigenspace of Hq, it holds that 



For the particular case we are interested in, it follows that 



' \s) 


if 


s' + 1 = s 


Ko\s) 


if 


s' = s 


K+\s) 


if 


s' = s + 1 





otherwise 





(3.80) 



II J< \s) 



n-2N ( ^ t 2 ^ . 



fl,U = 



n(N-s) 1 , , 

+ ^ ^2e \s) 



Q - 2s 



(3.81) 



The operator acting on \s) introduces scattering between the ideal fermions of the states in 
the s-eigenspace. The diagonalization of this operator is simplified by the blocking effect 
in much the same way as was the case for the original Richardson Hamiltonian. A simple 
way to consider the blocking effect here is to take it into account before doing the mapping, 
i.e. by restricting the fermion space from which we map to the subspace of full fermion 
Fock space in which no levels are singly occupied. This implies that in the mapped system 
we only work with states oH-Lbf that have no ideal fermions singly occupying levels. To 
verify the above statement, note that the image of aj^al^ creates no ideal fermion in the 
state yU without either also annihilating it or creating an accompanying ideal fermion in 
the state — /x as well. Therefore, because of the blocking effect, the operator in eq. 13.811 
only has to be diagonalised in the space spanned by states of the form 



{N-s)\ Vfc=i 



|0), 



(3.82) 



which is what we now do. 

The operator in eq. 13.811 treats the bosons as spectators. Furthermore, the term 
proportional to identity and the scaling factor ^^z^ in the first term are unimportant for 
the diagonalization procedure. This allows us to construct the eigenstates of the operator 
in eq. 13.811 as follows: First we find every eigenstate l-^) (with eigenvalue r}^) of the 
operator 

-f^O.eff 



Vt-2s + 2 



(3.83) 



in the space of states with 2s ideal fermions with no singly occupied levels and no 



bosons. The eigenstates of the operator in eq. 13.811 are then 



with eigenval- 
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ues 



n-2N 



N-s 



2Qe. The real problem to solve is diagonalizing the operator -fCo.efr in 



n-2s 'IV ' Q.-2s' 

eq. 13.831 in the space of 2s ideal fermions with no singly occupied levels and no bosons. 
This last operator is sufficiently similar to the Richardson Hamiltonian that we started 
with that its diagonalization proceeds with the same ease or difficulty as the diagonaliza- 
tion of the Richardson Hamiltonian for 2s electrons. What we gain, for low lying states 
at least, is that the number of ideal fermions (2s) for which the operator Kq ^^ in eq. 13.831 
must be diagonalized, is much smaller than the number of electrons (2A^) for which the 
Richardson Hamiltonian must be diagonalized. 

After playing around with the operator Xo^ofr, and for instance diagonalizing it in the 
space of a two-particle system, one starts to suspect that it has a lot in common with the 
original Richardson Hamiltonian. We therefore try to diagonalize it following the same 
procedure that was followed in [211 Appendix B]. This involves making an ansatz for the 
eigenstates of the operator Kcefr- We assume that every eigenstate I?/') has the form 



IKJIo)> 



(3.84) 



\k=l 



where J^^^p, = Z^^^i 2e -r?^ ^ f^i^ complex numbers rj^^k still have to be fixed. Because 

the ansatz for the form of the eigenstates uses the f3l operators, it is useful to collect a 
few facts about these operators, namely 

= 0, 



(3.85) 
(3.86) 
(3.87) 



and to rewrite the operator in eq. 13.831 in terms of the /?'s as 

2 



. 1 



(3.88) 



where we used the fact that e^, = e_ 
First we compute the commutator 



Q - 2s + 2'' 

With this form of the operator we check the ansatz. 



0,eS 



A. 



ip,k 



E 



I Vtl>,k 



\ 



Pp 1 /^f 1 Pp 



n-2s + 2 



^2e 

P— 2 V 



2er. 



Pi 



1 



n-2s + 2 



E 



2er. 



1 ^t/',fc 



=<5..p/3+[/3p,/3j]/ 



Pi \Pl , P^ 



- 2s + 2 



E 



2e. 



Pi , P, 
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(3.89) 



where, in the last hne, we used the relation ^/nS'^ = J2'' -i Pi- Now we compute 

p— y y 



0,eff 



Q , n 4,^ 



k=l 



|0) 



k=l \l=l 



Kl=k+1 



k=l l=l^k 



n-2s + 2 



k=i \i=i 



Pi , Pp 



n <i |0) , 

j=k+i J 

(3.90) 



where t]^ = J2k=i V^^k- In order to proceed, we note that 



piPp] n </io) 

l=k+l 



pi^pp ' n 



l=k+l 

/-I 

n 

l=k+l \m=k+l 



E ( n ^\,m 



io)+ n </[/?j,/3p]io) 

=-|o) 

Pl,Pp I n ) |0) 

\m=l+l 



- n </io)- 

l=k+l 

This last result we substitute into eq. l3.9Ul to find 



(3.91) 



0,eff 



k=l 



I 



V 



fe = l ( = fe + l \m=l7i:fe 



^264[/5t,/3, 
X L - 



2ep - 



n ^i,™ I |0) . (3.92) 



The trick now is to deal correctly with the indicated commutator. 



X 2e„ 



E 



P=2 



1 2ep-r7^, 



E 



2e. 



t (2ep - ri^^k) (2ep - r^v>,«) 



Vi>,k — Vtp,i 



(3.93) 
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where the second hne follows from the first if we notice that 



2e. 



1 



— I . (3.94) 

The last manipulation is not valid when both rj^^k and rj^^i are zero. We side-step this 
problem by assuming that none of the ?7^,fe, k = 1,2, ... ,s are zero, and later treat the 
case where some are. (The case where only one of the ?7's is zero could have been included 
in the present analysis, but fits more naturally with the treatment for the case where any 
number of the t^'s are zero.) By substituting eq. 13.931 into eq. 13.921 we find 



O.eff 



E 



2e. 



k=l 



n </io) 



Vip,k 



Q-2S + 2 
'Q-2s + 2 



^1 i^^i Vij,k Vi>,i 

E/ E/ ^'01^ 
k=i i=k+i ^.'^ ~ 



n Am 1 10) 



|0) 



(3.95) 



We now enumerate differently the first k, I summation by noting that for the set X of 
indices over which the summation runs, it holds that 



J, . . . , ^, 



s; 1 = k + l,k + 2. 



I = {{k,l) : k = 1,2, 

= {{k,l) : I = 1,2,..., s; k = l,2,...,l-l}. 

and rename some dummy-indices to find 

Ko,eS lip) = % 1^) 

1 



s} 



(3.96) 



-Vns^ E { 1 



k=l 



Q-2S + 2 



E 



2e„ 



1 2ep - 77^, 



+ 2E 



,_l V^,l ~ V^,k 

i o 



X n <n^|0)• 

m=ljtk 

(3.97) 



If we choose the rj^^k to satisfy 



1 



n-2s + 2 



E 



2e, 



2 E 



Vi>,i 



(3.98) 



for A; = 1, 2, . . . , s we have -ft'cefr 
can be rewritten as 

' 1 



= r]^ 1-0). Provided none of the rj^^k are zero^, this 



o = E 



+ 2 E 



k — 1, 2, . . . , s, 



(3.99) 



"^The derivation involves dividing the equation enumerated by k by ??^,/c. 
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which are just Richardson's equations in the limit where the coupling constant dominates 
the one-particle energies — > 0). Now we return to the case where one or more of the 
r] are assumed zero. Define an operator Aq = Z^^^i ^Pl^ which is just Alp^^, with r]^^^ set 
to zero. The commutator between Kq^cS and Aq reduces from that given in eq. 13.891 to 

[i^o.eff ,Ao] = v^^t |l - ^/^^^} • (3.100) 
If 1 2s — 2) is any state with 2s — 2 ideal fermions, then |l — S7^^7+ 9 } 1 2s — 2) = so that 

R\csAl |2s - 2) = AlKo,e{i |2s - 2) . (3.101) 

Therefore, the eigenstates of i^o,efr in the space with 2s — 2 ideal fermions, multiplied by 
Aq are the eigenstates of -ft'ces with 2s ideal fermions and at least one of the rj^^k zero. 
The implication is that r]^^k = solutions do exist and should be treated as follows: The 
solutions with t < s of the rj zero are obtained by finding the eigenstates of Xo,efr with 
2(s — t) fermions such that none of the rj^^k, k = 1, 2, . . . , s — t are zero. This is done by 
solving eq. 13.991 with s replaced by s — t. The corresponding eigenstates of i^o.efr with 2s 
fermions and t of the r] zero, are 

l^) = (4)*n<JO), (3.102) 

k=l 

with eigenvalues r]^ = J2k^i V^,k- 

We are now ready to write down the first order corrections to the energies, and the 
correct zero'th order eigenstates that the perturbation selects, for each degenerate s- 
eigenspace of Hq. For given s and t we enumerate the different solutions to 

= E ^ + 2 E ^ ' A; = l,2,...,s-i (3.103) 

with different ip. Corresponding to each ip is an eigenstate |'0,s,t) of Hq given by 

N-s 



1^, s, t) = {Aiy n t ^-4^ 1 10) • (3.104) 



These are the correct zero'th order eigenstates selected by the perturbation. To zero'th 
order the state 1-0, s,t) has an energy 

= -^(iV-s)(fi + l-iV-s). (3.105) 



To this is added a first order correction due to the perturbation 

^* - = TT^iJ E V*, HN-s) —21. (3.106) 

fe=l 
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From the eigenstates in eq. the subset that spans the physical subspace must now 

be selected. The procedure we employ to do this, is improvised specifically for this model. 
We shift all the one-particle energies by the same amount e. On the level of the original 
Richardson Hamiltonian, this amounts to adding a multiple of the total electron number 
operator, i.e. H — > H + eX]v=-j ^^ctv- This means that the energies of the 2N electron 
eigenstates of the Richardson Hamiltonian are shifted by an amount 2Ne. In the mapped 
system, only the physical eigenstates have anything to do with the original Richardson 
Hamiltonian, and so the eigenvalues of physical eigenvectors must still shift by an amount 
2Ne. States that do not lie in the physical subspace, do not have to behave this way. 
(Although, there is nothing that says they cannot shift their eigenvalues by an amount 
2Ne). Therefore, if we find an eigenstate of the mapped system whose eigenvalue does 
not shift by an amount 2Ne when each is shifted by an amount e, we know it is a ghost 
state. Let us therefore calculate the shift in the first order correction A^^l, j to the energy 
of the state |'0,s,t) when we replace each with e'^, = + e. The average one-particle 
energy shifts from e to e + e. The argument about the 77^^^ that are zero does not involve 
the values of e^. Thus, for given t, there are still solutions with only s — toi the rj non-zero 
or in other words the rj that were zero before the shift remains zero after the shift. We 
we will denote the new non-zero i] hj t]'^ f., k = 1, 2, . . . , s — t. They must satisfy eq. 
IH.IOHI with each replaced with Cp + e. It is trivial to see that if rj^^k satisfied eq. IH.IOHI 
before we made the shift, then ^7^^ = ?7^,a: + 2e satisfies it after the shift is made. Hence, 
the non-zero 77 all get shifted by 2e. This means that the the first-order correction to the 
energy of state 1-0, s,t) is 

= 7)3^S(')« + 2.) + (A'-.)jj3^2(.-+.) 

k=l 

= A^^l^ + 2Ne-2t^^^^. (3.107) 

To get a shift of 2Ne we have to set t = 0, which means that the states \ilj,s,t) with t 7^ 
are ghost states. 

The question now is whether these are all the ghost states. The answer is affirmative. 
To show this we count the number of eigenstates we are left with after removing the t ^ 
states, and compare this number to the dimension of the physical subspace which we know 
to be (J^^ . Let Us therefore be the number oi t = states in the s-eigenspace. According 
to eq. 13. 9 9^ these states are associated with finite solutions of Richardson's equations 
for 2s ideal fermions in the limit of infinite coupling. If we accept the result that was 
stated without proof in Section ESI that for given s there are x such solutions 

and therefore as many t = linearly independent eigenstates in the s-eigenspace, we can 
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compute the combined number of eigenstates left over in all the s-eigenspaces as follows: 

N 
s=0 




(3.108) 

This result says that the number of t = states just makes up the dimension of the 
physical subspace, which confirms that the t = states span the physical subspace. 

At this point, let us collect the main results together (for the benefit of readers who 
skipped over the preceding intricate arguments) and briefly discuss their implications. 
We exclude ghost-states from our discussion. The unnormalized zero'th order physical 
eigenstates of the mapped Hamiltonian, selected by the perturbation K, are 



St 



N-s 



s 



( 



{N-s) 



k=l 



Pi 



where the parameters r]^,k, k = 1, . . . , s refer to one of the (^j^^itM+i) finite solutions of 
the strong coupling limit Richardson equations 

0=E ^ + 2 E ' k = l,2,...,s, (3.110) 

'rii>,k — i=i^k ^'Z'.' ~ ''^V'.fc 



p= 



and different ■0 refer to different solutions. The quantum number s takes on values from 
for the ground state to for the most excited states. The energy of the state \ip,s) 
with the first order correction included is 

G r2 — 2N ^ 

E^,s = --{N-s){n + l-N-s) + {N-s) ^^3^26 + -^^—^ E V^,k ■ (3.111) 

. ' ^ , ' , fe=i 

(1) (2) 5^ 

The results for the non-degenerate ground state is obtained by setting s = in the 
the formulas for the eigenvector, eq. 13.1091 and energy, eq. 13.1111 and dropping all 
products and summations where the lower bound exceeds the upper bound. This is 
the same result as was found in jSHj. The way we arrived at it though, suggests an 
amusing new interpretation not discussed in [SH] , for the low-lying spectrum to this order 
in perturbation theory. We assume the number of one-particle levels Vt to be large, so 
that, for low-lying states, terms of order ^ may be neglected. In that case term (1) 
represents an equidistant spectrum with level-spacing G. It counts the number of bosons 
in the state \ip,s), assigning an energy —j:N{yL — A -|- 1) if there are A^ bosons, and 
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approximately increases the energy by an amount G for each boson less than A^. Thus, 
term (1) approximately constitutes an additive 'pairing' energy decrease of G per boson. 
Term (2) can be seen as the total one-particle energy contribution of the pairs of electrons 
that went into making the 'Cooper pair bosons', loosely speaking. This interpretation 
follows from the fact that the average one-particle energy of a single Cooper pair is 
(0| S'I]^=_j EyO^ayS'^ |0) = 2e together with the fact that the state 1-0, s) contains N — s 
bosons. So term (2) represents an additive 'internal' energy of 2e per boson. (The factor 
is close to one for small s.) For given s, term (3) represents the energy contribution 
of 2s ideal fermions, interacting through the Richardson Hamiltonian, where the energy 
scale is scaled by a factor of approximately 1 — ^ and the coupling constant G is sent 
to infinity while insisting that only finite energy eigenstates are accessible to the ideal 
fermions. This interpretation of term (3) is corroborated by the form of the eigenstate 
in eq.EHSl 

The fact that the ideal fermions are arranged in eigenstates of the original (hermitian) 
Richardson Hamiltonian with a strong pairing interaction, implies that they are mutually 
orthogonal. This is an amazing result: Because the unperturbed mapped Hamiltonian 
is hermitian, the various s-eigenspaces are orthogonal to each other. However, since 
the perturbation K and more specifically the component Kq that determines the zero'th 
order eigenstates, is non-hermitian, the zero'th order eigenbasis for an s-eigenspace that 
it selects, will not be orthogonal. Yet, for the mapped Richardson Hamiltonian the non- 
orthogonality disappears when the ghost states are removed. In other words, only the 
ghost states break the orthogonality of the basis selected by K, so that the orthogonality 
of the zero'th order physical eigenstates is preserved by the mapping. This makes the task 
of finding the left basis that goes along with the states |'0,s) simple. It should however 
be noted that these remarks do not imply that the boson-fermion mapping is unitary. 
Although the orthogonality of the states \il)iS) are preserved, the mapping still introduces 
stretching of the norms of these states. 

To proceed to higher orders in the perturbation expansion now 'only' involves the 
calculation of matrix elements of the perturbation K and the summation of terms con- 
taining products of these matrix elements. The calculations involved soon become suf- 
ficiently inelegant to convince one that one is only interested in the first order energy 
corrections and that the correct zero'th order eigenstates will do just fine. Yet, should 
one proceed to higher orders in perturbation theory by calculating the necessary matrix 
elements, the least unpleasant route to follow is, I believe, to work with the mapped 
system. Then one can express the H; operators in terms of free particle basis states of 

the form ^, nfc=i Pi, |0), for which the action of the perturbation is not too hard 

to calculate. Honesty compels me to mention at this point that the method of [21] has 
an advantage over the above procedure, namely that the calculation of matrix elements 
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is completely circumvented. However, the method is designed only for the Richardson 
Model. There might thus be models for which the only way to proceed is conventional 
perturbation theory, as we have done, in which case the Dyson mapping might be a useful 
tool. 

In the work we have done so far we have not yet used the full image of the Richardson 
Hamiltonian. We have only used the expressions for Hq and Kq. The expressions for 
and K_ are as yet untested. Therefore I calculate the second order correction to the 
ground state energy. This calculation is not too complicated and uses the full image of 
the Richardson Hamiltonian. 

To zero'th order, the ground state is |GS) = ^= [B'^\ |0). According to eq. 13. 781 the 



second order correction to the ground state energy is 

A (2) _ \^ (^^1 -^in^^-^l |GS) 

s=0 -C^o " 

where H^ is the projection operator that selects states with 2s ideal fermions and N — s 
bosons, i.e. states in the s-eigenspace of i/o, and Hi is the perturbation K. Again 
the decomposition K = + Kq + K_ is very useful. Because Kq leaves every s- 
eigenspace invariant, while K± takes states from the s eigenspace to the s ib 1 eigenspace, 
the expression for the second order correction to the ground state energy reduces to 

^(,)_ (GS|i^-H,=iK+|GS) 



We conveniently choose the orthonormal basis \u) = -^-^=^=pj^ |0) for the (s = 1)- 

eigenspace, so that the operator Ils=i can be expressed as H^^i = -i I/") (/^l- If we 
substitute this into the expression for the second order correction to the ground state 
energy, we see that we have managed to express the correction in terms of the matrix 
elements (GS| \^) and (//| X+ |GS). These matrix elements are not hard to calculate, 
the result being 

(GSIE^.AVI/^) = 26,^, 

ME^.^'.+ IGS) = 2(6.-6)^^^. (3.114) 

The energy difference in the denominator is Eq^^ — e[^^ = —G. Combining everything we 
find a second order correction to the ground state energy 
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/ 1 '^u 

where var(e) = — e^. Note here how we used a basis for the (s = 1)- 

eigenspace that mixes physical and ghost state. Again, the result agrees with the result 
obtained by (22] . (The method of the authors of [SH] allowed them to calculate this 
particular correction to the seventh order.) 



3.6 Time-Dependent Perturbation Theory 

The good news about the results of eq. 13.1111 and eq. 13.1151 obtained through time- 
independent perturbation theory is that they agree with the results obtained by 
Yuzbashyan and co-workers in [22]. In a sense this is, however, also the bad news: results 
already exist to compare with everything we calculated. As mentioned in the introduc- 
tion, the time-independent perturbation expansion that uses the boson-fermion mapping 
does not yield anything that is not in the scope of previous analyses. In an attempt 
to calculate something with the aid of the boson-fermion mapping that would otherwise 
be hard to calculate, we therefore consider a time-dependent perturbation expansion for 
transition amplitudes between states of the unperturbed system. The reason I believe 
this calculation to be simplified by the mapping will be highlighted at the appropriate 
point in the derivation of the expansion. 

We start by reviewing some basics of time-dependent perturbation theory. We have a 
quantum system with a Hamiltonian H = Hq+Hi where both Hq and Hi are hermitian. It 
is assumed that the orthonormal eigenstates of Hq are known. We denote these eigenstates 
la), with Ho\a) = Ea\a). Here a stands for all the quantum numbers necessary to specify 
a state uniquely, and it is immaterial whether the spectrum is degenerate. We assume 
that Hi is a small perturbation to Hq. For this system let U(t) be the time-evolution 
operator which satisfies ih^^^j^ = HJJ(t), with the boundary condition (j{0) = I. We 
define an interaction picture time-evolution operator 

Ui{t)=exp(^--^^U{t), (3.116) 

which is close to identity when Hi is small. Hence we have the hope that it is possible 
to expand Ui{t) in powers of Hi around identity. To find this expansion we go about as 
follows. First we note that Ui{t) obeys a differential equation similar to that obeyed by 
U{t), namely 

in^^^ = Hu{t)Uj{t), (3.117) 
at 

if we define Hu = exp (i^^^ Hi exp (^—i^^. This differential equation can be converted 
into an integral equation, 

Uj{t) = 1-^1^ dt' Hiiit') + [-'^ ' dt' dt" Hij{t')Hij{t") + ..., (3.118) 
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taking into account the boundary condition Ui{t = 0) = I. Our goal is to calculate 
the transition amplitude at time t between eigenstates of Hq, i.e. we are interested in 
a quantity U{t) \a). Using the definition of Ui{t) we can rewrite this in terms of a 
matrix element of Ui(t) as follows: 

{p\ U{t) \a) = exp (-«^) Ui{t) \a) . (3.119) 

The task that remains is to calculate the matrix element Uj{t) \a) to some order in 
Hij with the aid of the power expansion for Uj{t) in eq. 13.1181 This involves calculating 
matrix elements of operators like Hij{t')Hij{t") with respect to eigenstates of Hq. 

Let us suppose at this point that the Hamiltonian H that we started with is a fermion 
Hamiltonian for which a boson-fermion mapping is available, and let us further suppose 
that it is a simpler task to compute the above-mentioned matrix elements for the mapped 
system. In this section, we use the following notation: For the similarity transformed 
operators we will use the notation Q = TQT~^. We have here kept quiet about the 
fact that we have at our disposal extended images of operators and functionals. Rather 
untidily we will from now on understand i^o to denote the operator defined on the whole 
boson-fermion Fock space which, when restricted to the physical subspace, is equivalent 
to Hq. We can choose as bi-orthonormal basis for the whole boson-fermion Fock space, 
the left and right eigenstates of this operator, properly normalized: 

Hq M = Ey \vr) , 

il^L Wr) = S^,. (3.120) 

The luji) and |z/l) bases contain states outside the physical subspace, and thanks to the 
possible degeneracy of physical and ghost states may even mix physical and ghost states 
together. None the less, the identity in the whole boson-fermion Fock space, in terms of 
these states, is 

I = El'^i?)KI- (3-121) 
In terms of our mapped operators and states, we find to second order in H^i 
{f3\lJj{t)\a) = 5^p-'-f^dt' exp (^ iEp-^.)t' ^ {(3l\T~'H,T \an) 

+ ..., (3.122) 

We now apply these formulas to the mapped equivalent {H)j^ = K+Hq of the Richard- 
son Hamiltonian that we derived in Section 13.41 As usual, we take the unperturbed 
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Hamiltonian to be Hq that only includes the pairing effect, and treat the one-particle 
energy contribution, as the perturbation. We will calculate the ground state to ground 
state transition matrix element to second order in the perturbation K. If the normalized 
ground state of the original system is |GS) then we know that for the mapped ground 
state we have 



IGS 



R 



r|GS) 

C r— 



|o), 



(3.123) 



where c is a (real) number which we will not need to calculate. Similarly we have 



(GS, 



(GS|T-i 
(0| 



(3.124) 



The next step is to choose a convenient eigenbasis of the unperturbed Hamiltonian in 
which to express the identity operator. Because we only work to second order in the 
perturbation, and because of the frequently mentioned property that the perturbation 
takes states in the s-eigenspace to states with components only in the (s — 1)-, s- and 
(s + l)-eigenspaces, we need to specify this basis only up to the (s = l)-eigenspace for 
the ground state (s = 0) to ground state transition amplitude. We choose to work with 



the right eigenbasis consisting of |GSij) and |z/) 

corresponding left eigenbasis (GSi| and (z^| = (0| /5y 
choice eq. 13. 1221 reduces to 
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(3.125) 



Since, in every term where \GSji) appears, there is a corresponding {GSl\, the unknown 
proportionality factor c is always canceled by its inverse c~^. Forgetting about the c's then, 
we have already calculated all the matrix elements appearing in the above expression. It is 
here that the advantage of using the boson-fermion mapping lies. The states with respect 
to which the matrix elements are calculated for the mapped system are free particle 
states. This is to be contrasted with the strongly correlated eigenstates of Hq = —GS^S 
we would have had to use in the unmapped system. In the case of time-independent 
perturbation expansions Yuzbashyan got around the problem of calculating these matrix 
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elements thanks to the Richardson equations. However, for time- dependent perturbation 
expansions there is, as far as I can see, no similar trick to avoid the calculation of the 
matrix elements. This makes the mapping very useful. 
Evaluation of the double integral gives 

/ Gt'\ ( Gt"\ nt / -itG \ 

Substituting this, and the previously calculated expressions for the matrix elements of the 
perturbation into eq. 13.1251 we find 
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(3.127) 



For this truncation to be valid two conditions have to hold. Firstly, as with the time- 
independent perturbation expansion, the coupling constant G should dominate the one- 
particle energies. Secondly, the expansion is only valid for small enough times: the one- 
particle energies must be much smaller than j. We can use eq. 13.1271 to calculate the 
probability of finding the system in the unperturbed ground state at a time t -C if it 
was in the unperturbed ground state at time t = 0. The result is 

|(GS| U{t) \GS)f = 1 - 16iv|^^sin^ |. (3.128) 

The same result could also have been obtained by calculating the total amplitude for 
scattering to excited states. We do not perform the calculation but merely mention the 
following. The first non-zero term in the expansion for the amplitude for scattering to 
excited states is of the order ^. To get the total probability to order ^ for excitation to 
any state above the ground state, only the first (order ^) non-zero term in the excitation 
amplitude is therefore required. To this order, only excitation to the s = 1 level is possible. 
I contend that such an excitation represents the decay of one Cooper pair. The argument 
is as follows: We map back to the fermion system to find that the states in the first excited 
level (up to normalization) are given by 



/ 



aZa: 



E ^ _; I |0) , (3.129) 



where r]^ is a finite solution to = Z)^_i 27^^' (Actually, we map the state in eq. 
13.1291 in the forward direction and confirm that it is equal to one of the zero'th or- 
der eigenstates in the s = 1 level of the mapped Hamiltonian that is selected by the 
perturbation.) In such a state, the collective pair operator 5"+ creates a Cooper pair 
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with 'long-range' momentum space correlations, in the sense that all the pair opera- 
tors {—Y'^a'^atiy contribute with equal amplitude. It is this type of correlation which 
defines a Cooper pair. (The momentum we are talking about here is the momentum 
difference between the two particles in the pair. The total momentum of the pair is 
zero. Also, we assume that momentum correlations closely mimic Bloch wave-number 
correlations.) In contrast, if we set {u) = {—y^^a^a^^ |0) and examine the amplitude 
f4i{v) = (/^l S^-i — 2e -r]^ ~^ I*-*) ^ function of u, we find that it is sharply peaked. In 

this sense, the operator — 2e ~^ creates a two-electron state with only short- 
range momentum correlations, and which is well-localized in momentum space (or more 
correctly, Bloch wave-number space). Thus, the two particle state created by such an 
operator is 'close' to an eigenstate of the two-particle system with the pairing switched off 
and far from the type of correlation that marks a Cooper pair. A state in the s = 1 level, 
such as that in eq. IH.129| therefore contains — 1 Cooper pairs. The two remaining elec- 
trons are, to a good approximation, in eigenstates of the single particle energy operator 

Note that without doing a detailed calculation we could have known the probability to 
remain in the ground state is of the form 1 — c sin^ This form follows from the fact that, 
to second order, the only processes possible are (a) the system remains in its ground state 
and (b) a single Cooper pair decays and recombines, the energy involved in this virtual 
transition being G. All the detailed information about the non-zero matrix elements of 
the perturbation resides in the coefficient c, which in this case is equal to 16N^^^^^^^^. 
The fact that the energy dependence takes the form of a variance, guarantees that shifting 
all one-particle energies by the same amount does not infiuence the transition probability. 
Furthermore, for large Q, the expression grows linearly with the number of particles. 
In other words, 1 — |(GS| U{t) |GS)|^, which is the probability for one Cooper pair to 
have decayed by time t, doubles if we double the number of particles. This makes good 
sense. In our intuitive picture. Cooper pairs are independent entities that, to a good 
approximation, propagate without noticing each other. They have finite lifetimes because 
of the perturbation K. In this picture, the probability for one Cooper pair out of N 
to have decayed after a (short) time t, is extensive, as it is the sum of the independent 
probabilities for Cooper pairs 1, 2, ... or A^ to have decayed by time t. For smaller Q this 
independenc is clearly inhibited by the factor 

3.7 Conclusion 

Much of the analysis presented in this chapter is new, and contains contributions relevant 
to both the fields of generalised Dyson mappings, and the Richardson Hamiltonian. As 
far as the generalised Dyson mapping is concerned, we analysed what is perhaps the 
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simplest system for which the Hamiltonian cannot simply be written in terms of only 
linear combinations of products of operators in the even sector of the superalgebra. We 
also included states with more than one ideal fermion in the mapped system. To my 
knowledge, no boson expansion method type analyses has previously been applied in such 
a general setting. It was pointed out that, when there are Q ideal fermions in the mapped 
system, the similarity transform of Section ITTTI breaks down. We demonstrated how to do 
perturbation theory, both time-dependent and time- independent, for a mapped system. 
The ease with which a perturbation expansion can be made in the mapped system provides 
a motivation for doing a generalised Dyson mapping. 

As far as new insights into the Richardson model goes, it was demonstrated why 
doing a strong coupling expansion is feasible. The key insight here is that the one-body 
term ^uO-to-y only introduces scattering between adjacent energy levels of the infinite 
pairing system. This result could in principle be derived without doing a mapping, but 
is immediately obvious only in the mapped system. 

I conclude with a speculation on possible future developments involving the mapped 
version of the Richardson Hamiltonian. In presenting the results of time-dependent per- 
turbation theory, I have stressed that the aim is not to understand measurable effects in a 
superconductor. The reasons for having to make this cautionary remark are the following. 
Firstly, the strong coupling limit is not realised in a superconductor. Secondly, we have 
little control over the state the electrons are in at the start of the experiment. It is not 
possible to prepare the system in the ground state of the infinite pairing regime and then 
to let it evolve with a finite (large) pairing force. 

A setup much more suited to the kind of calculation we performed here, is that of ultra- 
cold Fermi gases |22I- There atoms that are believed to obey Fermi statistics are cooled 
in a trap, to a point where they condense into what is loosely called a single macroscopic 
wave-function. By means of the clever use of external magnetic fields, a situation can be 
created in which the interaction between the fermionic atoms can be adjusted. Depending 
on the type of interaction the magnetic field is tuned for, the macroscopic wave-function is 
either a condensate of Cooper pairs, or a Bose-Einstein condensate of diatomic molecules. 
On a completely practical level, it is good news that the pairing strength of the Cooper 
pairs can be adjusted. Therefore a system can be prepared in the infinite pairing regime. 
In a short time-interval the pairing force may then be decreased, and the resulting state 
of the system probed. An analysis of such a situation, that goes beyond the mean field 
theory of jSS] , is possible with the techniques developed in the section on time- dependent 
perturbation theory. The two challenges to this program are the computation of high 
order corrections and the identification of observable quantities that indicate the change 
the system undergoes. For instance, with reference to the result in eq. 13.1281 how might 
it be possible to see if a single Cooper pair has decayed? 
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On a more conceptual level, the whole question of when a system is a Bose- Einstein 
condensate of diatomic molecules and when it is a condensate of Cooper pairs, and what 
happens in between, may supply fertile ground for generalised Dyson boson-fermion map- 
pings. Since a textbook system representing a Bose-Einstein condensate has many differ- 
ent one-particle levels that the bosons can occupy, while the mapping we considered in 
this chapter piles all the bosons into the same level, more work is required for this analysis. 
A mapping needs to be devised that allows more one-particle levels for the bosons. The 
big stumbling block in this enterprise is that the similarity transformation of Section 11.71 
will then become much harder to deal with. 
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